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Some Phenomena of the Upper Atmosphere 


By SYDNEY CHAPMAN * 
Queen’s College, Oxford 


Presidential Address, delivered 2nd May 1951; MS. received 21st May 1951 


ABSTRACT. The height distribution of the atmosphere depends on the temperature and 
the mean molecular weight of the air; recent progress in our knowledge of these two 
quantities will be reviewed, the main sources of information being outlined. The motions 
(winds and tides) in the upper atmosphere will similarly be considered. 

The main gross uncertainty regarding composition at high levels relates to the proportion 
-of atomic nitrogen present. Besides the major constituents that chiefly determine the mean 
molecular weight, there are important rare constituents, such as ozone, sodium and hydroxyl, 
which reveal their presence by their absorption or emission spectra; and there are others 
whose presence can be inferred from theirs. 

The composition of the atmosphere changes by escape from the top (hydrogen and 
helium) and, by addition, from below (e.g. helium and carbon dioxide) and, as shown by the 
auroral spectrum, from above (hydrogen); cosmic rays also both add to the atmosphere and 
produce changes of composition by nuclear reactions. 

Near the magnetic equator there are important ionospheric phenomena not yet fully 
explored and explained. Among these one of the most interesting is the abnormal intensifi- 
cation of electric current flow, specially notable over Huancayo in Peru. 


§1. INTRODUCTION 

VEN a slight acquaintance with the history of thought suggests that the 

human mind generally thinks too simply of matters on which it has little 

knowledge; the phrase omne ignotum pro simplict might often describe 
the position. ‘This is well illustrated in the science of the atmosphere. After 
Torricelli had invented the barometer and Pascal had shown that the barometric 
pressure decreases with height according to the weight of the intervening air, 
even so great a man as Halley was ready to think that the atmosphere was sub- 
stantially understood. In the eighteenth and nineteenth centuries balloonists 
found that in the free air, as well as on mountains, the air temperature decreases 
with height. ‘The decrease was associated with the strong convection present 
in the air at low levels, being about half the adiabatic rate characteristic of thorough 
mixing. It was generally tacitly assumed that this decrease continued, so that 
the atmosphere would have a finite height. 


* Professor Chapman was elected President for the period 1949-51, but owing to his absence 
in America during 1950-51 he was only able to hold this office for one year, and the Presidential 


Address had to be postponed. 
PROC. PHYS. SOC. LXIV, 10—B ~ : Bur 
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But early this century Teisserenc de Bort * found by kite measurements 
that above about 10km. in our latitudes the upward decrease of temperature 
ceases, and may be followed by a slight upward increase. The layers below 
and above this level of change of temperature gradient were called the troposphere 
and stratosphere, and the level itself was called the tropopause. : 

When de Bort’s discovery gained general acceptance, opinion swung round. 
to the conception of a static stratosphere, in which the temperature was controlled | 
mainly by radiative processes, and remained constant at all greater heights. 
This conception led also to the inference that the gases of the atmosphere settled 
out by relative diffusion, so that the proportion of the heavier constituents; 
decreased steadily upwards, and ultimately the lightest gases predominated. 
Calculations on this basis were made as to the height distribution of the pressure : 
aud density in the air, both as a whole and for each constituent separately; but} 
by 1920 or so there were questioning voices, as regards the constancy of! 
temperature, and also as to the diffusive separation. 

There was indeed already ample evidence that in the atmosphere much goes; 
on that is not comprised in so simple a philosophy. 


§2. THE COMPLEXITY OF THE UPPERVATMOSPHERE 


Even in the eighteen-eighties Balfour Stewart had inferred that somewhere: 
high in the atmosphere the air is rendered electrically conducting by solar radiation, , 
and that electrical currents flow there, which are induced by the dynamo action} 
of the motion of the air across the earth’s magnetic field—whose variations : 
observed at the earth’s surface formed the basis of these remarkable inferences. 

Before the first World War Marconi’s success in radio telegraphy over great { 
distances led Heaviside and Kennelly to reach anew the conclusion that the upper | 
atmosphere is ionized; and in the nineteen-twenties Appleton, and Breit and} 
‘Tuve, by new methods of radio research, wonderfully confirmed these inferences, / 
and added numerical substance to them. They and their successors have since? 
contributed a wealth of complicated facts about the upper atmosphere, in the? 
region named by Watson-Watt the ionosphere—facts gained by that mostt 
penetrating of all our present means of atmospheric exploration, the radio: 
transmitter and receiver. 

But Nature unaided also shows that the upper atmosphere is not simple. | 
Meteors continually shoot into it and appear as transient falling stars; a few of] 
them leave trails that endure for minutes, or even for an hour or more, in rare4 
cases. ‘These trails almost always become contorted as they fade, and thus. 
indicate clearly the presence of non-uniform winds, at levels of 70 to 100km. 
and more. 

During the first World War the sound of gunfire in Flanders was heard here 
in London, and since then many experiments have been made with artificial] 
explosions to investigate such abnormal propagation of sound to zones beyondj 
the distance at which the sound usually ceases to be heard. It became clear that} 
the abnormal sound came down from above, and indicated a downward deflection! 
of waves originally travelling upwards. This implies the presence of a region, 
at a height estimated at about 40 km., in which sound travels more rapidly than i 


* At the time of printing, Professor Chapman’s papers had not arrived from America; it was} 
therefore impossible to insert bibliographic references in the text in the normal manner./ 
These references, however, are listed at the end of the article. 
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the air at ground level. One fantastic interpretation of this was that the air in this 
region consists of hydrogen; but the most complete quiescence in the stratosphere 
could not enable diffusion to render hydrogen predominant at that level. 

The now accepted alternative explanation is that the air in the said region is 
hotter than near the ground. But it was on the basis of a theory of meteors and 
their luminosity that Lindemann and Dobson first definitely concluded that 
there is an upward increase of temperature above de Bort’s isothermal layer. 

For a higher level, at about 100km., near the lower limit of aurorae, Vegard 
found evidence in the auroral band spectrum for a much lower temperature, of 
the same order asatthetropopause. Fora time this seemed doubtfully compatible 
with the ideas of Lindemann and Dobson, but support for both views came during 
the nineteen-thirties from an unexpected source—the study of the daily variations 
of the barometer. 


$3. ATMOSPHERIC TIDES AND THE TEMPERATURE DISTRIBUTION 

Laplace, who first developed Newton’s tidal theory in its application to the 
atmosphere, was also the first to seek to determine the lunar atmospheric tide 
from the barometric observations. After more than a century of effort in this 
direction, the distribution of the lunar air tide over the earth is now fairly well 
ascertained. Substantially it is in phase with the moon, that is, the air is heaped 
up twice daily, raising the barometric column, when the moon is on or near the 
meridian of observation or the opposite meridian. But the whole range of 
pressure thus produced is very small—only about a tenth of a millimetre of mercury 
—and, whereas the sea-tides are controlled by the moon and only modified by 
the sun, the moon makes only a minor contribution to the air tide. The chief 
daily barometric variation is solar semi-diurnal, with maxima at about 10 a.m. 
‘and 10 p.m., and minima at 4a.m. and4p.m.; this variation has a range of about 
2mm. at the equator, about fifteen times as great as the lunar air tide. 

The reason for this peculiar reversal of tidal predominance as between the 
sun and the moon, and the seas and the atmosphere, was first suggested by Kelvin 
in 1882, but it was not substantially established until half a century later. G. I. 
Taylor and C. L. Pekeris then confirmed Kelvin’s idea by showing that the 
atmosphere can have more than one mode of free travelling wave of tidal type, 
and that one of these can reasonably be supposed to have a period of almost 
exactly 12 solar hours. If this is so, the solar tide—to which the sun’s thermal 
action contributes about equally with its gravitational attraction—can be magnified 
about a hundred times by resonance, whereas the lunar tide, whose period is 
25 minutes longer, will be much less magnified. ‘Taylor’s proof that more 
than one free mode of such oscillation can exist enabled this resonant property 
to be reconciled with the observed speed of propagation of the great waves of air 
pressure sent round the world—more than once—by the volcanic explosion 
of Krakatoa, and with that of the smaller waves set up by the great Siberian 
meteorite. This reconciliation was possible only if the atmosphere had a certain 
type of thermal structure, involving a temperature maximum in the stratosphere, 
and above it an upward decrease to another minimum. 

Above the level of this minimum temperature the air is thought to become 
progressively hotter with increasing height, at least up to the F2 layer, where 
temperatures of the order 1,000° or 1,500°k. may exist. This view is based on 
several independent lines of radio evidence. 

3 1-2 
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Thus the current conception of the temperature distribution is as shown i 
Figure 1, in which the details are to be regarded as diagrammatic, not quantitative 
as regards both height and temperature. A nomenclature appropriate to thi 
distribution is also indicated (modifying slightly some suggestions of my own 
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Figure 1. Atmospheric nomenclature in relation to (a) the temperature distribution, (6) 1 
composition (cf. § 6), and (c) the ionization or the electron density. Pause refers to the j 
of a layer or region, peak to the level of maximum of temperature or electron densit inel 
to a region in which one or other of these quantities is increasing upwards, decline ee regy 
of upward decrease. ‘The nomenclature is provisional, and international a reement 4 
or some modification of it is being discussed. : ae 
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The height distribution of the air density and pressure is, by and larg 
determined by the equation of static equilibrium dp = —pg dh, hich connec 
the decrement of pressure dp with the increment of height dh, and with t 
density p and gravitational acceleration g. As p=knT, where k denot 
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Boltzmann’s constant and m the number-density (or number of molecules per 
unit volume), and p=nm, where m is the mean molecular mass, 

dp mg dh 

| pled vail? 
if H=RT/mg=RT/ Mg, where R denotes the gas-constant 8-3 x 10’, and M the 
mean (chemical) molecular weight. Thus, if H were constant, the variation of 
pressure with height would be given by p/p,=e ””, where py is the value of p 
at zero height; actually H, now called the ‘scale height’, is not constant, because 
T and g vary with height; so also does m in certain ranges of height. The 
variation of g is known: g/g)=a"/(a+h)?, where gy denotes the ground value of 
g, and a the earth’s radius; over the range of height now of interest in atmospheric 
physics, which is at least up to 300km., the variation of g is not negligible; the 
reduction at 300km. is 8-8%. 

To calculate p and p as functions of height it is further necessary to know T 
and m. Paneth has shown that at the ground m is constant all over the earth 
within very narrow limits, although in certain regions notable amounts of helium 
enter the atmosphere from the ground. The mixing action of the winds will 
maintain this uniformity of composition at least up to the tropopause. ‘Through- 
out the troposphere, and indeed throughout most of the stratosphere, the process 
of diffusive separation is very slow. 


wor LHE CHEMICAL COMPOSITION OF THE UPPER ATMOSPHERE 

Not long before the second World War E. Regener and Paneth obtained 
samples of air, by means of balloons, from the layer from 20 to 25km. height, 
in order to find whether there is any appreciable diffusive separation of the 
constituents in that region. Their results in different ways agreed in indicating 
a slight separation; Regener found a small reduction in the proportion of oxygen, ~ 
one of the heavier constituents, and Paneth found a small but larger increase 
in the proportion of helium. 

Since the war intense efforts have been made in the United States to extend 
the range of direct investigation of the upper atmosphere by means of rockets, 
partly German V2’s and partly new American types. On account especially 
of the high speed of the rockets and the low air density, the problems of measure- 
ment and investigation are extremely difficult, but the difficulties are being 
attacked with great energy and skill. Among the chief government agencies 
concerned in this work are the U.S. Naval Research Laboratory, the U.S. Army 
Signal Corps, and the Air Force Geophysical Research Directorates. One of 
the many projects undertaken is the collection of air samples from great heights. 
The University of Michigan, under contract with the U.S. Army Signal Corps, 
has obtained many such samples. Some of them, taken at a height of rather 
more than 60km., and analysed at the University of Durham by Prof. Paneth 
and his colleagues, have shown that the helium content of the air there is not 
appreciably different from that at the ground. ‘This and other evidence obtained 
in the same analyses is in conflict with the pre-war results from the 20-25 km. 
levels, which at the time were accepted as reliable; the discrepancy has not been 
explained, and unless this can be done it seems to call for further investigation of 
balloon samples of air, which should now be obtainable from 40 km. or even 
higher. As regards the high-level conclusions, they have been challenged by © 
some published results for the ratio of NN nitrogen molecules to the main 
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14N14N type, from the same samples; the determinations, made by mass spectro- 
graph, gave a perceptible reduction of the proportion of the rarer isotope, as: 
compared with the ratio at the ground. But further mass spectrograph deter- 
minations from the same samples do not confirm these published results (s 
Hagelbarger et al. 1951). | 
At present, therefore, the balance of evidence indicates that the mean molecular: 
weight M of the air remains constant up to at least 60 or 65km. Work is activel | 
proceeding to obtain samples by rockets from still greater heights. New problems: 
of sampling will arise at the levels, about 100 km. and above, where, there is good 
reason to believe, the oxygen is partly dissociated into atoms. ‘These may oxidize‘ 
the entry tube and containing vessel, and in any case the atomic oxygen cannot bec 
brought down to the ground in atomic form. It may be advisable to place inr 
the entry tube some substance that will absorb the atomic but not the molecular 
oxygen, so that they may be distinguishable and separately measurable when the 
vessel is recovered. Accurate values of the O,/O ratio at different levels in th 
ionosphere would be of great interest to atmospheric and radio physicists. ! 


§6. PHOTOCHEMICAL ACTION 

Besides the opposing influences of mixing and diffusion, respectively tendin 
to maintain and to alter the composition of the air, there is the important influenc 
of the solar radiation absorbed at different heights, which also tends to modi 
the composition. One constituent whose presence is due to such photochemicak 
action is ozone, which cuts off the solar spectrum at the ultra-violet end. Th 
total amount is of the order 3 atmo-millimetres, that is to say, if all the ozone i 
any vertical column of air were collected at the bottom of the column, at norm 
temperature and pressure, its thickness would be about 3mm. ‘This ozone is! 
distributed with a concentration which increases upwards from the ground u 
to about 35km., above which the concentration decreases. Though ozone is 
rather stable gas in the conditions of the free atmosphere, the chemical processes’ 
responsible for its formation and destruction are clearly more rapid than th 
mixing action of turbulence, which continually tries to blur out the gradients of 
concentration. Similarly in the ionosphere the ionized layers and the atomic 
oxygen region are able to maintain themselves against the tendency of both mixing 
and diffusion to modify their distribution. 

I have lately suggested the name turbosphere for the region from the ground 
upwards in which mixing is more powerful than diffusion. The upper boundary 
of the turbosphere, which will of course be only roughly defined, may similarly 
be called the turbopause. Above this level the composition of the air, and 
consequently the value of M, definitely changes; but owing to photochemical 
action (chiefly the dissociation of oxygen) M may change appreciably at levels 
below the turbopause. The level at which M ceases to be approximately constant 
may be called the homopause, the upper boundary of the homosphere, the regior 
in which M can be treated as constant. ‘The overlying region may be called the 
heterosphere. Heterospheric air is materially different from the homospheri¢ 
air we know, and not only as regards small constituents like ozone, which, thougt 
of great interest and importance in many respects, are too rare to modify M 
appreciably. ‘Theory is not yet able to assign levels to the turbopause and homo- 
pause, which must be determined by observation (directly or inferentially); sc 
also must we rely on observation for the value of M for heterospheric air. If the 


ye 
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ratio of nitrogen to total oxygen by mass remains the same in the heterosphere 
as at the ground, complete dissociation of the oxygen would reduce M from 
29-1 to 23-8. 4% 

Thus on account of M and g, the scale height H, which at the ground is about 
8km., may be increased at, say, 300 km. height by about one-third. The change 
of temperature with height, however, has much more influence upon H. 


§7. THE MEASUREMENT OF TEMPERATURE AT HIGH LEVELS 

It is very far from easy to measure the temperature in the upper air, even with 
instruments carried by balloons, on account of the low density of the air (and 
especially by day in the sunshine); but this can now be done, up to heights of 
over 40 km., by the new large neoprene balloons used by the U.S. Army Signal 
Corps. Above this level, measurements have been made by means of rockets, 
whose rapid motion and disturbing influence makes the determination of temper- 
ature from them a most difficult technical problem. Several different ways of 
doing it have been tried by our American colleagues; one of the most promising 
is by the measurement of the conical angle of the shock wave extending from the 
nose of the rocket, in conjunction with the Taylor—Maccoll theory of the wave. 
Another method uses grenades ejected from the rocket at successive heights. 
They explode in the free air, and the time of travel of the sound wave to the 
ground is measured. In this way the mean temperature of the air between the 
levels of successive explosions can be determined. 

These and other methods used are in fairly early stages of development, and 
it is still uncertain whether the differences between their results are due to real 
differences in the air temperature distribution on different occasions, or to 
undiscovered errors affecting one or more of the methods. ‘The results, which will 
not be given here, have been compared with the tentative ‘standard’ atmosphere 
specified a few years ago by the U.S. National Advisory Committee on Aero- 
nautics (N.A.C.A.); the main conclusion is that between about 30 and 70 km. the 
direct observations indicate rather lower temperatures. 

It would be highly advantageous if reliable methods were available to deter- 
mine the air temperature distribution up to great heights by ground-based 
observations, such as radio physicists use to explore the ionosphere. One 
such method, already mentioned, namely the observation of anomalous sound 
propagation from explosions, has been much used in America in recent years, 
also in India and elsewhere. But its height range is limited, very great pre- 
parations are required, and the results are affected by high-level winds and are 
difficult to interpret unless the propagation is known to be very nearly symmetrical 
in all directions. 

Another method of much promise, for places where the night sky is often 
cloudless and the air clear, is the observation of the light scattered from a powerful 
well focused modulated searchlight beam. ‘This method is being actively 
developed at the Geophysical Institute of the University of Alaska; the 
preliminary results give the temperature up to about 60km., though they must 
be critically examined and if possible checked by some independent method 
before final acceptance. 

The temperature of the stratopause, somewhat below 100km., has not yet 
been measured in situ with proven accuracy. Above that level radio measure- 
ments suggest values of the scale height and of the collision frequency which 
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imply that in the upper ionosphere the air is spread out in a manner corresponding 


to temperatures that increase to high values. 


The bulge in the temperature graph for the stratosphere is reasonably _ 


attributed to absorption of solar radiation there by ozone; but the theory of the 
heat balance in the upper atmosphere is still rudimentary. D. R. Bates, in a 


paper recently presented to this Society, has studied the energy balance in the ~ 


ionosphere; he finds it difficult to explain why the air there is as hot as it is generally 
supposed to be. He suggests that there may be an unsuspected supply in energy 
of ionization of some quickly recombining constituent, which does not contribute 


appreciably to the observed ionization (itself still in many respects not understood). | 


§8. WINDS AT HIGH LEVELS 


I have referred to winds in the upper atmosphere, as shown by the deformation | 
of meteor trails. Only rarely are these trails sufficiently extensively and reliably — 


observed to give values of the wind at various points. The most remarkable — 
case known to me is that of 23rd February 1909, for which J. E. Clark collected _ 


and discussed a large number of observations from both sides of the English 
Channel; it lasted for over two hours. The main difficulty in deriving wind 
results from such occurrences is their unexpectedness and the fewness of reliable 
and constant observers. D. R. Bates has suggested that it may be possible by 
means of rockets to lay a trail of sodium vapour in the upper atmosphere, which 
may be luminous and remain sufficiently long for the determination of the wind, 
giving, in fact, the equivalent of a long-enduring meteor trail at a predetermined 
place and time. It is hoped that this experiment will ere long be made in 
America. If successful, it should also throw important light on the photo- 
chemistry of the upper atmosphere. 

However, wind determinations from meteor trails have been made by radar 
methods by Manning and his colleagues at Stanford University, though the 
method gives (by Doppler effect) only the line of sight component of the wind. 
In Australia, in this country and in America it has been found possible to 
determine ionospheric winds by the motion of irregularities in the ionosphere. 

The lunar tidal currents in the various ionospheric layers have also been 
extensively studied, following the pioneer work of Appleton and Weekes. Both 
the vertical tidal motion and the associated changes of electron density have been 
considered, but the data do not lend themselves easily to the determination of the 
horizontal tidal motion. 


§9. ATMOSPHERIC EMISSION SPECTRA 
Recent years have seen great progress in the exploration and interpretation 
of the spectra of the airglow (formerly known as the light of the night or twilight 
sky) and the aurora. Besides the main lines and bands contributed by atomic 
oxygen and molecular nitrogen, there are others that reveal the presence of rare 
but interesting constituents such as atomic sodium, seen in the night airglow 
(and enhanced at twilight, apparently by ultra-violet excitation); its continued 


luminosity throughout the night is still not yet definitely explained. Last year 


saw a very important addition to our knowledge of these rare constituents, when 
Meinel of the Yerkes Observatory showed, with the aid of Herzberg, that the 
infra-red radiation near 10,400. (which is of exceptionally high intensity, and 
was formerly supposed by some writers to be due to atomic nitrogen), is emitted 
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by hydroxyl, OH. An important study of the reactions associated with the 
OH molecules in the upper atmosphere has been made by D. R. Bates and 
M. Nicolet, who have estimated the amount and distribution of this gas and of 
the associated gases H, H;, and H,O at high levels. 
The degree of dissociation of nitrogen in the ionosphere is at present a 
‘controversial question; my own opinion continues to be that there is very little 
atomic nitrogen at least up to and including the F2 peak. This opinion is based 
on the observation of the band spectrum of the nitrogen molecule in sunlit 
aurorae up to much higher levels, 800 or even 1,000 km. 


S10" THE ENTRY OF SOLAR HYDROGEN INTO THE ATMOSPHERE 

The past year has brought another spectral discovery of intense interest, 
also by Meinel, though foreshadowed by Gartlein and Vegard. This is the 
observation of highly Doppler-displaced lines of atomic hydrogen in the auroral 
spectrum, when looking nearly along the auroral rays. ‘The maximum dis- 
placements are to the ultra-violet, and indicate the downward passage of hydrogen 
atoms with speeds exceeding 3,000km/sec. This affords the most direct proof 
yet available for the corpuscular theory of magnetic storms and aurorae, for which 
the previous evidence, though strong, was indirect. ‘There is still much obscurity 
as to auroral theory, but it seems probable that the hydrogen atoms thus observed 
to enter the atmosphere travelled from the sun with a decidedly lower speed, 1,000 
or 1,500 km/sec., and in the form of stripped atoms, namely protons, in company 
with equal numbers of electrons and smaller numbers of atoms and ions of the 
minor constituents of the solar atmosphere (as first proposed in this connection by 
the present Lord Cherwell). Near the earth the protons are in some way 
accelerated and guided along the earth’s lines of force to enter the atmosphere 
in auroral latitudes. In their passage through the air the protons ionize and 
dissociate and excite the air molecules and atoms, and occasionally themselves 
capture an electron, which completes the hydrogen atom, and enables it to make 
its presence and motion known by its radiation and Doppler effect. ‘This 
discovery certainly opens an exciting and very important new chapter in auroral 
research. 

§11. NUCLEAR CHEMISTRY AND THE ATMOSPHERE 

Another new phase of atmospheric physics now being developed is the 
influence of cosmic rays on the atmosphere. ‘The ionization they produce is 
unimportant compared with that due to solar radiation, but the rays cause nuclear 
reactions and changes of isotopic ratio in the atmospheric gases which are of 
great interest. 

A somewhat analogous topic is the addition of gases to the atmosphere from 
the ground, particularly the products of radioactive decay. ‘T'atel has recently 
estimated the age of the atmosphere on the basis of the abundance of the radio- 
active potassium isotope 4°K in the earth, and the amount in the atmosphere of 
the decay product *°A, an isotope of argon. 


§12. IONOSPHERIC VARIATIONS 
Turning again to higher levels, radio exploration has shown that the ionosphere 
has a normal régime of daily varying ionization in layers ionized mainly by solar 
radiation. But occasionally this régime is disturbed, perhaps only for an hour 
or so, as in radio fadeouts, or for a day or a few days, as in the much more intense 
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disturbances known as ionospheric storms. In extreme storms the = ave | 
rapidly rises and seems to disappear into outer space, whereafter ee Re | 
may be re-created at the normal level within an hour or two. itt < simp 
interpretation of the radio records is correct, there may be loss of air to free space 
outside the earth at a rate which, though operating only for short times at long , 
intervals, may be more important than the continual slow leakage of fast-moving 
atoms and molecules usually alone considered in the theory of escape of planetary | 
atmospheres. The latter, however, suffices to explain the absence of a vo 
layer at great heights, and also the loss of most of the helium that has entere | 
the atmosphere from radioactive transformations within the earth’s crust during | 


its long history. 


§13. ELECTRIC CURRENTS: THE AURORAL ELECTROJET 

I do not possess the knowledge required to give even an outline of the course | 

of an ionospheric storm; instead I will devote my remaining minutes to some » 

remarks on the electric currents that flow in the ionosphere as revealed by the » 
geomagnetic variations at the earth’s surface. 


Figure 2. Atmospheric electric current systems responsible for the daily magnetic variations: fer ° 
the equinoxes (a) and for the June solstice (b). In each case the left-hand diagram | 
shows the current system as viewed from the sun, and therefore flowing over the sunlit © 
hemisphere; and the right hand diagram shows the current flow over the night hemisphere. 
The earth is shown with the north pole uppermost. The lines shown are lines of current 
flow, and between adjacent lines the flow is 10,000 amperes. 

There is a regular régime for these currents, as shown by the regular daily 
magnetic variations, in which both the sun and moon play a part. The magnetic : 
data clearly indicate the form and intensity, but not the height, of the current- - 
sheet. Figure 2 shows the system of currents associated with the solar daily - 
magnetic variation, as viewed from the sun or, for the night hemisphere, from . 
an imaginary anti-sun. At times of magnetic and ionospheric storms, however, 
the normal régime is disturbed by the addition of a much more intense system | 
of currents very differently distributed over the earth, as shown in Figure 3, , 
as seen from the sun and from above the north pole. ‘This system differs from | 
that of Figure 2, (1) in being as intense over the night as over the day hemisphere, 
(11) in being most intense in the polar regions instead of over the middle belt of ' 
the earth, and (ili) in being exceptionally concentrated along the line of the auroral | 
zone. By analogy with the jet streams discovered in the upper atmosphere by | 


meteorologists in recent years, I suggest that such a concentration of electric 
current be called an electrojet. 


§14. THE EQUATORIAL ELECTROJET | 

But this auroral electrojet is no longer the only one known. The regular. 
régime on quiet days also includes an electrojet, though it is not shown on 
Figure 2, which is based on an analysis of magnetic data made in 1918. Four 
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years later the Department of Terrestrial Magnetism of the Carnegie Institution 
of Washington founded a new magnetic observatory at Huancayo in Peru, on the 
magnetic equator, which is there 12° south of the geographical equator. It soon 
appeared that at this station the daily variation of horizontal force H was abnormally 
great, with a range of the order 100y instead of 30 or 40y as elsewhere near the 
equator. In 1948 the International Association for Terrestrial Magnetism and 
Electricity appointed a Committee to investigate over what range of latitude 
this abnormality extended. The daily variation of H has now been measured 
at several low latitudes along meridians in Africa, India and America. The most 
extensive set of observations is due to Dr. A. A. Giesecke, Director of the 
Huancayo Observatory. His results show that the abnormality covers a range 


Figure 3. The atmospheric current system that could produce the surface magnetic field during 
amagnetic storm. As in Figure 2, the lines are lines of current flow spaced at 10,000 amperes 
apart. The left-hand diagram shows the system as viewed from the sun; the flow over the 
night hemisphere is similar. The right-hand diagram shows the system as viewed from above 
the north pole. The part of the system consisting of the auroral electrojets and the flow over 
the polar caps within the auroral zones is really situated in the atmosphere, the part between 
the zones may correspond to currents not actually flowing in the atmosphere. 


of only about 10° of latitude. It is not difficult to infer from this that the 
abnormality at Huancayo—paralleled also elsewhere on or near the magnetic 
equator—is due to an equatorial electrojet, which grows in intensity from dawn 
to about 11 a.m. local time, and then decreases again. ‘This is superposed on the 


- normal current flow already illustrated. ‘The electrojet must have a return 


flow, probably both to the north and south, and it will also induce currents in the 
earth. ‘The field observed at the earth’s surface includes the effects of all these 
currents, and is certain to be much more spread out than the electrojet itself, 
which probably has a width of only about 2° or 3°, like the auroral electrojet. 
It must also cause abnormalities in the daily variations of the vertical force and 
declination at stations to the north and south of those that show the abnormality 
in H, and the study of these further irregularities will enable us to determine the 
height and intensity of the electrojet. It seems to me that both Bombay and 
Manila show this abnormality in vertical force. ‘lhe Magnetic Association meets 
again this year, and it is to be expected that plans for further work on this 
interesting phenomenon of the upper atmosphere will be made. 


§15. CONCLUSION 
In conclusion, I hope I have shown by this discursive survey of many 
phenomena, ranging widely both in time and subject, that the field of upper 


atmospheric physics is one that rests on many techniques of observation, touches 


many branches of physics, chemistry and mathematics, and is in the full tide of 


active progress. 
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ABSTRACT. This paper contains a general discussion of electron drift in gases according 
to the method of free paths in which a fallacy contained in some earlier treatments is avoided. 
A general formula is derived for drift velocity in an electric field which is an arbitrary 
function of the time, and from it formulae for drift velocities in the following special cases 
of importance are deduced: drift in a constant electric field, in an alternating sinusoidal 
field, in a constant electric field with an applied magnetic field, in an alternating electric field 
with an applied magnetic field. Formulae for the electrical conductivity of an electron-laden 


gas follow from those for the drift velocities. Applications of the formulae to ionospheric 
and other studies are mentioned. 
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§1. INTRODUCTION 
N this paper a general formula is derived for the drift velocity of the centroid 
] of a group of electrons moving through a gas in a uniform electric field which 
is an arbitrary function of the time. In order to extend the generality and 
usefulness of the treatment it is supposed that a constant and uniform magnetic 
field is also present. 

The current density and conductivity of the gas are then obtained immediately 
from the drift velocity. 

Special cases of the general formula possess a wide range of application, as for 
instance in the theory of the electrical conductivity of metals and semiconductors, 
of the conductivity and permittivity of an electron-laden gas, of the propagation of 
electromagnetic waves in an ionized atmosphere and of the Hall effect. The 
general formula, and the special cases implicit in it, are derived by a correct 
application of the method of free paths which avoids an error which is often made 
in the derivation of formulae for drift velocities by this method. 


§2. NOMENCLATURE AND SUMMARY OF FORMULAE 
It is convenient in this section to define the chief symbols employed throughout 
the paper and to give a summary of some of the formulae to be derived in later 


sections. 
2.1. List of Chief Symbols 


A bar above a symbol or expression denotes its mean value taken with respect 
to the agitational speeds U. W =drift velocity, 1,, i,, i; = coordinate unit vectors, 
U=Uu-=agitation velocity of an electron, /=mean free path of an electron 
(supposed to be independent of U), p=e/m=the specific charge of an electron 
(e is treated in formulae as a positive quantity, that is to say the value 
e= —48x10-E.s.u. would be used in substitution). By=a uniform and 
constant magnetic field, w = —ypB, (y is a constant determined by the system of 
electric units employed. Its value in the chief three systems is given in the table 
in the Appendix; since p is negative, w is positive for electrons), 8&2 = —ypB, (the 
gyro-angular velocity, |Q|=w), p=an angular velocity of an alternating electric _ 
field, v= U/J=collisional frequency at speed U, v= U/l, t=time measured from 
the instant of application of the electric field E(t). 

Other symbols are introduced in the course of the discussion. 


2.2. Formulae for Drift Velocities 


It will be shown later that the velocity of drift of the centroid of a group of a 
large number of electrons in an electric field E(t) =i,X(t) +i, Y(t) +i,Z(¢) which 
as an arbitrary function of the time is expressible in the form 


W= 518,48, +8, +i) ae (1) 


in which the quantities S;,(k=1, 2, 3, or 4) are displacements which are defined in 
the sequel. It proves to be the case that the displacements S, and S, contribute 
only terms that are evanescent with increasing ¢ and that the persisting terms in W 
are derived from S, or S, or both. We are therefore usually chiefly concerned 
with S, and S, and it is convenient to state the formulae for them only. 

Let the electric field be E(t) =i_X(t) + i, Y(t) + 1,Z(¢) and let the applied mag- 
netic field be B, =i,Bo, that is to say its direction is that of the Oz coordinate axis. 
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Let the components of E(¢) normal to By be combined into a complex quantity 
E(t) =X(t) + Y(d), then as will be shown, the components of the displacements 
S, and S, are to be obtained from the following formulae: 


SiatiSy=or | exp(—on){ [| ["" Bxlt—e-H) exp (jut dt | de 
ae (2) 


Sear i, exp (— vn) { I, | ee g—h) ah | as| dy tlre (3) 


t t—n 7 ; cn—g 
Sve+jSoy=p¥ | exp(—) j— “PICs Ex(r-+-g —h) 
/0 0 0 3 0 


x exp (jwh) ah | as ar | ah 4+ koe Geeeeeee (4) 
Sete [ exp (— 17) lilt [ ["Ze+n _g-h) ah | az ar | de 
ee (5) 


In these formulae the symbols p, v, w, t, Ex and Z have the meanings given to 
them in §§2.1 and 2.2, whereas h, g, 7 and 7 are dummy variables. 


It is to be noticed that, according to equations (3) and (5) the components of 


W parallel to By are independent of By since w does not appear in these formulae, 
and that (3) and (5) are special cases of (2) and (4) with w=0. 

It is important to note that in deriving these formulae it is assumed that the 
mean speed of agitation U of the electrons and the distribution of the speeds U 
do not change appreciably during the interval¢. Thus the formulae are applicable 
to the following situations: drift in a field such that the ratio (field strength/gas 
pressure) is small so that the electrons effectively retain their thermal agitation 
energies; drift in a steady field E(t) =constant after a time ¢ such that the mean 
agitational energy }mU? of the electrons is constant; drift under an alternating 
field whose frequency p/27 is sufficiently great in comparison with the time 
constant of decay of the energy 4mU? (when the field is removed) that the 
fluctuation of 4mU? about its mean value during a cycle of the alternating field is 
relatively negligible. 

Because the general formulation is of less immediate practical use than some of 
the special cases that are implicit in it, we proceed to give in summary the formulae 


for the drift velocities W in some of these special cases which have a practical | 


interest. 
2.3. Summary of Particular Instances of Electron Drift 

The following special formulae are derived from equations (1) to (5). 

(a) Drift in a uniform and constant electric field E(t) =constant =E with no applied 
magnetic field Bo. 

If in the general equations we put E(t) =constant, w=0, it will be seen that 

the terms that are not evanescent ast increases give (see the end of sub-section 3.2), 
W= 3Epl(U), aaveealG) 

a well-known formula (Lorentz 1916, Becker 1933, Townsend 1936). 

The formula, more often than not, is given incorrectly in the literature, with 


the factor 2/3 replaced by unity. The factor (U-) is the mean of the reciprocals 
of the speeds U. 
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If the number of electrons per unit volume is , the current density J and 
electrical conductivity are given by J=cE= 2(ne2/m)l(U-)E. 
(6) Drift in a uniform sinusoidal electric field with By =0 (Huxley 1937 b) 

Put w=0, X=E,cospt, Y=E,)sinpt, whence Ey=E,exp(jpt). Then 
W=W,+jW, 


rd v 5 i . f 
= 308 ep wip i rte te * aaah |r 
By aa ' 
Be AH) Nlexp (ine) nen, (7) 


where x=v/p and 


x 1 ’ 
d(x) = Lox E + i. ; (x) = Faye 3 E 1 = SoG (8) 
It follows that W,= : = [A(x) cosptt+p(x)sinpi]. =... (9) 


When v Sf, x >1, h(x) > (1/x) , v=o iS), (x) > 0, and the formula for W,, 
agrees with formula (6) as it should. When v<p, x<1, (x) > 0, (x) > 3/2, 
We (pi aip\ suv pt ame i ihe {ok ss (10) 
a formula which is derived immediately from the equation of motion of a free 
electron and the condition (U) =0. 
(c) Drift in a uniform and constant electric field E with an applied magnetic field By 
(Huxley 1937 a, Townsend 1937, 1947). 
Put B, =1,By, w = —pBy, Ey =X+jY=constant. Then the general formulae 


give Ex — = 
WtjWi=4p— (A) +A wee (11) 


in which y=v/w and ¢(y) and #(y) are the same as the functions that appear in 
equation (7) but with x replaced by y. It is evident from the general equations 
that W, is obtained from equation (6) with || replaced by Z. 

The formulae for W,, W,, and W, are therefore 


= H(p/e)L X86) — VIO) | 


Wy = Helo XM) + VSO) fo waves (12) 
W, = 3pZ(T) | 
(d) Drift in a rotating electric field with an applied magnetic field normal to the 


electric field. 
Let B,=i1,B), E=1,X+i,Y, X=, cospt, Y=£,sinpt, Ey=E, exp(jpt). 
With these postulates, the general formulae give 


1 J(w—p) : 
W,+jW, = 3 pk | 1- Ae \ tas oti. 13 
J = §p v—j(w—p) v—j(w—p) exp (jpt) ( ) 
That is to say 


WatiW = t08s| soe Bena 
Hit - soe | exp(jpt). weve. (14) 


848 De Ge, Phesley | 


It is convenient to introduce the symbols 


i v (o—p)* \ 
Mo) I= eo pP ( * P(o-pP 
a Chet 18 2 ot 
a aA era ies) a Bape 
then W,,+jW, = 2pEy{Fl(w —p), vy) +ifl(@ —p), vi} exp (spt). . +... (16) 
If instead of an electric field Ex = E, exp (jpt) which rotates about B, in a positive | 
sense, we consider the oppositely rotating field Zy=E,exp(—ypt), then it is | 
necessary only to replace p by —p in (14), (15) and (16) in order to find 
(W,,+jW,); that is to say, in this field, | 
W..+jW, = BeBe Fl(o +p), >) +ifl(o+P), 1} exp(— jp). s+... (17) 
(e) Drift in a uniform sinusoidal electric field with a magnetic field By (Huxley 1937b 
(incomplete treatment), 1940). : 
B,=i;By), E=i,X, X=E,cospt, Z=Z,cos(pt+f). The field X can be 
resolved into a pair of oppositely rotating fields as follows: 
E=(E,/2)[exp (spt) + exp (—jpt)], 
whence, from (16) and (17) 
W+jW, = $s Bl Fl(w—p), v) + fl —p), vl} exp (pt) : 
+{Fl(w +p), -}+7fl(w +p), vifexp(—spt)]. ....-. (18) 
In order to obtain compact expressions for W,, W, and W, we further extend the : 
symbolism as follows: 


Slo, p, vy] = F[(w +p), vy] + Fl(w—p), v] | 
s[w, p, vy] =fl(w +p), v] +f[(w —p), v] 


| 
| 
{ 
| 
\4 
| 


D[w, p, vy] = Fl(w+p),y]-Fl(w—p) vy) fo om) 
d[w, p, vy] =f{(w +p), vy] —f[(w —p), v]. J 
It follows from (18) and (7) that 
W,=4pE,[S cos pt+dsin pt] 
W=%4pE|scospi—Dsinpt| (20) 
W, = ¥pZLF(b, ») cos (pt + 8) +f(P, »)sin (pt +f)]. | 


In the more general case in which X=X,cospt, Y=Y,cos(pt+«), 
Z = 7,003 ptt f) and are represented in the notation of complex quantities as 
X=X,exp(jpt), Y= Yoexp(jpt +), Z= Z, exp (jpt + B), then 
W,=3p{[S—jd]X —[s+jD]Y} 
W, =}olls +jDIX+[S—ja]Y} 
W,= 3pZ{F(p, v)—If(b, »)} J 
the real parts of W , W, and W, being the physical components of the drift velocity. 
The matrix of the complex conductivity of the gas is 


Ong Ony 0 
lol] = lo, oH oes (22) 
0 0) o 
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in which Oxy = (ne?/3m)(S —jd)=o,, 
Cyn = — OF gy = (ne?/3m)(s +jD) Cameal ea} 
o,, = (2ne*/3m)[F(p, v) —if(P, v)]- 
‘The current density is te ay 
Oy ee ere Lee ke ee (24) 
J, Z | 


It remains to deduce the general formulae (1) to (5) of sub-section 2.2. 


§3. THE METHOD OF FREE PATHS AND THE THEORY 
OF DRIFT VELOCITY OF ELECTRONS 


3.1. General Assumptions 


In the more rigorous discussions of the motions of free electrons among 
molecules of a gas it is best to suppose that electrons and molecules interact as point 
centres of force. ‘The forces are supposed to be of relatively short range so that 
for the greater portion of its life an electron moves in rectilinear motion, except 
when it approaches a molecule which deflects it through some angle whose magni- 
tude is determined by the precise circumstances of the encounter. Formulae for 
the diffusion coefficient and drift velocity of electrons based on this picture of 
electron motion, although theoretically satisfying, suffer in practice from analytical 
complexity and an element of uncertainty concerning the precise laws of force that 
operate between electrons and molecules of various gases. 

In the method of free paths a simplified situation is envisaged in which the 
electrons move among the molecules in random agitational motion along paths 
which comprise a succession of straight segments (when the magnetic field B, 
is absent) with discontinuous changes of direction between adjacent segments. 
The lengths of the straight segments or free paths are distributed at random 
according to the law that the proportion of a large number of free paths that 
exceeds s in length is exp (—s//), where / is the mean free path. 

The directions of the free paths are supposed to be distributed at random in 
space, that is, the agitational velocities of the electrons are distributed at random 
with respect to direction. ‘The mathematical expression of the assumption of 
chaotic agitational motion is the following: 

Let the velocity of agitation of an electron at any instant be 

CET Ons Ob Eee) 


Zz? 


then U2=U,2=U2=U23 | 
ee EST Wee oe aoc 5 (25) 
l= eae U0), 


where the average is taken over any large number of free paths described at speed U. 

Thus electrons and molecules are regarded as particles with dimensions 
negligible in comparison with the mean free path. ‘They collide as rigid bodies 
bounded by fixed geometrical surfaces. 

In order to preserve the chaotic character of the agitational motion it is not 
necessary to postulate any particular shape for the molecules, ‘They may be- 
regarded as spheres, regular or irregular polyhedra or as bounded by smooth 
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surfaces of arbitrary form. Provided the axes of inertia of the irregular bodies are « 
directed at random the chaotic agitational motion of the electrons will be preserved. 
It is obvious that this indifference of the method of free paths to the form of the : 
molecules renders it impossible for the method to provide information about the | 
mean energy lost by an electron in collisions with molecules. The method however - 
when used to interpret the results of experiments on the diffusion and drift of | 
electrons permits this quantity to be assessed. 

If the mean cross section presented by a molecule to an electron is A and ifn 
is the number of molecules per unit volume, then the mean free path is /=1/nA. | 
It is often convenient to consider the molecules to be spheres of radius a=(A/7z)"?.__ 

According to the method of free paths the mechanism of electron drift through 
a gas under an applied electric field E is the following: Since the field E exerts a | 
force Ee on each electron the free paths are, in general, no longer straight because — 
each electron receives a small displacement in the course of a free path. The 
vector sum of such displacements for all the electrons of a group taken over many _ 
free paths is not zero, but is in general a function of the time. The centroid of the 
group is therefore displaced in time, and the time rate of change of the vector 
joining centroid and some fixed point is the required drift velocity. Thus the 
procedure for finding the drift velocity W requires an accurate vector summation of 
the elementary vector displacements caused by the field E in the course of the free 
paths. 


3.2. Derivation of Formulae for Drift Velocities—General Method 


It is convenient first to suppose the magnetic field B, to be absent. In the 
absence of an applied electric or magnetic field the path of an electron between 
successive collisions with molecules is a straight line traversed at constant speed U. 
Ina time of flight » the distance travelled iss =U. When the electron moves in an 
electric field E its path in general is curved and ina time of free flight 7, not termin- 
ated in a collision, the distance traversed is different from s=7U. Thus, in the 


f 


Figure, OP =s represents the vector displacement s =U when E the applied field — 
is zero. When E is not zero the path is curved and in time y the electron reaches | 
the point Q. Thus the additional displacement effected by the field E in time n | 
is As,. It is easy to obtain a general expression for As, when the electric field is _ 
an arbitrary function E(f) of the time. 

The vector equation of motion of the electron is d2r/dt2 =pE(t), where r is the | 
vector displacement from O and p=e/m. | 
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Double integration of the equation of motion gives 


OO=nn=nU+p le | i E(h) ah A ge eat (26) 


in which 7 is the time at which the electron passes through the origin at velocity U 
in the direction OP. Since OP=s=yU it follows that 


PQ=As,=p | ‘ [ | E(h) ah | de. So,” NRE (27) 


It is again stressed that As, is the displacement effected by E(t) in a specified time 
of flight. ‘This, however, is not the only type of elementary displacement effected 
by E(t) that we are required to consider. 

It is important to note that the effect of the field E is to curve the free paths but 
not to alter the mean value of the free paths described at speed U nor the distri- 
bution of their lengths about the mean which still conforms to the law that the 
proportion of a large group which exceeds a length s is exp(—s/l). 

Thus, we are also concerned with the problem of finding the displacement As, 
due to E(t) when the length of the free path is previously specified to be s. 

Thus, whereas in the calculation of As, the interval 7 is the same for the dis- 
torted and undistorted paths, in the calculation of As, the lengths of the distorted 
and undistorted free paths are the same, but the time of flight is different for each. 

We shall suppose in all cases that the displacements As, and As, are small 
compared with the undistorted displacements s=OP. 

Reference to the Figure shows that the displacement As, is PR which for small 
displacements is effectively normal to OP. ‘The distorted free path is OR. 

There is an evident geometrical relationship between As,, As, and U=Uu 
‘(u is unit vector parallel to U): 

ASI=O8;—AWaAsS, Monae 4) 7S ae (28) 
Thus, in the example shown in the Figure if the time of flight along the free path 
OP =s when E=0 is n, then that along the curved path OR of length s when E40 
is (7— Ay), where An = RQ/U=PS/U=(u. As,)/U. 

Thus the effect of E is to alter the speed of the particle in the direction of U 
without adding to its displacement in this direction since that is already prescribed 
to be s to the first order of small quantities. ‘The only displacement produced by 
E(t) is As, as shown. Failure to recognize the distinction between As, and As, 
has led to erroneous formulae in many previous treatments. 

Of interest for what follows, is the mean value of As, taken over all possible 
directions of the velocity with free paths of prescribed and fixed lengths. Let the 
vectors As,, As,, U and u be expressed in terms of their components and the unit 
vectors i,, i, and i, with respect to a cartesian system of coordinates, as follows: 

As, =1,Ax,+i,Ay,+i,Az,, As, =1,Ax,+1i,Ay, +1,A20, 
U=i,U,+1,U,+1,U,, uni? +i +i, Zt, 
also, E=i,X+i,Y+i,Z. In terms of components, equation (28) becomes 
As, = 1,Ax, + i,Ay, +i,Az, 
=As,—(U,Ax,+ U,Ay,+ UAz,)(i;U,.+1,U,+1,U,)/U? 


: Dee U2 UUs Oy Ve eset UU, | 
=i,] Ax (1- =) a So | +g E& (1 - ) my | 


: DA OOLFOU, 
+1, Az, 1- U2 <5, . 


Bes 


| 
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In averaging As, over all possible values of U,, Uy and U, with fixed U we use. 
equations (25) and obtain | 
(As,) = $[iAx, +igAy, +igAzy]=$As, = «++. (29) 


nyhencey using (27),  \(As,)= ep | bh | | E(h) ah | 3° oe (30) 


where 7 =s/U. 

We proceed to derive a general expression for the drift velocity Wunder afield E | 
in the absence of amagnetic field. Let the speeds U of the electrons be distributed 
according to some law that the number of electrons in a large group NV whose speeds ; 
are not less than U but do not exceed (U+dU)at any instant is Ny dU= Nf(U) aU. , 
The mean collisional frequency of an electron is 


eames 


y= Uilwherel = | moyau and? We i Uf(U) aU. 
0 


It is convenient to neglect the dependence of mean free path upon the speed U! 
and to suppose that /is the same as/. We write, therefore, 


fe | +AU) au = | »f(U)dU, where »=Ujl. 
0 J0 
The number of collisions made by the N electrons per second is 
Ni= { yNf(U) dU = { N,vaU. 
0 0 


Thus the group of N electrons may be considered to comprise sub-groups + 
N,dU=Nf(U) dU, in which the electrons move at a fixed speed U and with af 
collisional frequency v= U/l. : 

The number of collisions within an elementary interval of time dr is ¥ 


Ni dr =dr ie NyvdU, consequently the number of free paths originating in this § 


interval in which the initial speed lies between U and (U+dU) is drvNy dU. . 
Suppose the field E to be applied at ttme7=0. We proceed to find the sum of the : 
eeberasts due to E along all free paths that exist during an extended interval | 
Sir Se 
Of the vN,dUdr free paths that begin within the elementary interval dr at ¢ 
time 7, some are not terminated within the chosen interval, that is they continue for : 
time 7 such that r+7 >t. The remainder terminate in collisions which occur att 
or before time ¢, that is to say r+ <t. | 
First consider the former type of free path, for which r+ >t. The number 4 

of such free paths is yNydUdrexp[—v(t—7)]. The mean displacement pro- }) 
duced by the applied field E in these unterminated paths at time tf is given by | 
equation (28), that is, 


As, =(As,) =p [ EE E(h) ah | ue Theis (31) 


The sum of these displacements for all possible initial intervals dr that fall] 
within the large interval 0 <r <t is 


S,(Ny dU) =(Ny aU) | As, exp[—v(t—7)] dr 


=(NydU)pv [ exp —v(t—r)] : f | He E(h) ah | as Orin Gian } . (32) 
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where S, is the coefficient of (Ny dU). 

We require also the sum of the displacements along free paths of the latter 
type, that is along free paths originating in an interval dr at time 7 but terminating 
before time + =#, that is to say the duration 7 of the free path is such that 7 +7 <t. 

According to equations (29) and (30) the mean displacement along such a 


ath i Fie 
al (Ass) = ge | | “E(i) dh | Con; re 22 (33) 


The number of free paths starting within dr at time 7 and persisting for times 
which exceed » but do not exceed 7+ dn is Ny dUv* exp (— vy) dy dr. 

Since terminated free paths of duration 7 are distributed at random over the 
interval 0 <7 <t—y, it follows that the sum of the displacements due to E in 


such free paths is ae 
Ny dUv exp (—vn) dy | (As,) dr, 
0 


where (As,) is defined in equation (33). 
It remains to sum over the range 0 <7 <¢ of permissible durations 7 that 
fall within the interval 0<7<t. The sum of these displacements is 


S,(N,dU)=(Ny dU) 3p" i ; exp(—™) he { | ye ( \ E(h) ah) as ar | ie 
eg (34) 


The equation serves to define S,. 

It remains to consider displacements along free paths that start before time 
+=( but persist into or beyond the interval 0 <7 <i. 

The number of free paths that begin within an interval dr’ at time 7’ earlier 
than time s=0 but are not terminated before time t=0 is 


NydU | exp(—r1') dr’ =Ny aU. 
YO 


The number of free paths being traversed at time 7 =0 that continue beyond the 
end of the interval 7 =tis N, dU exp (-— vt), and the sum of the displacements due 
to E in these paths at time ¢ is 
N, dU (As,)exp(—vt)=S,(NydU), ~ oe aee. (35) 

where (As,) is as defined in equation (31) but with 7=0. 

The number of the N,,dU paths whose durations are not less than y but do 
not exceed 7+ dy, when 7 is such that 0 <7 <t, is Ny dUexp(— vm) dn, and the 
sum of the displacements along all such paths terminated before time 7 =¢ is 


i S 
Ny dUv | (As,) exp(—vy)dn=S,(NygdU), se (36) 
0 
where (As,) is defined in equation (33). 
It follows that the sum of the displacements along free paths of all types 
traversed at speed U up to time ¢ is 


(N, dU)S=(N,aU)[S,+8,+8,+8,],  ...... (37) 
where S,, S,, S, and S, are defined in equations (32), (34), (35) and (36). 


The sum of all the displacements due to E for all free paths and the whole 


range of speeds Uis_ .,, ie cs . 
| SN, dU=N i Si ViGSiey 7) ee (38) 
0 0 


. 
a 
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This equation serves to define the mean displacements S of an electron of the | 
group. If at time =0 the centroid of the group of N electrons is at the origin — 
of coordinates, then at time ¢ the centroid is at the end of the vector 


R=s- ee (39) 

and the drift velocity W at time ¢ is ) 
_ dRy _ dS i 2 a ae es | 

ae [S, +S, +8, + ee Se (40) | 

‘The terms d(S;)/dt and d(S,)/dt are evanescent and become unimportant as t | 


becomes large. 

The first two terms in the square brackets in equation (40) not only yield — 
evanescent contributions to W but also give the steady state value of W, which | 
is what is most commonly required. Since the formulae for the displacements | 
S,, (k=1, 2, 3, 4) are distributed within the discussion it is useful to display them | 


aS a group: 


S.=pr | expl—o(e—7){ | [ [BC a as dr, 

Sees i, ee) Rae (f E(h) dh) as ar | pe 
S,=pexp(— vt) fe Be B(i) dh | de, 

S,= gor | Sena { I. l i E(h) an) | az 8 


In practice the formulae are given a more convenient form by making suitable | 
changes in the dummy variables so that every lower limit of integration is zero. 
They are thus equivalent to 


S,-py| exp(—rn){ || |” "Ee-g—h) dh | ds | dn | : 
: 


S,= or" | exp(— mm) Bees ( |" B-+7-g—Wah) as dr | ae 


S=pexp(= 9) i Bh B(h) dh | de, 
S.= for] exp(—m){ || [EG an | dg} dn, 


Formulae (41) are the same as formulae (2), (3), (4) and (5) when w=0. 
To illustrate their use, put E(¢)=constant in equations (41). Then 


P= lv @ fexp(—vt) 
S,=0E | 5-3 at: 
2 ps = he cn 
S,= $B | 1 a " exp ( ay _ raat exp ( a) 


v dv v 


E?? 
S,= a exp (— v2), 


s _ pvE d® exp (—vt)—1 
a 9 dv v E 


/ 


| 
: 
: 
: 
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The terms independent of ¢ vanish in equation (40) and those containing the 
factor exp (— v#) still retain it after differentiation, and are evanescent as t increases. 
The only term in W which is not evanescent derives from S,. Thus when 
tis large, W—> §pE(1/v) = $pEU~?) as stated in equation (6). 
It remains to derive the general formulae for drift in the presence of a 
magnetic field Bo. 


3.3. Electron Drift in the Presence of a Steady Magnetic Field 


It is important to obtain formulae for the drift velocity W when electrons 
move through a gas in the presence of a constant and uniform magnetic field By. 

It is necessary first to derive expressions for the elementary displacements 
As, (in a specified time 7) and As, (in a specified length of path s) as defined in 
§3.2. ; 

The equation of motion of an electron in free flight at velocity V is 
mdV /dt = e[E++(V XB,)] which may be written 


adit (Vi X82) = pee 8 (42) 


in which p=e/m, |Q|)=w=—~ypB, (e positive); that is to say, w is positive for 
electrons. The constant y is defined in the Appendix. 
Choose the direction of By to be that of the axis Oz; let 


E(¢) =i, X(t) +i, Y(t) +i,2(2). 


Equation (42) resolves into 


STaOye=pA,|  V-ar=pY, S=pZ. “3 ...0. (43) 
If y=(x+ jy) then it follows from equations (43) that 
vejuve=piX+jY|=plyy S=pZ.< ....2: (44) 
Equations (44) are equivalent to 
djdt(x exp(—jwt)]=plexp(—jwt)]Exy, %z=pZ. ...... (45) 


At time t=7 we suppose that y=0, z=0, x(7)=)U,+/)U,=Uyexp (ia), 
and z=,U,. Integrate each side of equation (45) between the limits 7 and g 
to find i 
ig) =Ve+iVy=exp (jug) |" [exp (—jeoh) 2h) dh | 


+ Uyexpj{w(g—7) +}, 
V,=8=p| Z(h)dh+0U,. | 


It is convenient to change the limits of integration and to write equations (46) 
in the equivalent form 


ig)=Ve+iV,=p | exp (iol Ex(g—Mah | 
+Uyexpj{fo(g—t)+a}, © = fF — sevaee (47) 
(g)=V.=p | Z(g—h)dh+U, | 
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‘Thus the orbital velocities at time g are ; 


V,.=Uycos[w(g—7)+a]+p [1x@-#) cos wh— Y(g—h)sinwh]dh=U,+v,, 
0 
G=Tt 

V,=Uysin[o(g—7)+a]+p 'p [X(g—h) sin wh + Y(g—h) coswh]dh=U,+2,, 


ryt 
V.=0U.+p | Z(g—h) dh=U,+0, 


where v=i,v,,+i,v, +i,v, is the additional velocity due to the applied field E. 
A further integration of equations (45) gives the following expressions for the 
displacements of the electron from time 7 to (r +7). These are 


ttn rg-7 U. hs iy E 
xrt+n)=e | rf exp (joh)Bx(g—h)dh | dg-+ “PCCD Nesp), 


=p ie | Ss (jh) Ex(7 +9 —8—h) ah | dg+ Un exp (ja)Lexp Gon) AG | 
0 Jo ep 


(r+n)=p [| [” "Zer-+n-g—A)dh | dink. . (49) 


Since all values of the angle « are equally probable and the mean value U, is zero, 
it follows that the mean displacement along a free path in a specified interval of © 
time from 7 to (7 +7) is As, =i,As,,+1,As,,+i,As,,, where 


- = ¥ 
Ase | [ {X(7 +n —g—h) coswh— V(r+9—g—2)sin oh) dh | ae, | 


on f- pn—9 

Asse | [’ [X(r+9-g—M)sinwh-+ ¥(r+q—g—h) cos wh} dh | dg, (50) 
a n-g 

ds.=p || | Ur+n-g—h) dh | dg 


rn 1-9 
whence (As,,+jAs,,)=p | ‘ | i Ex(7 +n —g—h) exp (jwh) ah | dg, .s.\oa 


where Ey=(X+jY). In order to derive formulae for the displacements { 
As, =1,As,,, +i,As,,+i,As,, which occur during a free path of specified lengths — 
it is first necessary to obtain an expression for the alteration in the time of flight © 
along the path s caused by the field E. The tangential speed ds/dt of an electron _ 
along a free path is given by 


ds\2 
(5) =Vi+V 2+ V2=(U,+0,)2+(U,+0,)2+(U,+0,)2 = U2+2U.v+e, 


in which the velocities U and v are defined in equations (48). 1 
It is assumed that the field components .X, Y and Zare not of sufficient strength 
to make the velocity v comparable at any point of a free path with the initial — 
agitational velocity U, consequently the term v2 in equation (52) is discarded as 
a small quantity of the second order of magnitude. With this restriction it 


follows from (52) that dsjdt=U+(U.v)/U.. nee (53) |i 


The length of free path traversed from the instant of collision at time 7 up to time | 
(7+7) is ) 


Conductivity of an Electron-Laden Gas 857 


When the field E is zero the time required to describe a helical free path of length s 
at constant speed U is 7»=s/U, but when E is not zero the path of length s is 
described in a different time, 7’=7—Ay. An approximate expression for Ay 
follows immediately from equation (54), which may be written 


? s pe r(z+8/0) Ty 


— a B=9-An, ie. An = | art Ren (55) 


Thus, whereas the displacement As,, whose components are given in equations (50) 
and (51), is that which is effected in time 7 =s/U, the displacement over the path 
length s is effected in time 7—Ay. ‘Thus As, =As,— Ay d(As,)/dy — Ayn U(t +7). 
The term Ay d(As,)/dy depends upon squares and products of the field 
components X, Y and Z and is therefore of the second order of small quantities. 
It follows that at time g=7+7 


T+ 
As, =As,—AnU(r +) =As,— ae 


Uivider pss (56) 


The mean displacement (As,) effected along free paths s for all initial directions 
of the velocity U is 


es Caine) 
Rens olay U.vde |. ae (57) 


The velocities U and v are defined in equations (48). Thus, at time g 

U,= Uy cos[a(g—7) +a] =)U, cos w(g—7) — )U,sin(g—7), 

U, = Uy sin [w(g—7) +«] = )U, sin a(g—T) + )U, cos w(g—7), 

if = Fue 
It follows that the second term in (57) is the sum of terms containing the products 
2 Uy 0U x 0Uz, oUy 0 U, as coefficients which vanish in the mean. There are 
also coefficients (JU,/U)?, (9U,/U)?, (9U,/U)* each of which has, according to 
equation (11), the same average value of 1/3. 

When these facts are noted it will be found that equations (57) and (47) and 

(50) lead to the following expressions for the components of the vector displace- 


ment As, : 
=, rn rn—g 
As._=p| | { (1-528) | [X(t + —g—h) cos wh 
0 0 
j r= 9 
= ¥(r-+9—g—h) sin oh} di} dg — 4 I [X(r + —g—h) sinwh 
0 
+¥(+9-g—h) cosh] dil dy |, 
As», =p [PE Px +928) cosh — Vor+y—g—h) sine dh 
0 0 


S 
we ¥(r-+n-g—h) cos wh] di as | , 


alle (1 na co528) J" X(r+1-g-h) sin oh 
0 


— n n-g 
As..—$p | oe Ur+n-g-N)dh | dg. 
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The vectors S,, S,, S, and S, are derived from As, and As, by the same procedure | 
as that employed in §3.2. We find, corresponding to equations (41), that the | 
components of S, and S, are given by equations (2) to (5) in §2.2. The formulae — 
for the displacements S; and S, are not stated because they contribute evanescent } 


terms only to W. 
Although only elementary analysis is required, the derivation of the formulae 


for the particular instances from the general formulae, as discussed in §2.3, is | 

somewhat tedious. 
. §4. APPLICATIONS OF FORMULAE 

4.1. Magnetic Deflection of an Electron Stream | 

In formulae (12) put Y=Z=0; then W, = (2p/3w)X4(y); Wy = (2p/3w)Xp(y); 

thus W,/W,=(y)/¢(y). When y? =(v/w? >1, o(y)>(1/y), Hy) > 31/9") and > 

W,,|W,,—> $o(1/v?)/(1/v) = — 2(yBo/X)W(U~)/(U)?, where W is given by (6). 

A method for testing these formulae has been described (Huxley and Zaazou 1949). 

4.2. Hall Coefficient 


Let acurrent flow in a rectangular lamina lying in the xOy plane and whose edges _ 
are parallel to the axes Ox and Oy. Suppose the current to flow parallel toOxand | 
let By be directed along Oz. Since no steady current can flow in the direction | 
+ Oy, it follows that W, in equations (12) is zero. Thus Y= —X¢(y)/¢(y) and E 
W ,,=(2p/3w)X4(y)[1 + (4(y)/¢(y))?]. The current density is J=neW and the 
Hall coefficient TT 

u 3wh(y) 


FT epi) + $0) 1 
When, as in a metal, y=v/w >1, this formula shows that 
R->0-75[(U-2)/(U-)]yB,|/ne. 


When the speeds U are distributed according to Maxwell’s law then 
R—-1-175yBo/ne, a standard formula. 


4.3. The Mean Power communicated by a High-Frequency Electric Field to an 
Electron in a Steady State of Motion in a Gas 
The mean power given to an electron is w=}e Re[XW,* + YW,*+ZW,*], 
where asterisks denote the conjugate of the complex quantity. It follows from 
equations (21) that (Huxley 1949) 
w=(e?/6m)[(Xo2+ Yo)S+2X,YoDsina+2F(p, v)Z2]. ...... (59) 
For instance, if the electric field rotates with constant amplitude about the Oz axis 
(parallel to By) then X»=Yy)=Ey; Z,=0; sina=+1. Thus from (59), with 
sina= +1, ok 2 
Word = ee [S+ D]= 2 moo E(w +p), v], 
and with sin« = —1, 
oa OE 
extraord 3m 
If Y=0 and 8 =0 (plane polarized field) w = (e?/3m)[4SX,2 + F(p, v)Z,2]. When, | 
in addition, By=0 (w=0), w=(e?/3m)F(p, v)[Xo2+Z,2]. These results are of | 
interest for wave propagation in the ionosphere. 


[S—D]= ; “BF —p), v]. 
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4.4. The Refraction and Absorption of Electromagnetic Waves in an Ionized Medium 


In order to illustrate how the theoretical formulae for the conductivity of an 
electron gas in a high-frequency field finds an application in the study of the pro- 
pagation of electromagnetic waves it suffices to consider the simple problem in 

which the applied magnetic field By is absent. (A more general: discussion has 
been given by Huxley (1940).) 

According to equations (7) the conductivity with no constant magnetic field is 

a complex quantity of the form 


CESS ee ee eo (60) 
where, according to equations (8) and (9), 
2Ne@l-—— ON Cw eam 
SSi35 na SO i Ermey Heat Snagit (61) 


in which WN is the number of electrons per unit volume. Since we require 
Maxwell’s equations of the electromagnetic field they are here stated in general 
terms valid for any particular system of self-consistent units. In this notation the 
equations read 


culE=—yB,  curlH=y(D+J)=(D+cE), | 
div D=p, div B=0, j 

math B=p,K,H; D=<,K.E. 
- The values of the constants y, €) and jx) are determined by the particular system 


of units adopted, as is explainedin the Appendix. K,,is the magnetic permeability 
and K, the dielectric constant. In this notation the Poynting vector is 


NOSE I Ser Re amen eee Ci © (63) 
In (62) put K,,=1; p=0; J=cE, and suppose E and H to be the fields of a plane 
sinusoidal wave train. It follows that 


curl H =(jypeg) (1 +2] E=ye ( — z) E =ye K.E, ales ois (64) 


where Koes lajolpes)s ae lw ttare we (65) 


is an equivalent dielectric constant. 
The phase velocity in a vacuum is c = 1/[y(€o449)”] whereas that in the medium 
is v=1/[y(eou,K.’)"2]. The refractive index of the medium is therefore 


nan—jcp=(K. YY P=C—jolegp)¥?. inane (66) 
Since (9)! =(e,K,)"£, it follows that 
Pieenttey ICRP), § 9) | -vmenes (67) 


The absorption coefficient is x, consequently if E = Ey exp (—«s) exp [7p(t —ns/c)] 
then H = Ey €o/to)™ Bs +c payee exp (—x«s) exp|[jp(t —ns/c) —ja], where 
tana=(cx/np). Put o=(a—yjb) in (66) to obtain 


—(cx/p)? =1—b/eop, 2nce/p=Afegp. = saceee (68) 


When «A/2z is small compared with unity, that is to say the wave is not strongly 
absorbed, it follows from (68) that n2=1-—5/e 9p. Furthermore, if the wave is 
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travelling in a region in which the refractive index m is not greatly different from | 


unity, that is the region does not reflect the wave, then «=a/(2e9c). It follows 
from the definitions of a and b given in equations (61) that in these circumstances 


n® =1—(2Ne%/3m)(1/p2)b(a), «= (Ne*/3megc)(1/p)4(#) 


when B=0/p<1; yx) 3/2; p(x) > 2% =(2/p)v =2U/p. 
The refractive index is then given by the standard formula 
n?=1—Ne?/(egmp?), but «2(2vNe"/3eqmcp). ...... (69) 


The usual formula for « is, however, « = Ne?v/(2emcp*), that is to say it makes k too | 
small by a factor of 3/4. Formulae (69) apply to the F region of the ionosphere — 


when p/27 ~3 x 108 c/s. or greater. 


Next suppose x=1/p=1. (x) is a maximum at x =0-78, its value also being _ 
0-78. Since ¢(1)=0-75 and ¢(0-5) =0-72, the maximum is flat and 4(x) does not — 
change markedly in the range 0-5 <x <1. In this range « is virtually independent _ 


of x and is approximately represented by 
k=0°25.Ne?/(egmp). «9 ea 2 | eee (70) 


The usual formula adopted for « is x =(Ne?/2egmcp)[x/(1+x?)]. The function 
x/(1 +x?) attains its maximum value of 0-5 at x =1, consequently, according to this 


formula, near x =1 « is approximately represented by x =0-25.Ne?/(<gmp), whichis _ 


the same as (70). It may be concluded that the usual formulae for « are not 
seriously in error when x = v/p ~ 1, as for instance in the E region of the ionosphere. 


4.5. The Complex Dielectric Constant of an Ionized Gas 
Equation (65) is the formula for the complex dielectric constant of an electron- 


laden gas. ‘The formula is equivalent to 
K, =1 + (b/peo) —j(4/p<o) 


_— ——. 


=1+ (2Ne?/3m)(1/p*€o)b(x) —j[(2Ne?/3m)(1/preo)}b(x), .--- (71) 
in which x=v/p=U/(Ip) and ¢ and % are defined in equation (8). 


§5. CONCLUSION 
The theoretical investigations presented in this paper, of the drift velocities of 
electrons in various circumstances, reveal an underlying unity in the subject which 
was hitherto lacking. What were previously a collection of apparently indepen- 
dent formulae are now seen to be particular examples of one comprehensive general 
formula covering the subject of drift velocities of electrons in direct or sinusoidal 
electric fields accompanied or unaccompanied by a permanent magnetic field. 


APPENDIX 


The constants €9, 9 and y enter electrical theory as follows: The two coulomb | 


laws (ey and ji») are 
P= q199/(47€)K 7°), F=mym,/ (4p) Kp?), 


in which K, is the dielectric constant and K,, the magnetic permeability. It is 
the custom to write « =«9K,, =) K 
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The constant y enters through law of the equivalent shell: y=yut=ypoK yi, 
where x is the strength of the uniform shell equivalent to the current loop 7. 

Asingle relationship exists between the constants in any absolute self-consistent 
system of units: y*e9u9 =1/c?, where c is the velocity of light ina vacuum. Thus, 
if so preferred, y may be eliminated in favour of c. 

The values of €9, #9 and y in the three principal systems of units are the 
following: 


System €o Ho y  c(approx.) yo Yeo 

Gaussian 1/47 1/42 4mfe 3 x 101° An/|c? 1/c 

E.M.U. 1/42rc? 1/47 4a 3x10 4ar 1 

Rationalized M.K.s. 10°/47rc? Ax IOs? ~ 1 SeLO® dor 104" © 4a 10" 
=1/(367 x 10°) 


It should be noted that induction B in general or in M.K.S. units is not the same 
physical quantity as induction B’ (say) in gaussian or electromagnetic units. 
The former is B=»H=,,K,,H; the latter is B’=K,,H, thus B=p,B’. Since 
#49 in general possesses physical dimensions, B and B’ are not the same physical 
quantity. Similarly D=«E=«e,K,E=«,D’. The gyro-angular velocity Q =i,w, 
employed in the representation of drift in the presence of an applied magnetic 
field By, is 22= —ypBp, i.e. w= —ypBy. The general formula for the Lorentz ° 
force is F=yelV x B)]. 
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ABSTRACT. The refractive indices of air and its principal constituents have been ik 


measured at a frequency of 24,000 Mc/s. with a precision comparable with that obtained in 
the optical range. 
cavity resonator first when it is filled with the gas and then when it is evacuated. The source — 
of oscillations used is a Pound-stabilized velocity-modulated oscillator, and its frequency is _ 
measured by reference to a high harmonic of a quartz standard, with a precision of 1 part in | 
108, The cavity is provided with a tuning plunger, which changes the resonant frequency — 
through a range of about 10 Mc/s., and is calibrated to an accuracy of within 1 ke/s. | 

Most of the measurements were made by the frequency-change method and the plunger — 


was used in a narrow region only for a precise setting to resonance. It is possible bya | 


larger movement of the plunger to compensate for the whole frequency change due to the 
removal of the gas and thus to work at a fixed frequency. 

The following results were obtained for (z—1)10® where n is the refractive index at 
0°c., 760 mm.Hg: dry CO,-free air, 288:15+0-1; nitrogen, 294:1+0-1; oxygen, 
266:440:2; argon, 277:8+0-2; carbon dioxide, 494+1; and the value for water vapour 
at 20°c., 10 mm. Hg pressure was 60°7+ 0-1. 

The dielectric constants can be calculated from the relationship pe=n?, where p is the 
magnetic permeability and « the dielectric constant, the values of (u—1)10® being taken as 
0-4 for air, 1-9 for oxygen and zero for the other gases. 

Accurate formulae are given for obtaining the refractive index of moist air at different 
atmospheric conditions, and are reduced to the following simple formula which is applicable 
for normal atmospheric conditions : 


103-49 pen Wigs eo 
T Pi “Se Pe T a Ps, 


(14, »—1)10°= 


where Pj, pe, p3 are the partial pressures of dry air, carbon dioxide and water vapour, t is the 
temperature in degrees c., and T=273+t is the absolute temperature. A value of 
1-839 + 0-002 x 10-18 E.s.u. was derived for the dipole moment of the water vapour molecule. 


§1. INTRODUCTION 
HERE are several different branches of work in the radio-frequency field 
| which require a more accurate value of the refractive index of air than 
is so far available. 

Precision radar, for example, has reached a stage of development at which 
the accuracy is limited by uncertainties in the value of the correction which 
must be applied to allow for the effect of the atmosphere. An accurate value 
is also required for a microwave interferometer which is being developed at the 
National Physical Laboratory for the measurement of length. The present 


operating frequency of the instrument is 24,000 Mc/s. and the measurements | 


described here have therefore been made at this frequency, although the method | 
could be used at any frequency in the microwave region. | 
No previous determinations appear to have been made at 24,000 Me/s.,. | 


although there are a number at lower microwave frequencies and at low radio | 


: ; me il 
frequencies. As shown in the discussion later, none of these has the precision | 


required for the present application and, moreover, they diverge by considerably | I 


The method is based on the measurement of the resonant frequency ofa | 
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more than the limits given by the authors. The values of refractive index obtained 
for dry air by extrapolation from the dispersion formulae established by the 
most reliable optical measurements also differ appreciably, and it is not certain 
that the effect of the resonant absorption by oxygen in the millimetre band can be 
altogether ignored, and therefore that extrapolation is completely valid. 

The most convenient method at microwave frequencies is the measurement 
of the change in the resonant frequency of a cavity when filled with the gas under 
test and then evacuated. Then for gases which do not absorb microwaves 


heavily (k<n) n—1=(f,—f,)lfe=(ue)2—-1 meen a 


where 7 is the refractive index of the gas, » the magnetic permeability, « the 
dielectric constant, k the absorption coefficient, f, the resonant frequency of the 
cavity when evacuated and f, the resonant frequency of the cavity when filled 
with the gas. 

The original aim was to measure m—1 with an accuracy of at least 1 part in 
1000 and since n—1 itself is only 3 x 10-4 the precision of the resonant setting 
must be about 1 part in 10’, although the absolute value of resonant frequency 
is not so important. One method giving the necessary sensitivity has been 
used by Crain (1948) and others. T'wo cavities, one containing the gas, control 
the frequencies of two oscillators in the manner described by Pound (1947), 
and the heterodyne beat frequency between them is measured. The cavity 
containing the gas is then evacuated and the change in the beat frequency 
determined. This technique though beautifully simple is open to several 
objections. In the first place, any change in the frequency of the comparison 
oscillator during the measurements causes a direct error in f, —f,, and to overcome 
this uncertainty in the work described here the comparison oscillator was 
dispensed with and the frequency measurements made absolutely in terms of a 
quartz standard. A second and probably more serious source of error arises 
from the fact that a Pound-stabilized oscillator is not controlled at exactly the 
resonant frequency of the cavity and the point of control varies with the setting 
of other components in the circuit. The operation may not be at the same point 
for the two different frequencies f, andf,. ‘The absolute frequency measurements 
made here enabled the effect of circuit changes to be studied in some detail. 
It was found that the modulating crystal diode used in the control circuit was 
particularly sensitive to frequency. Fora fixed setting of the cavity the frequency 
of control could be changed by as much as 1 part in 104 for extreme circuit 
conditions and direct measurements of refractive index by this method indicated 
that errors of 1 part in 10° (in m) could readily occur. It was decided therefore 
to connect the cavity in a separate circuit and to use the Pound circuit merely 
as a source of oscillations. 

A number of other very sensitive methods of using cavity resonators have 
been designed for detecting small changes of refractive index but they have not 
been used for absolute measurements with the required accuracy. 


§2. EXPERIMENTAL ARRANGEMENT 
The 24,000 Mc/s. (K band) oscillator in Figure 1 is controlled in a Pound 
circuit by a cavity resonator (not shown) which is provided with a micrometer 
adjustment for making small changes of frequency. The hybrid junction 
(‘magic T’) A is a part of this circuit and the usable power output is transmitted 
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through an attenuator towards a second hybrid junction B. At Bit is divided into | 
two equal parts, one going towards the measuring cavity resonator and the other * 
towards a load and matching stub. The fully absorbing load is equivalent to a 
perfectly matched cavity at resonance. The tuning stub can be made to represent + 
the mismatch of the cavity including the reflection from the window, which | 


Communication 
Receiver and 
Loudspeaker 
(with crystal 
calibrator) 


Frequency 
Control Unit 
Attenuator 


Pound | 
| 


Oscillator 


250 K-band 
i Mc/s. Oscillator 
K-band c.24,000 Mc/s | 
Attenuator Waveguide 
Compensating Arm Gas! 
ZZ, Measuring 
Cavity 
5 Mc/s. Comparison 
with 100 k¢/s 
Quartz Crystal beh 
clock 


Galvanometer 
Figure 1. Diagram of cavity bridge and frequency comparator. 


isolates it from the waveguide circuit and enables it to be evacuated or filled f 
with the gas under test. The waves returning from the stub and from the cavity # 
are again divided at the junction, a part returning towards the oscillator and a # 
part entering the detector arm. Hence the output from the crystal detector is | 
proportional to the vector sum of the two reflected waves. The matching stub 7 
is first fully withdrawn and the cavity tuned by adjustment of a small plunger () 
attached to a micrometer to give a minimum deflection; the penetration of the 7 
stub is then gradually increased and its position adjusted until the minimum } 
obtained by tuning the cavity is very near to zero. This corresponds to the true 2 
resonant condition of the cavity. Other conditions for zero deflection can be ‘ 
obtained with further penetration of the stub but the true resonant condition is < 
likely to give the greatest stability and reproducibility of setting. The stub is: 
situated at nearly the same electrical distance from the junction as the cavity) 
(owing to the properties of the hybrid junction the distances should in fact‘) 
differ by 4/8 for the zero current condition) and the balance is therefore obtained!) 
at the true resonant setting of the cavity for a considerable range of frequencies. | 

The waveguide bridge arrangement enables a very high sensitivity to bec 
obtained. The power entering junction B is such that with the cavity detuned! 
a galvanometer deflection equivalent to 10,000 scale divisions is obtained and aj 
deviation from the minimum of 0-1 division can be detected. ! 

A small amount of radiation, which enters the fourth arm of junction A due(. 
to a slight mismatch of the padding attenuators, is used for the frequency) 
measurement. It is mixed in a crystal diode with an output of 3,000 Mc/s.) 
obtained as shown from a harmonic of a 250 Mc/s. output derived by multiplication} 
from a 5 Mc/s. quartz oscillator. The frequency of the quartz oscillator is set!) 


e 
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exactly at its nominal value by reference to the N.P.L. standard quartz clock. 
The power at 250 Mc/s. is about 2 watts and this is applied to a crystal diode 
(a selected CV 247) which is connected across a tunable length of waveguide. 
The harmonic at 3,000 Mc/s. is selected and taken from the guide through a short 
length of coaxial cable to a second crystal diode (a selected VX 4026) connected 
across the waveguide carrying the small output from the K-band oscillator. 
This crystal acts as both harmonic producer and mixer. The heterodyne beat 
at about 9 Mc/s. produced between the eighth harmonic of the 3,000 Mc/s. 
output and the K-band oscillator is detected in a sensitive communication 
receiver. ‘The receiver possesses an internal quartz calibrator which can be 
used for the measurement of the beat frequency or for the setting of this 
frequency to an exact multiple of 0-1 Mc/s. by adjustment of the K-band oscillator. 

If the frequency multiplying stages are carefully designed to be free from 
phase modulation the purity of the beat note heard in the receiver gives an 
indication of the degree of stabilization attained by the Pound controlled circuit. 
It was found that the note could be set to below 200 c/s.—or better than 1 x 10-® 
of the operating frequency—and that any gradual drift, due mainly to the 
temperature change of the cavity in the Pound circuit, could be readily followed. 
It was thus possible to maintain the frequency of the oscillator to within 1 x 10-° 
of the desired value during a set of measurements. 


§3. DESIGN AND CALIBRATION OF THE CAVITY 
- One of the two cavities used for measuring the gas is shown in Figure 2. 
They are made from invar, the inner walls being silver or copper plated. ‘The 
Hy,; mode of resonance is used in the larger cavity and a groove is cut in the 
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Figure 2. Figure 3. Resonance curve. 


top wall to shift the frequency of the E,,, mode, which, in a perfect cylinder 
coincides with the Hy, mode, outside the working range of the instrument. 
A cavity of fixed dimensions could theoretically be used, the setting to resonance 
being made by adjustment of the oscillator frequency, but there are a number 
of practical advantages in making the cavity tunable over a small frequency 
range. If the setting is made by varying the oscillator frequency there is always 
PROC. PHYS. SOC. LXIV, I10—B 3L 
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a danger that the power and phase of the oscillations will also be changed, thereb ) 
reducing the accuracy of setting. Moreover, in a cavity of this size, there arg 
many possible modes of resonance and it is useful to be able to vary the resonant) 
frequency through a small range to make sure that none of these other mode ss 
falls within the frequency range used in the experiments. As will be seen in tha 
next section, the tuning facility also simplified the measurement of frequenc 
change on evacuation. An additional advantage, which is not used in this wor 
but which would be valuable in some experimental arrangements, is that if the) 
tuning range is greater than the change on evacuation, this change can b i 
compensated by retuning the cavity and the measurements can be made at a) 
fixed frequency. The value of f,—f, is then read from the calibration curve of 
the instrument. If this procedure is used, no particular care need be taken i 
ensure that the components of the circuit in which the cavity is connected are 
frequency insensitive and satisfactory operation would, for example, be obtained( 
in the arrangement used by Crain (1948). | 

An axial plunger was therefore fitted as shown in Figure 2. The cross bar ati 
the top of the micrometer is gripped by a spring clamp attached to a rotatingy) 
sealed joint when the cavity is mounted in its vacuum chamber and thus enables 
the adjustment to be made from the outside. : 

The two cavities used, designated by the letters A and B, have the dimensions: 
given in Table 1. 


Table 1. Details of Cavities 


Diameter of Diameter 


plunger of plunger ee plungeg 


, i h 
Cavity Mace Diameter Lengt 


(cm.) (cm.) material movement 1 

hole (cm.) (cm.) of 2 | 
A Ellente 6-660 4-491 js 1-09 silver 18 Mc/s. 
B leling 4-125 4-030 0-63 0-31 copper 9 Mc/s. 


It was found that with both cavities the setting to a minimum could be made? 
with a precision of a few kilocycles per second as may be judged by the bottom 
portion of a typical resonance curve reproduced in Figure 3. The micrometer) 
readings for the resonant condition were taken at intervals of 50kc/s. and the’ 
complete calibration curve for cavity A is reproduced in Figure 4. This curve: 
corresponds to a plunger movement of nearly one-half wavelength in the cavity”) 
and would have continued in a similar manner for a further penetration of the: 
plunger. However, the nearly straight part of the curve gave the required range) 
and sensitivity for measurements on air. The calibration was made with the: 
cavity filled with air, evacuated, and with different settings of the tuning stub!) 
in the waveguide bridge; and it was found that although the absolute frequency’ 
of resonance was, of course, modified, the law of frequency change apaifistl 
plunger position was the same in each case within 1 ke/s. for the complete range} 
of 18 Me/s. | 

The method of attaching the cavity to the vacuum system and of sealin 
cavity A is shown in Figure 5. It will be noticed that any bending of the babe 
plate of the containing chamber upon evacuation will not be communicated to 
the base of the cavity. Cavity B had a mica window 0-003 in. thick in place of 
the distrene plug. ‘Tests showed that the window was pushed in at the centre 
by 0-0015 in. on evacuation and, although the effect of this appeared from roug 

| 
| 
| 
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calculations to be negligible, it was decided that the solid plug one wavelength 
thick was to be preferred. Most of the measurements were therefore made with 
cavity A, which was also the more convenient because it happened to be an 
almost perfect match to the waveguide. Consequently, the matching stub was 
nearly fully withdrawn. It was, however, thought to be worth while checking 
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the results with cavity B, which differed in size, surface finish, and calibration. 
Very close agreement was obtained between the results thus helping to establish 
that some of the possible sources of error, discussed in more detail in § 5, had a 
negligible effect. 

§4. PROCEDURE OF OPERATION 

The method used was essentially the measurement of the frequency change 
produced on evacuating the cavity, but the following procedure was found to 
be most convenient in practice. The frequency of the oscillator was set 
approximately to the resonant frequency of the full cavity and was then adjusted 
slightly so that the beat frequency in the receiver coincided with the nearest 
harmonic of the 100 kc/s. calibration oscillator incorporated init. (This 100 ke/s. 
oscillator was itself set to give a zero beat with the N.P.L. standard.) ‘The cavity 
was then adjusted to resonance by a small movement of the tuning plunger. 
The gas was removed and the process repeated. The frequency change was thus 
obtained to the nearest 100kc/s. from the calibration of the receiver dial, the 
residual change being obtained from the difference in the plunger settings, and 
the calibration curve. In each case, the mean of six settings of the tuning plunger 
was taken and experience showed that the repetition accuracy of setting was then 
1kc/s. or 1 part in 24 x 10%. Since only a small plunger movement was required 
it was possible to work on a strictly linear part of the calibration, and thus avoid 
the need of referring to the actual curve. 

The cavity temperature was measured by a thermometer reading to + 0-01°c., 
and for the gases other than water vapour the pressures were obtained from a 
barometer reading to 0:02mm. Hg. Both of these instruments were calibrated 
shortly before and after the experiments against standards maintained at the 
National Physical Laboratory. The water vapour was admitted to the evacuated 
cavity from a vessel containing distilled water at a temperature at least 4° c. 

3 L-2 


868 L. Essen and K. D. Froome 
below that of the cavity. Its pressure was measured directly on a mercury j 
manometer made from glass tubing having an internal diameter ofs 175 cm, | 
the height of the mercury meniscus being read to 0:005mm. on a travelling 4 
microscope. The gases were dried by being passed through tubes containing 2 
glass wool impregnated with phosphorus pentoxide. In preliminary experiments ¢ 


Table 2. Summary of Measurements on Dry Gases at 0°c., 760mm. Hg 
Frequency of 
Gas (n—1)10® (e—1) 10° Silko Bie 0th Author Date 
AIR 
Ted 287:8+0-1 575-7+0-2 * Barrellt 1951 
576 0:5 Mc/s. Watson 1934 
567:0+1:0 1:0 Me/s. Hector & Woernley 1946 © 
572 9000 Mc/s. Crain 1948 | 
577:0+1:0 9000 Mc/s. Lyons, Birnbaum & Kryder 1948 | 
608 3000 Mc/s. Philips 1950 
575:-44+1:-4 9000 Mc/s. Birnbaum, Kryder& Lyons 1951 | 
288:15+0-10 576:0+0:2 24000 Mc/s. Essen & Froome 1951 
NITROGEN 
290°6 581:3 = Scheelt 1907 
587 0:5 Mc/s. Watson 1934 
57/9)-6 ta Oe One Se Hector & Woernley 1946 
587+ 2:0 Lyons, Birnbaum & Kryder 1948 
586-9-4.2-9 [one MCs “Bicnbaum, Keyder 6 Lone 
294:1+0:1 588:34+0:2 24000 Mc/s. Esscn & Froome 1951 
OXYGEN 
266°3 532°7 “ Rentschlert 1907 
531 0:5 Mc/s Watson 1934 
5233+ 1-0-0 Mels Hector & Woernley 1946 © 
SoZ Deel Ou lc OMVIC/S Jelatis 1948 @ 
SSL st LAD Lyons, Birnbaum & Kryder 1948 © 
S30-0i4-90f so, IG *-\pienbaune Kryder & Lyons 1951 — 
266:4+0:2 531:0+0-4 24000 Mc/s. Essen & Froome 1951 
ARGON 
279:-2 558°5 . Burtont 1907 — 
277°3 554-7 Cuthbertson & Cuthbertson 1932 
550 0:5 Mc/s Watson 1934 7 
545:140°5 1:0 Mc/s Hector & Woernley 1946 7 
554:2+0:9 1-0 Me/s. Jelatis 1948 
277°8+0:2 555:74+0:4 24000 Mc/s. Essen & Froome 1951 | 
| 
CarRBON DIOXIDE 
989 0:5 Mc/s. Watson 1934 
987:5+2:0 1:0 Mc/s. Hector & Woernley 1946 
988+ 2-0 Lyons, Birnbaum & Kryder 1948 | 
reese 60) poo Orleis: Birnbaum, Kryder & Lyons 1951 | 
4944+1:0 98842:-0 24000 Mc/s. Essen & Froome 1951 | 


* Optical determination, extrapolated to infinite wavelength. 


t Quoted by Kaye and Laby (1948). 


t From a revision of the most accurate determinations, 


Note. Our values of (e—1) are obtained from equation (1) on the assumption that the value 
of (%7—1)10° are 0-4 for air, 1-9 for oxygen and zero for the other gases and for water vapour. The 
results of other workers are given as published by them although in the microwave determinationg 
the quantity actually measured will have been refractive index and not dielectric constant. 


permeability has been allowed for in the values given by Birnbaum, Kryder and Lyons (1951)i 
but possibly not by the other workers. 
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silica gel was used to remove most of the moisture, thus lengthening the life of the 
phosphorus pentoxide, but there was some evidence that the gel first absorbed 
and later emitted other gases, especially CO,, and thus affected the purity of the 
gas under test. It was also feared that such absorption might be differential for 
the constituents of dry air; and although there was no direct evidence of this it is 
perhaps significant that the standard deviation of the results was reduced by 
one-half after the gel had been removed from the system. At first higher values 
of refractive index were obtained when the cavity was filled slowly than when it 
was filled quickly. ‘This was found to be due to moisture given off from the walls. 
of the Pyrex glass cylinder surrounding the cavity or from the glass tubing in the 
system. When a tray of pentoxide was placed inside the Pyrex container and 
near the cavity the same value was obtained from measurements made with 
rapid and slow filling and from measurements made a few minutes and 
twelve hours after filling. Carbon dioxide was removed from the gases by passing 
them through a tube containing soda lime. 


S52 RESWELS SN DeACCURACY, 

The results of the measurements which have been briefly announced (Essen 
and Froome 1951) are given in Tables 2 and 3, together with those of a number of 
other workers, which will be discussed in $7. 

Table 3. Values of m for Unsaturated Water Vapour at 20°c., 10mm. Hg 


Dipole moment 


(n—1) 10° Sangisee Frequency Author Date 

2-94 z Barrell & Sears 1939 
62-7 1-842 + 0-008 1 Mc/s. Sanger 1930 
61:3+0-4 1:831+ 0-006 0-5 Mc/s. Stranathan 1935 
61-3 1-84 9000 Mc/s. Crain 1948 
62-4 a= 3000 Mc/s. Philips 1950 
60°7+0-1 1-839 + 0-002 24000 Mc/s. Essen & Froome 1951 


* Optical determination, extrapolated to infinite wavelength. 


The air, nitrogen, oxygen and argon were free from carbon dioxide 
and the refractive indices are given for 0°c. and 760mm. Hg, the pressure 
being measured in millimetres of mercury at 0°c. and standard gravity 
(760 mm. Hg = 1,013,250 dyne/cm*?). The measurements were made at a mean 
temperature of 20°c. and a mean pressure of 760 mm. Hg, and the results have 
been reduced to the standard condition by the extrapolation formula (2), given 
later in the paper. ‘The measurements on water vapour were made at 
temperatures between 15°c. and 25°c., and with pressures between 6mm. Hg 
and 14mm. Hg, the result given being obtained from formula (8) which is based 
on the mean of the experimental values. The standard conditions have in this 
case been taken as those normally pertaining in a room because water vapour 
is a polar gas and its extrapolation formula cannot be given with the same 
confidence as that of the other gases. 

A correction of +3-0 x 10-7 has been added to all the measured values except . 
that for water vapour to allow for the slight increase in the dimensions of the 
cavity when the container is evacuated, due to the hydrostatic expansion of the 
invar walls and end plates. 

For the gases other than air and water vapour, ‘tank’ gases were used. 
It was found essential to flush the equipment at least eight times when changing 
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from one gas to another, before taking measurements; and spectroscopically f 


pure gases are not available in large enough containers for this procedure. 
The limits of error given in the tables are derived as follows: 


(a) Air. 
Random error: standard deviation of experimental results 
when reduced to N.T.P. +4x1I03 


Estimated maximum systematic errors: 
thermometer calibration +2 10-8; barometer calibration +410 
Possible phase shift errors inherent in circuit +4x104 
Invar expansion correction error +3 x104 
Statistical total = +8 x10-* 


Consequently we have adopted an error of +1X10~’, and this should be 
regarded as being equivalent to a standard deviation. 

(b) Nitrogen was stated by the suppliers to have less than 1 x 10-® impurity 
and the result for this gas is therefore as accurate as that for dry air. 

(c) Oxygen was stated to have 0-5°% nitrogen as an impurity and wider limits 
have therefore been given to allow for an error of about 100° in the estimation 
of this impurity. ; 

(d) Argon was stated to have 0-2°% air as an impurity, so wider limits are again 
given. 

(e) Carbon dioxide was stated to have 0-5°% air as an impurity and for this 
gas the standard deviation was also higher than for the other gases, being about 
+2x10-’%. Wider limits are therefore given. 

(f) For water vapour there is no thermometer error, or invar correction to 
take into account. Neither is there a systematic pressure error since variations 
in the setting on the manometer meniscus are included in the experimental 
spread of the results. The error given is therefore completely represented by the 
standard deviation. Possible effects due to adsorbed layers of moisture on the 
walls of the cavity are considered in §7 and it is concluded that they had no 
appreciable effect on our results. 

No mention has yet been made of possible errors due to imperfections in the 
hybrid junction and to spurious reflections from the waveguide arms and joints 
constituting the cavity bridge. It has been tacitly assumed that this system 
is perfect and that the crystal detector current (Figure 1) is simply the vector 
sum of the wave reflected from the cavity and that from the equally spaced 
matching stub in the compensating arm of the bridge. In addition to these two 
components, however, there may be others due to the direct coupling from the 
input to the detector arm of the junction and to various mismatches along the 
cavity arm, and along the compensating arm of the bridge. Thus it is seen that 
the minimum galvanometer deflection corresponding to the resonant setting of 
the gas measuring cavity is obtained when the vector sum of these five components 
is zero, ‘Uhe spurious components will not cause any error so long as they remain 
constant; but since they arise from reflections which are at different distances 
from the detector they are likely to alter in phase, relative to the wave from the 
cavity, when the frequency is changed. The minimum galvanometer deflection 
will then be obtained at a setting of the cavity which is slightly different from the 
true resonant condition. 

The magnitude of this effect was checked for cavity A, which was so well 
matched at resonance that the compensating stub could be fully withdrawn. 
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In this case, therefore, there were only four components in the detector arm—the 
wanted one from the cavity and three unwanted ones. The phase of the reflection 
from the cavity relative to the others was altered either by inserting near the 
cavity a short spacer or by leaving there a small air gap. It was estimated that the 
phase shift produced by these means corresponded to that which would be 
produced by a frequency change of 200 Mc/s. but there was only a 10-20kc/s. 
change in the resonant setting of the cavity. Since in practice the frequency 
changes were only about 7 Mc/s., it was safe to assume that errors arising from 
this cause were almost entirely negligible. The allowance included in the 
estimation of errors is certainly more than adequate. 


; §6. EXTRAPOLATION FORMULA 
In practice the value of the refractive index of air is required at the prevailing 
conditions of temperature, pressure and humidity, and the results are therefore 
most usefully expressed in the form of an extrapolation formula. It is sufficient 
for most purposes to assume that the refractive index varies proportionally with 
density, that dry air obeys the ideal gas laws and water vapour a simple Debye 


equation. But if the values obtained are to be accurate to within 1 x 10-7 over 
a wide range of conditions, more accurate formulae must be used to allow for 


observed deviations from these simple laws. Dry air, carbon dioxide and water 


- vapour will first be treated separately. 


6.1. Dry Air 

The equation used for reducing our results at a pressure pmm.Hg and 

temperature ¢°c. to standard conditions was 
Be ie = =, 760 LP ; -6 

(n,,—1)= 760°606 (1 +0-0036610) [1+ (1:049—0-01572)107%p] .. 1. (2) 

It was found by Barrell and Sears (1939) from optical refractivity measurements 
on dry air, and it holds good to a very high degree of accuracy over a wide range 
of temperature and pressure. It is expected to hold with almost equal accuracy 
for nitrogen, oxygen and argon. 


6.2. Carbon Dioxide 
Carbon dioxide may be assumed to obey the ideal gas laws to the accuracy 
of our results for this gas, and the equation is simply 
ie thi eee (3) 
where T is the absolute temperature. Absolute zero is assumed to be —273°c. 
6.3. Water Vapour 


Water vapour is a polar gas with an electric dipole moment and therefore 
obeys Debye’s (1929) equation of the form 


<-1=9' (74 7) ie (4) 


where p’ is the pressure that the vapour exerts assuming the ‘ideal gas’ laws 
to be obeyed at any fixed temperature. The A term represents the contribution 
to the dielectric constant of the atomic and electronic polarization, and the 
B term the contribution of the dipole moment. Our measurements were made 
over only a small temperature range and do not therefore enable the two constants 
to be accurately determined. Other measurements made at radio frequencies 
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are scarcely accurate enough to give a reliable value for A because this te } 
is only about one-twentieth of the B term at these frequencies. It was therefore: 
decided to use the values of A obtained by extrapolation of the optical dispersioni 


formula to infinite wavelengths. > |i 
For p=10mm. Hg and T7=293°x., Barrell and Sears (1939) give n—1 fon 


water vapour as 2-944 x 10-6 or 1 


e—1=5°888 x 10*=p Aj293. ee (5) 
It is also shown in their paper that” | 
p =p(L+24x%10%). sss 2 ort eee (6) | 


This formula is an averaged one; but it is accurate over the range of conditions: 

we stipulate later; thus from (6) and (5) we obtain 
A=1725 x10) 8 ee (7) 
This value is believed to be accurate to 0-005. Using it in (4) for each of| 
our measured values of «—1, values of B are calculated. The average value of 
B was 0:9913. Putting these values of A and B in (4), we obtain from (4) and (6), ; 
after some simplification, the working formula: 


(n—1) 10° = ore (1+ “T ) (+24 x 10-4). Hep: (8) 


The mean measured result given in Table 2 is obtained by putting T=293} 
and p=10 in this equation. It may be noted that any error in the value assumed | 
for A is largely compensated by the different value of B that is obtained as a result. # 
Thus A could be in error by 40° which is inconceivable, without affecting the 7 
validity of the equation within the required accuracy for temperatures up te jj 
60° c. and unsaturated vapour pressures up to 100mm. Hg. The main interest 
in using the most accurate obtainable value of A is that, from the value of B, . 
it is possible to calculate the dipole moment of the water molecule. 


6.4. The Dipole Moment of the Water Molecule 
Using Debye’s theoretical expression for the B term in (4) and the ideal | 
gas laws, we obtain, converting the pressure p to dynes per cm?, 


_ Arp? 101325 108 9 
= Rape (9) | 
where k is Boltzmann’s constant (1-3805 + 0:0003 x 10-16 erg deg) and , is the + 
dipole moment. 


With B=0-9913 we obtain 
pe =(1-8395.0-002)10- Bc. Sa eee (10) 


To change this value by 0-002 would need an error of 4°, in the A term of f 
Debye’s equation compared with the estimated uncertainty of 0-3%. | 


- —— 


| 


6.5. The Extrapolation Equation for Moist Air 1 
. Combining equations (2) and (3) and (8), and inserting the values from | 
lables 2 and 3, the refractive index of air at a temperature ¢°c. and atmospheric 
pressure p becomes 


0-37884p, Z 1-749, 
73.0-00366I7 1 + (1-049 —0-01572)10-*p,] + — 


273 +t 
8624p; 5748 ; 
27st (1 273 =) (1-5 2:4.% 1062s) 5 aan (11) 


(n, »—1)10°= 
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where /, is the partial pressure of dry, CO,-free air, p, is the partial pressure 
of COs, ps is the partial pressure of water vapour and p=p,+p.+p3. All pressures 
are expressed in millimetres of mercury at 0°c. at standard gravity. 

Equation (11) can be used for all atmospheric conditions with temperatures 
between —20°c. and +60°c. and unsaturated water vapour pressures less than 
about 100mm. Hg, without introducing additional errors in m as great as that 
of the actual experimental results. It may also be expected to hold for all radio 
wavelengths above 7mm. where the strong oxygen absorption line at 5mm. 
may begin to have an effect. Water vapour has a weak absorption at 1:35cm. 
but it appears from Van Vleck’s (1942) theoretical work that this should have 
an entirely negligible effect upon refractive index. The equation (8) for water 
vapour should in fact hold for wavelengths down to 3mm. when the strong 
absorption in the region of 1mm. becomes effective. The absorption in this 
region is over 1,000 times stronger than that at 1:35cm. 

For most practical work, such as radar applications, where the very highest 
precision is not required, equation (11) can be simplified by assuming that dry 
air behaves as an ideal gas and that water vapour behaves as ideal at any one 
temperature. Adjusting the constants to give the correct measured values, 
equation (11) becomes 


103-49 177- : 
(m »—1)108 = eae ‘ vee A (1 + Tr) es ete (12) 
where 7 =273+1. 

This equation is still fully accurate at the temperature and pressure of the 
measurements but introduces an error in m amounting to about 5 x 10-7 at the 
extremes of the specified range of temperatures and with normal water vapour 
pressures. For most purposes the effect of the presence of CO, can be neglected, 
since the percentage present in open air is only 0-03 by volume, but ‘indoors’ 
it may be several times this amount. Putting p,=0, we obtain the familiar 
form of equation. ‘The values are about 4x 10° lower than those obtained 
from the equation given by Smith-Rose and Stickland (1943), for normal 
atmospheric conditions, and 15 x 10-® lower than those obtained from another 
commonly used formula. * 

The use of these extrapolation formulae involves the measurement of the 
water vapour pressure of the atmosphere. Some tests were made to check 
the accuracy with which this could be done by using the ordinary electrically 
operated ‘wet and dry bulb’ hygrometer. In these tests it was necessary to 
avoid the rapid absorption of moisture by the glass walls of the equipment, and 
rubber tubing was therefore connected directly to holes in the cavity resonator 
and a continuous slow-moving stream of air taken from the vicinity of the 
hygrometer was drawn through it until a few seconds before the measurements 
were made. It was estimated that the accuracy of measuring the vapour pressure 
was 0:05 mm. Hg, giving, in a total pressure of 10mm., an uncertainty of 0:5% 
in the water vapour contribution to refractive index. ‘This corresponds to 
3 x 10-7 in the total refractive index; and it was found that the measured values 
agreed with the calculated values to this accuracy. ‘Thus if precautions are taken 
to prevent the absorption of the vapour by the walls of the equipment, it is rather 
more accurate to measure the refractive index of the actual sample than to 
calculate it from humidity measurements made with an Assmann hygrometer. 


* See Foreword to Meteorological Factors in Radio-Wave Propagation (London: Physical Society ) 
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$7. DISCUSSION, OF RUE RESULTS 


There have been many previous measurements of the refractive index and 
dielectric constant of dry air by optical and electrical methods, and there is no 
need to discuss them all here. The optical values tabulated by Bender (1938) 
show a total spread in m of 2 10-6 but quite recently, since the completion of 
the work described in this paper, Barrell (1951) has derived a new extrapolated 
value based on the most recent optical measurements and this value only is 
given in Table 2. Hector and Woernley (1946) give a list of values of dielectric 
constants obtained by electrical methods, and the spread of these values is 6 x 10™. 
Many of these earlier results are clearly unreliable and only those of Watson are 
included in the Table. Hector and Woernley’s own results seem to have a rather 
large systematic error. 

It will be seen that Barrell’s extrapolated value of m for air is about 4 x 107? 
lower than our measured value. The difference is in part due to the magnetic 
permeability of oxygen, the effect of which is removed in the derived values of 
dielectric constant. The discrepancy is then only 3x10~? (equivalent to 
1-5 x 10-7 for m) and this is within the combined experimental error. ‘The radio- 
frequency value should in any case be slightly greater than the extrapolated 
optical value, even when allowance has been made for permeability because of | 
the magnetic dipole resonant absorption by oxygen at wavelengths near 5mm. 
Van Vleck estimates that the effect of this on the dielectric constant of oxygen is 
to increase it by 2 x 10~’ giving an increase of 4 x 10-8 in the value for air. 

In the extrapolation formula for dry air it has been assumed that the 
contribution due to permeability varies with temperature in the same way as 
that due to dielectric constant because, although in fact there will also be a 
1/T? term, neglecting this will introduce an error of only 2 x 10-5 at the extremes 
of temperature given. The error for pure oxygen will, of course, be five times as 
great, but in either case the errors will be negligible near 20° c. 

Van Vleck also estimates that the effect of the water vapour absorption at 
1-35 mm. is 1 x 10-4 of the total polarization. The actual effect could be much 
greater than this before it would invalidate the use of the extrapolation formula 
over the stated frequency range. 

In the radio-frequency range our results agree with those of Watson (1934) 
obtained at 0-5 Mc/s. within the accuracy of his measurements for all the gases 
except argon. ‘They also agree well with those of Lyons, Birnbaum and 
Kryder (1948) and Birnbaum, Kryder and Lyons (1951) obtained by a cavity 
method at 9,000 Mc/s. and with those of Jelatis (1948) at 1 Me/s. for oxygen 
and argon. 

The fairly close agreement between the various values for water vapour is 
satisfactory in view of the extra difficulties met with in the measurement of 
a vapour. Stranathan (1935), who used a frequency change method at 0-5 Mc/s. 
discusses the difficulties due to adsorption as the vapour pressure approaches 
saturation. ‘The refractive index first increases linearly with pressure but 
then much more rapidly, due, it is assumed, to the sudden formation of a 
layer of water 100 molecules thick. But Newbound (1949), working in the 
optical region, states that the thickness of the molecular layer is proportional to 
vapour pressure and is 100 molecules thick at a pressure of 10mm. Hg and a 
temperature of 20°c. Newbound’s result is not substantiated by the optical 
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work of Barrell and Sears (1939) and Cuthbertson and Cuthbertson (1932). 
In the cavity apparatus the effect of a layer would be much less than in the case 
of a condenser with a small gap between the plates such as used by Stranathan 
(1935). ‘The cavity results are therefore likely to be the more reliable and the 
good agreement obtained suggests that the effect of adsorbed layers, if any, must 
be very small, and that the values given by the present investigation are accurate 
within the limits stated. 
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ABSTRACT. ‘The measurement of the alpha-particle activity of a ‘thick’ radioactive 
source containing the natural radioactive series is discussed critically. | 

Two methods of determining the thorium—uranium ratio are considered for use in ore or 
mineral assays. The energy discrimination method is simple in procedure and can take 
account of lack of secular equilibrium between members of a series in the radioactive | 
sample. The beta—alpha coincidence method, making use of the difference in the disinte- 
gration constants of radium C’ and thorium C’, is more complicated. 

Initial experimental verification of the theory has been obtained. The estimated | 
statistical error in determining the thorium—uranium ratio by the energy discrimination 
method is +5% in the region of Th/U=1, with a 4% sample (U3O0, equivalent concen- 
tration) in 4 minutes counting. A similar accuracy is obtainable in the coincidence method 
with a counting period 10 times greater. : 


§1. INTRODUCTION 
INCE the range of alpha-particles in naturally occurring radioactive materials 
is extremely small, of the order of 30 microns, it has been customary te 
make measurements of activity and alpha-particle range using thin and 
finely powered samples. In this condition the effect of self-absorption in the sample — 
is thereby made as small as is practically possible so that the collection of the | 
particles by the chamber, fluorescent screen or photographic plate is dependent — 
mainly upon the geometry of the system. Reduction of self-absorption also serves | 
to increase the discrimination of alpha-particle range measurements. 
The application of thin-source technique to the routine assay of radioactive © 
minerals or ores raises two serious objections to the method: (i) use of very thin | 
sources implies low count rates, which for samples of low concentration will often / 
be less than the detector background rate; (ii) the preparation of very thin uniform | 
samples of complex ores, usually in powdered form, is in practice very difficult. ; 
In view of these objections the activity of the ‘ thick’ source has been investi- | 
gated. A thick source of alpha-particles is defined as one of thickness greater | 
than the internal range of the most energetic particle emitted. 


§2. ALPHA-PARTICLE ASSAY 

The alpha-particle activity of a radioactive source has been calculated by Finney » 
and Evans (1935). ‘The main steps in the calculation are repeated in order to | 
indicate restrictions imposed by the dependence upon particle energy of the stopping | 
power otf the source and, also, to extend the argument to include some cases of | 
disequilibrium in a radioactive series. The practical effect of self-absorption and | 
disequilibrium upon the measured activity of a source is discussed. | 


2.1. Calculation of Alpha-Particle Activity 


In Figure 1, dA represents an elementary portion of collecting area at a distance }} 
d from the surface of a semi-infinite radioactive source containing a mono-energetic : 
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alpha-particle emitter. The area dA will receive only particles emitted by that 
part of the source within a total distance equivalent to the range in air R, of the 
alpha-particle. 


Figure 1. Collection geometry for thick source. 
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Then the number of particles per second incident upon dA from an elementary 
volume dV at (r, ¢, @) is given by (cos 6/47r?) dA p, NdV, where N is the intrinsic 
activity of source in particles/gm. sec. and p, is the density of the source. 

The number of particles per second (count rate) incident upon dA from the 
whole source is given by 


dc 


a ee | tie ba (1) 


An Jy r2 
It can be shown after integration that the total count rate over the whole 
collecting area A is given by 


c=(NAR,p,/4R,) (Ree aye —D), (kr), 


where R, is the range of alpha-particles in the source and D the range in air corre- 
sponding to the minimum energy detectable by the collector and associated 


apparatus. 
Extension of the calculation to the case of an ore containing members of a 


tadioactive series emitting particles of air range R,,, Ry». . R,,, gives for the rate of 
collection of particles from the nth member 


< N,A (Rig Dd) 
CREE OT: 


D bs 
where S=R,/R,p, is defined as the mass stopping power of the source. 


an 
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For the whole series 
Nn 


(Aa 


Removal of the factor S from within the summation sign is strictly justified only | 


if S is independent of alpha-particle energy. Also, in the integration of (1), 


the upper limit for 7 involves the ratio R,/R,(=1/pyS) which is assumed to be 


independent of the residual energy of the alpha-particle as it emerges from the 


surface of the source. Beharrel (1949) has found that calculation of S by the | 
Bragg—Kleeman (1905) rule, which takes no account of energy dependence, gives. | 


results that are 10° too low for the common rock-forming minerals. 
If the series is in equilibrium 


Ni=N,=. = N= (2) 
n NA (Ra—D—d)? NA | 
and Ce ao = 4S 2 ay ETS hoa = 4S F(D, a). ee eeee (3) 


The latter general form (3) can apply to source-collector geometries other than — 


the specific case considered above. 
In general, a radioactive ore will contain members of both the uranium and 


thorium series in intimate admixture. The rate of particle collection from such 
a source is given by the sum of two expressions derived from (3). Provided that 
. the uranium and thorium are closely associated or that the respective particles 
are otherwise subjected to the same stopping power, then 


NyUA NynThA Fn,(D, @) 


C= Cut Cm = ae Fy (D, d) + 4S 


A 
ae [Ny Fy(D,d)U+ No Fm(D,d)Th), ...... (4) 


where N,,, Nip, are the specific activities per gramme per second of uranium and — 
thorium and U, Th are the concentrations in gm. per gm. of uranium and thorium 7 


in the source. 


Theoretical values of NF(D,d) for each series are plotted in Figure 2 asa | 
function of d for the idealized case D=0. A family of curves may be so con- — 


structed for varying D. 
2.2. Practical Considerations 

2.2.1. Self-absorption. It is apparent from equation (4) that a straightforward 
alpha-particle assay of a complex radioactive source is dependent not only upon 
the relative concentrations of uranium and thorium but also upon the particle 
collection geometry and mass stopping power, i.e. self-absorption. Knowledge 
of the variation of the last factor with density, atomic weight and grain size is 
incomplete. A direct method of alpha-particle assay, even for ores containing 
the uranium or thorium series only, therefore proves unreliable except: as an 
activity comparison method for range of ore samples of similar constitution. 


Attempts have been made by Finney and Evans (1935) and Nogami and Hurley | 


(1948) to assess the mass stopping power for various types of radioactive ore. 


These are based upon the absorption rule of Bragg and Kleeman (1905) which | 


states that the product of mass stopping power and the square root of the atomic 
weight is constant. Application of this rule implies a knowledge of the con- 
stituents, in their relative proportions, of the ore sample and would be impractic- 
able in the routine assay of ore samples. 
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It is desirable that the accuracy of radioactivity assay should be independent 
of the variation of such absorption factors. 

2.2.2. Diffusion. ‘The alpha particle activity of a complex source as calculated 
in §2.1 may be modified by lack of radioactive equilibrium. Owing to the loss 
by diffusion of the gaseous members (radon, actinon and thoron) of the three 
radioactive series a different state of equilibrium is reached when the net rate of 
formation of each is balanced by the loss by diffusion through the surface. A 
severe physical disturbance such as occurs in grinding will upset equilibrium, 
creating a transient condition that develops into a state of practical equilibrium 
after several radioactive half-lives of the diffusing gas. It is estimated that the 
loss of radon and its products in this manner may cause an error of 10°% in the 
calculated activity of uranium source, although knowledge of the coefficient of 
diffusion of radon is uncertain. Diffusion of actinon and thoron produced little 
effect, owing to their much shorter half-lives. 

It is possible to demonstrate quantitatively the effect of diffusion upon the 
count rate—distance relation. In Figure 1 the depth of the elementary volume 
below the surface of the source is (rcos@—d). It can be shown that the relative 
deficiency of gas at this depth is given by exp {—(A/f)!(rcos@—d)} under the 
new equilibrium condition, A being the radioactive disintegration constant and 
f the coefficient of diffusion (per unit concentration gradient per unit area). The 
integral (1) is thus modified to. give for the diffusing member 


_ psNdA | cos @ 
J 


An r 


dc [1 —exp {—(A/f)!? (r cos 0—d)}] dV. 

This may be integrated as before: a correction factor is derived for application to 
the calculated count rate of the diffusing gas. Factors for the descendents of the 
diffusing member may be obtained by substituting the appropriate values for R,. 
Thus for a uranium ore the loss of activity of the radon descendents will be 
governed by the radon disintegration constant but dependent upon individual 
alpha-ranges. The activity calculated for the whole ore is made up of the 
unmodified contributions of members preceding radon in the series plus the radon 
and succeeding activities corrected for diffusion. The magnitude of the diffusion 
effect, for radon only, is illustrated in Figure 3 by plotting the correction factor 
against source—collector distance. 

2.2.3. Other disturbances of equilibrium. Chemical or physical changes in 
past geological time may have upset the natural equilibrium conditions of radio- 
active ores. For example, the leeching and dissolution of primary uranium 
minerals by circulating ground waters and the deposition of secondary compounds 
after transportation in solution may result in considerable disturbances in secular 
equilibrium. In the case of a secondary deposit with excess radium, provided 
sufficient time has elapsed to establish a modified state of equilibrium, the 
condition (2) is changed to 


where N, represents the activity due to the transported member, and provided 
that the time elapsed is much longer than the half-life governing build-up from the 
rth to nth members in the series. 
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Then (3) corrected for disequilibrium becomes 


cS - [NF(D;d) EN’ FAD, 2), scacQnaa (5) 
” (Ri, —D—d)? 
where FAD, da) oe 5 in : 


§3. ENERGY DISCRIMINATION METHOD OF MEASURING 
THORIUM-URANIUM RATIO (Th/U) 


3.1. Theory 


The influence upon the measured activities of the mass stopping power 


restricts the use of an alpha-particle assay as a measure of the total active content 


of an ore or mineral sample (cf. §2.2.1). For an estimation of the thorium—_ 


uranium ratio this difficulty may be removed. For convenience equation (4) may 


be written C=(A/4,S)[XU+YTh] where X=NyF,(D,d), Y=NmFr(D, @). 


By the choice of two values of either of the two parameters D or d 


A ; 
hs (XLU+Y, Thi (G6 


Cy X,U+Y,Th 6) 
Cai Xacti VoD oe onda 


A 
C= 45 [X,U+ Y,Th],; Cy= 


and hence 


Since for a source containing the uranium (7h =0) or thorium (U =0) series only, 
(C,/Co)y =X4/Xo and (C,/Co)an = Y1/Vo respectively, then equation (6) becomes, | 


after substitution, 


ie ACOs » uw TG <= eee 
(C/Co)an — Cy/Co Xo VY1—-X1 Yo Koki Ages U 

re 1 GG he, ei a 
Cy/Co—(Ci/Co)u = Xo¥i-XiYo Xo V¥1—X Vo Th 

At Cy/Co— (Ci/Co)u = Yo Th 


(Cy/Co)m—CyCy X, U" j 


The equations (7), which are similar to that obtained by Finney and Evans (1935), t 
represent alternative forms of the same relation, each is linear in Th/U (or U/ Th) | 
and use of a particular form depends upon the availability of thorium or uranium _ 


standards or both. 


Methods of measuring 7Th/U may be based upon this relation, thereby | 


discriminating between the alpha-particle range or energy characteristics of the 
thorium and uranium series. For very weak activities Finney and Evans (1935) 
used an ionization chamber as detector with absorbers of various thickness: 


Curie (1946) used a photographic plate and discrimination by individual esti- 
mation of path length for this purpose. In the present work, which is primarily _ 


directed towards devising a method suitable for routine ore analysis, the 


scintillation counter is used as the detector. Discrimination is introduced by | 


varying the minimum output pulse amplitude accepted by the counting circuit, 


1.e. discrimination effectively in terms of the parameter D, the minimum | 
detectable alpha-particle range. In practice this is achieved by simple adjustment | 
of a pulse-amplitude discriminator control and represents a considerable gain | 


in speed and convenience compared with the previous methods. 
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It should be noted that a determination of Th/U by the method outlined may 
only be attempted when the unknown sample is known to be in equilibrium. 
Alternatively, it is possible in a range of samples to measure disequilibrium of a 
given type provided 7h/U is constant.. Proceeding as before in the simple case 
of an uranium series sample with no thorium present (7h =0) the disequilibrium 
equation (5) gives three alternative equations similar to (7) of the form 


Cy/ Co = (Cr/Co)u aK M 

(Cy/Co)m — Cy/Co U 
where ™ is the concentration gm. per gm. of radium (say) excess or deficit 
compared with the natural equilibrium amount, (C,/C)), is the count ratio for 
a sample containing radium and its descendents only and K is a constant. 
Equations (7) and (8) may be interpreted qualitatively for uranium ore samples 
that are either out of equilibrium or contain thorium but not both. Thus, if the 
count ratio C,/C) is greater than that of a standard sample containing the uranium 
series known to be in equilibrium, then the unknown sample contains thorium or 
an excess of radium. If C,/C, is smaller than for the standard then the unknown 
has a radium deficit. 

In the general case of course, the unknown sample will be out of equilibrium 
and also contain thorium. It is possible that by taking a third value of either of 
the parameters D or d, the ratios Th/U and M/U might be determined uniquely 
in the general case. In the present work, however, measurements of disequili- 
brium and 7h/U are considered separately and simultaneous determinations of 
the two ratios in a single ore have not been attempted. 


omen eS PERE (3) 


3.2. Apparatus and Procedure 

The apparatus used in making T7h/U measurements by this method comprises 
a scintillation counter, linear amplifier, pulse-amplitude discriminator and 
scaling unit. 

The scintillation counter consists of a silver-activated zinc sulphide (Levy 
and West G86) screen and an RCA type 931A photo-multiplier operated at a 
final anode voltage of 1,000. ‘This is followed by a linear amplifier (AERE type 
1008) operated between the gain ranges 160 to 1,600, having a gain stability of 
about 0-5°%; the pulse-amplitude discriminator (AERE type 1028) is set in 
conjunction with the amplifier gain control. ‘The scaling unit used was a modified 

Cinema-Television type UC2A. 

The discriminator is employed as a means of effectively varying the minimum 
detectable alpha-particle range D. Direct conversion of discriminator readings 
in terms of D and the computation of the constants of (7) are undesirable and 
unnecessary. Such computation, as carried out by Finney and Evans (1935), 
involves the evaluation of the function F(D, d) for each yalue used of D and d and 
demands a precise knowledge of the geometry of the alpha-particle counter. 
The constants of (7) are best determined by calibration with sources of known 
Th/U. 

With collecting and source areas of 20cm? and a source-collector distance of 
2mm. the overall sensitivity of the apparatus is such that 1% uranium source 
will give about 1,500 counts per minute: under these conditions the photo- 
multiplier ‘noise’ count is about 100 per minute. Alternatively by reducing 
the amplifier gain so that the noise count becomes negligible the sensitivity is 
reduced to 500 counts per minute for a 1° uranium source. 
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3.3. Experimental Results and Discussion 
A series of ore samples was prepared from pitchblende (70% U3O,) and 


thorium oxide in varying proportions. ‘These samples, diluted with powdered | 
quartz to assay roughly 4% U;O, equivalent concentration provided a range of * 


Th/U values at a uniform density. ‘The count ratio functions are plotted against H 


Th/U in Figure 4, for two settings of the pulse-amplitude discriminator, and | 


demonstrate initial verification of equation (7). 
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Figure 4. Count ratio function plotted against thorium—uranium ratio. 


Another series of ore samples was prepared from separate pitchblende and 
monazite (9% ThO,, 0-35% U3;0,) sources diluted with varying proportions of — 
powdered lead and magnesium oxides to give a range of density variation at — 


values of Th/U of 0 and approximately 25. According to equation (6) the count — 


ratio C,/C, should be independent of density. With the ore samples used — 
C,/Cy was found to increase with decrease of bulk density at both Th/U values. 


In particular this was due to decrease in Co, as predicted by equation (3); but a | 


disproportionately smaller decrease in Cj. 

Closer inspection of these samples showed a range of grain size in the inactive 
materials of from 4 to 80 microns. ‘The layers in and near the surface, as in the 
bulk of the sample, consist of mixed active and inactive grains. Since the internal 
range of the most energetic alpha-particle is about 50 microns, particles from 
submerged active grains will be wholly or in some cases partly absorbed by 


adjacent inactive grains. ‘The partial absorption is determined by the average 
stopping power of the two types of inactive grain. Because of the range of grain 
sizes this absorption may act selectively such that the particles of longer range 
as measured by C; possibly suffer proportionately less absorption than those of 
shorter range which are included in the measurement of Cy. The decrease of 


bulk density is achieved by increasing the proportion of magnesium oxide, | 


thereby increasing the average mass stopping power of the inactive dilutant > 


grains. It is conceivable that the increase in mass stopping power will account 
for the observed variation of C,/Cy with density. 


This explanation has not yet been tested experimentally, but if valid, it follows 
that the minimum grain size in a heterogeneous ore sample should be greater | 
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than 50 microns, i.e. each grain, active and inactive, should represent a thick 
source with respect to the alpha-particles emitted within the sample. 

As a further check of a possible density effect a careful comparison was made of 
pitchblende (70% U Og) and davidite (9% U,0,) samples known to contain 
negligible amounts of thorium. The intrinsic specific gravities of these minerals. 
are respectively 8 and 4-5. The count ratios C,/C, of the two samples were 
identical over the practical range of pulse-amplitude discriminator setting, and 
hence independent of the markedly different densities and concentrations of the 
two samples. 

The count ratio C,/Cy was next shown to be independent, for all practical 
pulse-amplitude discriminator settings, of the grain size in samples of davidite 
and torbernite (c.60°% U;O,, zero ThO,). Ranges of grain sizes used were: 
below 76 microns, 76-104 microns, 104-152 microns, extending, in the case of 
davidite, up to 850 microns. 

Count ratio measurements upon a variety of minerals are presented in Table 1. 
Values of the ratio C,/Cy were obtained with a statistical accuracy of +1-5%%. 
Taking the pitchblende sample as standard the results are given as a simple 
difference of count ratios. As the thorium—uranium ratio had not been determined 
for these ores by chemical or independent radiometric methods quantitative 
comparison is impossible. The negative count ratio difference for uraninite 
and torbernite indicates an apparent deficiency of radium. 


Table 1. Count Ratio Measurements upon Natural Mineral Samples 


oo Mineral a r (2), Remarks 

265 Pitchblende 0:0647 0) Uranium 

266 Uraninite 0:0625 —0-0022 Uranium 

267 Thorite 0-1130 0-0483 Thorium, possible 
traces of uranium 

252 ‘Torbernite 0-0582 —0-0065 Uranium, probable 
radium deficiency 

249 Uranothorite ~° 0-0679 0-0032 Thorium, some 
uranium 

238 Uranothorianite 0:0799 0:0152 Uranium and thorium 


in similar quantities 


Grain size of each sample 104 to 152 microns. 


As an estimate of the statistical accuracy of the energy discrimination method, 
the thorium—uranium ratio of a sample of activity equivalent to 4°% U,0, may 
be measured with error of less than + 5° in a counting time of 4 minutes. This 
excludes time devoted to counting any standard samples, and applies in the 
region of 7h/U=1. 


§4. THE BETA-ALPHA COINCIDENCE METHOD OF MEASURING Th/U 
4.1. Theory of Method 
An inspection of the various and widely spaced half-lives (1 psec. to 101° years) 
of the uranium and thorium series suggests a means of discriminating between 
the two series. For practical reasons the shorter half-lives are most amenable to 
treatment by electronic circuits: the disintegrations considered here will be those 
3 M-2 
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of radium C’ (140 sec.) and thorium C’ (0-3 sec.). By virtue of the Geiger— | 
Nuttall relation this discrimination by half-life may be considered complementary j 
to the method of discrimination by energy or range. | 

Each of the above disintegrations produces an alpha-particle and is preceded } 
by a beta-disintegration, from radium C and thorium C respectively. Assuming } 
the source to be presented simultaneously to alpha- and beta-particle counters, 1 
it is therefore necessary to design two gate circuits to be initiated by the beta- 4 
particle pulses and to accept alpha-particle pulses whilst in the open condition, 
The open periods or gate widths are fixed to enable mainly radium C’ alpha- 
particles to be counted through the ‘wide’ gate and mainly thorium C’ alpha- 
particles through the ‘narrow’ gate and in magnitude will be similar to the } 
respective half-lives. Thus each gate admits both types of related coincidences } 
together with unrelated or accidental coincidences, the latter being calculated | 
from the total alpha- and beta-particle counts. | 

The related coincidence counts may be calculated by considering the pro- - 
babilities of ThC’ and RaC’ disintegrations occurring within a time dr. These » 
are respectively Ap, exp (—Ayz) dr and Ayexp(—Ayt) dr where Aq, and Ay are : 
the respective disintegration constants. The number of ThC’ alpha-particle : 
counts in a counting period ¢, each occurring within 7 of a related beta-particle » 
count is 


mT hth, I Amexp(—Aqut) dr = mT hth (1 = expe eee (9) 


where m is the mass of the source (gm.) and ky, is an instrumental constant ; 
dependent upon the beta-particle activity of the ThC and the counting efficiency * 
for ThC beta- and ThC’ alpha-particles. It is assumed that the ThC present / 
is in equilibrium with a mass mTh of thorium. Similarly the number of RaC’ * 
alpha-particle disintegrations in a counting period ¢ each occurring within 7 of | 
a related beta-particle disintegration is 


mUtky{l—exp(—Agr)},  .  . .  veaeaes (10) 


disintegrations. 

Each coincidence gate admits related coincidences arising from both ThC’ ” 
and RaC’ disintegrations plus unrelated or accidental coincidences. If P is: 
the total alpha-particle count and Q the total beta-particle count the number of f. 
accidental coincidences admitted by the two gates during a counting period t/ 


are respectivel 
ote POrsjt: and. sPOsit cement ee (11) 


where 7, is the narrow gate width and 7, the wide gate width. By combination i 
of (9), (10) and (11) the total narrow coincidence count V and wide concent 
count W are given by 


v= mUtky{ 1 = exp (- AyT1)} aia mThtky 1 = exp (- Arn71)} == PQ7,/t, ] 12 II 
W = mUtky {1 — exp (—Ayry)} + mThtky,{1 — exp (—Aqyts)} + POr,/t. { a | 


{ 


The ratio of total coincidence counts, each corrected for accidentals, is given by | 


V' _V=PQr,/t _ ky{l—exp(—Ayrt,)}U+ Ref — exp (—Aqn7y)}TH 


W’ W- POQOrs/t A ky {1 Ta e49) (—AyT.)}U+ ky fl — exp (—Apnt2)}Th j 
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As in equation (6), this may be written in the form 
Vio 4,0 +95Th 


W' x,.U+y_Th’ 
where x=ky{l—exp(—Apr)} and y=ky,{1—exp(—A_yr)} with appropriate: 
suffixes. 

Continuing as in $3.1, 

LW ways YiY2 ‘Th 
ase en es G6 
(VW) —V"]W"  xtxyi- Hn Mar — Ho U ne 

The coincidence count ratio function is thus linearly related to Th/U. 

To minimize counting errors, optimum gate widths may be calculated to give 
the smallest standard deviations in V’ and W’. The optimum values of 7, and 
72; may be expressed in terms of the ratios of related to accidental coincidences 
V't/PQ7,, W't/PQr.. The magnitude of these ratios will depend upon experi- 
mental conditions, in particular upon the source—-counter geometry. 


4.2 Apparatus and Procedure 


The alpha-particle counter consists of a rectangular screen (area 32 cm?) 
coated with a layer of activated zinc sulphide and backed by three 931A photo- 
multiplier tubes connected in parallel. The ore samples, evenly spread on 
shallow copper trays are inserted immediately beneath the screen. The beta- 
particles are counted through the 0-005 in. base of the tray by two Twentieth- 
Century type B6 counter tubes situated underneath the tray. The combined 
counter is mounted in a pee. box. 
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Figure 5. Schematic diagram of beta—alpha coincidence counting equipment. 


As illustrated in Figure 5, the resultant alpha-particle pulses pass through a 
linear amplifier (type 1008) and pulse-amplitude discriminator (type 1028) and 
thence to the two gate circuits and the alpha-particle scaling unit. The beta- 
particle pulses pass via a quenching circuit (type 1014) to the beta-particle scaler 
and to the gate circuits. The narrow gate consists of a coincidence pair, width 
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being determined by a delay line, and is followed by an amplitude discriminator — 
and Post Office register. Since this type of circuit may be initiated by pulses | 
from both alpha- and beta-particles, the resolving time for accidental coinci- — 
dences is twice the gate width. Equation (12) should be interpreted accordingly. © 
The wide gate circuit consists of a square wave generator initiated by beta-particle _ 
pulses controlling a gating valve which admits alpha-particle pulses to a scaling _ 
unit during open periods. The particular experimental conditions gave optimum 
values of 0-43 and 115 usec. for the narrow and wide gate widths respectively. 
The narrow gate width was increased to 1 or 2 usec. to lessen the effect of variable 
rise time in the alpha-scintillation counter pulses, the amplitudes of which are | 
dependent upon the energy of the incident alpha-particles. The gate widths — 
could be approximately calibrated by an oscilloscope or more accurately by the 

conventional method using separate sources of alpha- and beta-particles. A | 
time-lag between the passage of an initiating beta-particle and the operation of a | 
beta-counter tube may modify the practical definition of the gate widths. This 


modification corresponds to lower limits, in the integrals leading to (9) and (10), — 


equal to the average time-lag. The effect of counter tube time-lags in coincidence 
counting has been discussed by Rotblat (1941): in the case of the B6 counter 
tubes the lag may be about 0-1 sec. and is accounted for by calibrating the 
apparatus with known standard ore samples. 

The equations (12) indicate that the errors in coincidence counting may be 
decreased by reducing the ratio of accidental to related coincidences. Ideally 
this ratio approaches a minimum when the source thickness becomes indefinitely 
small. In practice this degree of thinness is not attainable and the ore samples 
are spread as finely as possible on the trays. 

Except for very active samples the coincidence count rates are low: as long 
counting times are required for reasonable statistical accuracy, it is imperative 
to maintain stable conditions throughout the coincidence counting equipment. 
Most trouble was experienced in the alpha channel owing to a drifting count 
rate and was traced to variations in the type CV188 neon tube used for stabilization 
in the amplitude discriminator. The tube was replaced by a Mullard type 85A1, 
drift in alpha-particle counting rate being reduced to 1% per day. As an 
additional precaution the alpha-particle counting rate was continuously monitored 
by a stable counting ratemeter and pen recorder. 


4.3. Experimental Results and Discussion 


Equation (13) was tested with a range of ore samples containing pitchblende 
and thorium oxide in varying proportions. The experimental results obtained 
are illustrated by the graph of Figure 6, in which the coincidence count ratio 
function is plotted against the uranium-thorium ratio. The detailed coinci- 


dence counts are listed in Table 2. The range of samples was similar to but | 


richer (to provide higher counting rates) than those used in the energy discrimi- 


nation measurement of §3.3. In consequence there is an appreciable and | 


progressive bulk density variation from 1-8 to 2:2 gm/cm# over the range of samples. 
In addition to possible variable density effects in alpha-particle assay, similar 


to those encountered in the energy discrimination method, the variation of | 


beta-particle self-absorption may present similar difficulties on account of the 
different maximum energies of the related beta-particles from RaC (3-15 Mev.) 
and 'ThC (2:2Mev.). Fortunately, sample thicknesses are necessarily small 
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compared with the beta-particle ranges and this effect may introduce no serious 
errors in beta-particle counting. 

With the present experimental arrangement the statistical counting accuracy 
is such that the thori'um—uranium ratio may be measured to within +5°% for a 
sample of equivalent uranium concentration of 4°% after 50 minutes counting 
(these figures apply in the region of Th/U=1). An earlier, more sensitive, counter 
arrangement (consisting of six 931A photo-multiplier tubes in an outer ring 
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Figure 6. Coincidence count ratio function plotted against uranium-—thorium ratio. 


concentric with a cylindrical scintillation screen which surrounded the sample 
deposited upon tracing cloth fitted over a single inner cylindrical beta-particle 
counter tube) was rejected as too inconvenient and too liable to contamination 
for routine use. A more practicable system inverts the present counter 
arrangement; the beta-particle counter tubes are placed above the open tray 
which has a transparent Perspex base and contains sample and scintillation screen, 
both being considered expendable. The E.M.I. photo-multiplier with the larger 
cathode area registers the alpha-particle scintillations through the base of the tray. 


Table 2. Coincidence Count Ratio Functions 


Sample U ; ; V’ ; 
Tees. 4 is ‘i a ee 
= Pa _ os 3 eb 400900 

5 Oe 359 OS, L909 1814 0:1094 5-63 

6 4-0 429 10048 260 2083 O12 50 Nes 50 

8 LS 823 12604 616 2848 0:2163 Wert 

9 1-0 807 12491 602 2821 0-2134 1:73 

11 0-43 1246 15953 990 3873 O25 57 1:54 

14 0-25 1741 16511 1499 5081 0-2949 1-43 

15 0-11 2599 17427 2383 7237 0-32.92 1:38 


* By extrapolation. 
Counting time t=50 minutes; gate widths: 27,=—2:39+ 0-06 x 10~® seconds; 
T2=115:0+0:-4x10-* seconds; O=(V’/W)/[V’/W—(V’/W)y]. 


It is estimated that such an arrangement should measure 1°% samples to the 
same accuracy and within the same time as above. A further improvement 
in statistical accuracy is possible by delaying the operation of the wide gate by 
a few microseconds so that no related coincidences from ThC’ disintegrations are 


admitted. 
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Despite the initial satisfactory results reported above further development of © 


the beta-alpha coincidence method has been discontinued as it was felt the 
method and apparatus are too complicated for the purposes of routine geological | 
assay. Any future work would include consideration of a possible density 


effect, extension of the measurements to natural mineral samples and simplifi- _ 


cation and compression of the circuit and counting arrangements. 


§5. CONCLUSION 


Two methods of thorium—uranium ratio measurement in radioactive minerals | 
and ores have been developed and tested experimentally. Both methods provide | 
satisfactory counting accuracies using less than 1gm. of sample, comparing 
favourably in this respect with beta- and gamma-methods. Comparing the two, | 
the energy discrimination method is much simpler and more straightforward | 


than the beta—alpha coincidence method, which involves greater complication 


of circuits and apparatus to an extent that probably precludes use in routine assay. _ 


Systematic errors may arise due to variations in sample density, to diffusion 
and to equilibrium disturbance. The effect of density variations is likely to be 
most marked in the case of heterogeneous ores. Investigation of this effect is 
incomplete but for the moment errors in estimating the thorium—uranium ratio 
in this type of ore are reduced by calibrating with standard samples of similar 


density. Diffusion errors are minimized by similar treatment of unknown and © 


standard samples during grinding and preparation. ‘The effect of other dis- 


turbances of equilibrium on the energy discrimination method has been indicated ~ 


and in certain cases it is possible to measure the amount of the disturbance. 


Estimated statistical errors in the determinations of Th/U are: energy discri- — 
mination method: +5% in a 4% (U,O, equivalent concentration) sample in — 
4 minutes counting; beta—alpha coincidence method: +5% in a 4% (U,O, © 
equivalent concentration) sample in 50 minutes counting. These figures apply — 


in the region of Th/U=1. 
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ABSTRACT. An expression is derived for the flow field in a two-dimensional jet of ideal 
fluid impinging normally upon a plate situated at an arbitrary distance from a parallel sided 
orifice. Numerical data are given when the orifice is at an infinite distance and in three other 
cases. In each example the trajectories of solid particles carried by the fluid are calculated 
by a stepwise method, and the conditions determining whether a particle will strike the plate 
or not are found. 


§1. INTRODUCTION 

HE problem discussed in this paper arose in connection with the study of 

| dust-sampling instruments in which particles are deposited from dusty 

air, on account of their mass, when an air jet is directed on to a plate. 

An experimental investigation has been made using different kinds of particles, 

from 4 to 10 microns in diameter, and jets ranging from low-velocity slits, several 

millimetres wide, down to small orifices in which the air approached the speed 
of sound. 

The experiments enabled practical information to be collected about the 
design of the jet, the adhesion of dust to the precipitation plate and the break-up 
of particles on impact, as well as other factors. ‘To assist in interpreting the 
experimental results it was considered worth while to explore, theoretically, the 
performance of an idealized impingement system; this would make clear and 
definite the purely dynamical factors involved and would set up a criterion of 
efficiency for the collection of particles against which real apparatus could be 
compared. 

The equations of motion of a massive particle moving in a fluid have been 
solved by G. I. Taylor (1940) for conditions near a stagnation point where the 
streamlines of fluid flow approximate to rectangular hyperbolae. He followed 
the idea of Albrecht (1931) in considering the fluid ideal, that is, non-viscous, as 
far as the flow field was concerned, but taking it as viscous relative to the particles 
so that the drag upon them could be calculated by Stokes’ law. The particles 
were assumed to be massive points, creating negligible disturbance, whose size 
only came into the expression for their drag. 

In the dust-sampling instrument, however, the air jet has finite width, while 
the distance of the jet orifice from the plate is variable and exerts an important 
influence upon performance. 

The two-dimensional flow field for an ideal fluid emerging from a parallel- 
sided opening at an arbitrary distance from a plate can be calculated by a 
conformal representation method. Since the jets in dust-sampling instruments 
are frequently slit-shaped, it seemed that this method would give as close an 
approach to real conditions as could be expected, while preserving a rigid 
dynamic definition. The complexity of the expressions derived in this way 
for the flow field made it clear that no direct solution of the particle equations 
would be possible on the lines of Taylor’s attack on the simpler problem. 
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Accordingly, the particle trajectories were calculated step by step after the flow , 
fields had been fully computed from theoretical formulae. ' | 

Discrepancies between the computed and the real flows must arise on account | 
of the viscosity and compressibility of air, the contour of the jet intake and its : 
spread due to turbulence. In the theoretical work the air stream 1s supposed 
to be laminar and to flow with uniform velocity U, between parallel planes which | 
terminate at the slit orifice at distance d from the plate (Figure 1); the stream 
then spreads out symmetrically against the plate and rises to velocity Vo at right ¢ 
angles to the original flow. 

In practice, the velocity at the walls of the jet and on the surface of the > 
impingement plate is zero. A velocity gradient therefore exists adjacent to the > 
solid boundaries, and the velocities at infinity are not constant and equal to) 
U, and V, across the incident and deflected streams as assumed. | 

The Reynolds’ numbers of the jets in dust-sampling instruments are usually / 
high, and in most cases the stream approaches the orifice between converging | 
walls. These circumstances favour the existence of a fairly flat velocity contour ~ 
across the jet section together with a thin boundary layer. It is considered that t 
ideal fluid theory should provide a reasonable approximation in this respect. . 
The existence of a boundary layer along the plate is probably of little significance. . 

Converging airflow also tends to inhibit turbulence in the stream approaching ¢ 
the orifice. When the jet emerges, the boundary layer forms a region of { 
discontinuity between the rapidly flowing jet and the stagnant air outside. , 
Ultimately, therefore, it will roll up into vortices and give rise to an increasing 7 
zone of turbulence in the outer parts of the jet. This behaviour was observed |) 
in one type of instrument, but seemed to have little effect on the dynamics of * 
particle impingement. 

Other possible causes of difference between theory and experiment, which 
are related to fluid mechanics, are the compressibility of air and the existence of © 
a contracted vein outside the orifice. In general, however, it was found that the 
most striking deviations of the experimental results were due to particles 
shattering on impact with the plate or else failing to adhere. : 

The assumption of Stokes’ law of resistance in calculating the particle 
trajectories is a close approximation unless the relative motion of particle and | 
stream becomes too great. A correction can be introduced by using the known | 
relationship between Reynolds’ number and drag coefficient for spheres; since » 
the velocity difference is not constant the drag will, in fact, not be given accurately, | 
but for spheres representing dust particles error on this account should be small. . 
Diffusion of particles is quite negligible when they exceed diameters of about : 
0-2 micron in systems like those considered, and no account need be taken of it. H 

In the theoretical work the particle was assumed to start with the stream | 
some distance upstream of the orifice. Actually, heavy particles lag behind the | 
stream owing to the acceleration of the air approaching the jet, but discussion | 
of such points is beyond the scope of the present paper. | 


§2. THE FLOW FIELD IN A TWO-DIMENSIONAL JET OF IDEAL 

FLUID IMPINGING NORMALLY UPON A PLATE | 

Let z=x'+iy’. Then the boundaries of the parallel-sided jet AB, A’B’ | 

and the plate DED’ are shown on the g plane in Figure 1 with the origin at E 

and the x’ axis along EF, which is also an axis of symmetry. The width of the jet 
is AA’ = BB’ =h, and the distance of the orifice from the plate is d. 
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The jet emerges and spreads sideways over the plate, being bounded by the 
free streamlines BC, B’C’ along which the.pressure, and therefore the velocity, 
is constant. AFA’, CD and C’D’ are at infinity. | 

At AA’ the stream velocity is — U), at CD it is V, and at C’D’ it is —V,. 
Let V,/U,=«. From continuity, CD =C’D’ =h/2«. 

If the pressure at AA’ exceeds that outside the free streamline by an amount P, 
then equation of the potential energy at AA’ to the total energy at CD, C’D’ gives 
PhU,) =3p(V o°h/x— U,*h), where p is the density of the fluid, so that 
P=tpU,7(«@7—1). Since P20, «21. When d is infinite, «.=1 and BB’ 
coincides with AA’; otherwise, for a finite value of d, «>1. 

If U’ and V’ are the resolutes of the velocity q parallel to x’ and y’, then 
U’=q cos 6, V’=qsin@. If ¢ is the velocity potential and s the stream function, 
then U’=—d¢/0x’, V’=—dd¢d/ey’, and U'=—ddb/dy’, V'=db/ex'. The 
assignment of ideal properties to the fluid, implying volume continuity and 
zero vorticity, leads to V76=V?4=0. Thus ¢ and # are conjugate functions 
obeying the Cauchy—Riemann conditions, so that we can put w=6 +7 =f(x’ +1y’). 

We now transform the boundaries in the z plane to the w plane. Let 6=0 
at E, then ¢ tends to oo at AA’ and to — o at CD, C’D’. Let Y=0 along FED, 
FED’, then 4 = U)h/2 along ABC and — U,h/2 along A’B’C’. The w plane is 
shown in Figure 2. 


hy 4 
E Fy! 
_& 
a 
pyc! 
Figure 1. 2 plane. Figure 2. w plane. 


Since the boundaries of flow are now straight-lines they constitute a polygon 
with vertices ABCDED’C’B’A’. The interior of this polygon can be mapped 
on the positive half of a ¢ plane, so that the vertices lie at corresponding points 
along the real axis of ¢, by applying the Schwarz—Christoffel theorem. 
Proceeding, let E fall at €=0, C’D’ at €=1, CD at (= —1, etc., as shown in 
Figure 3. Then dw/d€=K(€+1)-%¢(¢—1)1, which, from a consideration of 
Figures 2 and 3, integrates to 

pee exp (dy), | gee (1) 
where a=27/U,h. 


A -In(Ya)- 
Figure 3. ¢ plane. Figure 4. @ plane. 
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By means of the Kirchhoff transformation a relation can be obtained between — 
dw/dz and ¢. 

Let Q=In(U deidw). tv 1 PS (2) 
Now it can be shown that —dw/dz=U’—iV’=Q=qe~”, where Q=U'+iV’", | 
so that Q=In (U,/q)+1(8—7). This relationship enables the boundaries to be | 
plotted on the Q plane (Figure 4), where they again appear as straight lines ! 
since g is constant along BC, B’C’. A further application of the Schwarz— | 
Christoffel theorem then maps the interior of the polygon in the Q plane on to | 
the ¢ plane with the vertices as before. Hence 


dQ|dl = K(C+)9(C4 1) 22E-U(6 — 1126-4”, 


from which, by integration, Q=In {(f —b?)1? +0(@ —1)¥?}—In {a£(b? — 1)1?}, 

where b=(2+1)/(a2—1), for the field A’B’C’D’E, where ¢€>0; since | 

1<b<a, «>1. ) 
Eliminating Q between this expression and (2) we find 


Udz/dw = (a? — 1) {(2 — b?)1? + b(@@ —1)4?}/2022. (3) 
Also, from (1), dw/d€=2¢/a(¢?—1). Eliminating dw, 
Upz = {(a* — I) Jaa }{ f bde/ (2? — 1)¥? + f (2 — bye) (? — 1)} 
U actz (6+(C = = ae ee i — Beye 4 4 
1 


or a2 + 1 ach © Ded) (2@—B)R—L 
an C4 i(€2 — b2)12 — Qak /(a? —1) 
oF aL CEL eea i Cesiien ne: (t 


At B’, x’ =d, y’ =h/2 and ¢=8, so that, from equation (4), d can be expressed | 
as a function of «, 


h {(o2+1)In(a+1)/(a—1 
i= 5 {eae er) # i} . (5) 


Since — dw/dz =U’ —1V’ = O, equation (3) gives 
UplO=(1—a2){(1 82/2)" +01 -1/P)}2e2 cane (6) 


Elimination of ¢ between equations (4) and (6) will therefore lead to a relationship _ 
between zand Q. This elimination is carried out by substituting t=(1—22/2)¥2_ 
into equations (4) and (6) so that the latter reduces to t= (a2 Oe Q/ U,)/(1—22). | 
ae this value of t into equation (4) and separating real and imaginary parts | 
we have 


2x! +1, (a2 + 1)4 — 2(a? + 1)?(402 + A? cos 27) + of 


A 


= = N- OOOOD— 
ho Ano [a2 +1)? — (0? + 1) f2(o2 + 4a2 + A® cos 27) PP + oP 
w—1, A*—2(o2—1)2A? cos2r+(a2—1)! 1 40.4 
ate > hh Ne elds asta 
4na® " TA®—2(2 —1)A cost + (a1 | ma 2" ga age’ 
fy 2y' eeoote fee (a? + 1) {2(0? — 4a? — A? cos 27) }1/2 a 
h Tr. (a? + 1)? —o? 


a1 2g? = Asta aed 4 
+ tan) = heh 
Daroe (2-1 A® * Ine Sigpa ee, | 
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where A cos t= —cos 0(a2U,/q—9/ Up), Asin t= —sin 6(a2U,/q+q/Up), 
o* =16a*+ 8a2A? cos 27+ A4. If «=1, these expressions reduce to 


Sa =(1 U,/¢+94/U)+2 cos 0 _, —2q 08 a } 


x= 4 ln —*+_+ _* ——_ + tan 


h a\" Uo/¢g+q/U)—2 cos 6 1—q@/U,? * 
_ 2y" _2 tan-? — 24 8in 8/ Uo tn Uolg t+ a/Uy +2 sin G | hag) iets 
ey h vin 1—q¢/U; 2 U,/g+q/U,—2siné > ] 


which represent the flow field when the jet is at an infinite distance from the plate. 
In this case symmetry exists about lines through the origin making angles 7/4 
with the axes. 

Along the free streamline B’C’, q is equal to V) and r=7/2. Equations (7) 
then reduce to the equations of the free streamline as follows: 


ee bl noe tl 2x cos 8 1 | 
; 2ra2 = o2+14+2xcos@ «’ 
F F ee Ae Se (9) 
gett Basle i Ho" + ey niasing [ 
2702 mo 2asind ma 1+sin@” | 


Along A’B’ the velocity rises from — U, to — Vy and we have 


ee i eee tent 
Ore Ge i (4a 2 A*) Os den? eal A ae A fe 


while from F to E the velocity falls from — U5 to zero and 


ee! 4 a? + 1+ (402+ A?) gece _ a?—1+4+A eile a, 4aA 
— Dara® ~— 02 —-14+ (402+ A2)"2 © Qero® ~ —0? +14A ° orm A? — 402 
fy = 0. 
Na ice (11) 
Finally, from E to D’ the velocity rises from 0 to — Vp) and 
x=0 
Beene gl yee 1 ate pe (12) 
ei (A? — 402)? © ara? cee at 2 A 2a 


The same order of the terms has been preserved in all these expressions. 


§3. COMPUTATION OF FLOW FIELDS 
The formulae derived above relate to the field A’B’C’D’EF. Alternative 
quadrants can be represented by a suitable selection of signs at an early stage 
ofthe argument. ‘Table 1 gives values of the velocity and its direction at important 
points on the flow boundaries. 


Table 1. Values of U’, V’, g and 6 at Salient Points on the Boundaries of the 


Flow Field 
Point A’ Be Cs |DY E F 
(Gi —U, —V, 0 0 0 —U,) 
Va 0) 0 — V5 —V, 0 0 
q Uy Vo Vo Vo 0 Uy 
6 7 7 3277/2 32/2 32/2 7 


For «=1 equations (8) were computed directly for a range of values of 
tan @ and q/U), and the results are exhibited in Table 2(a)-(d). By reflecting 
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this table about the bisector of the angle between the coordinate axes the 
remainder of the flow field can be obtained. Table 3 contains the coordinate 


of points along the free streamline. 


Table 3. Free Streamline. «=1 
| 1-001 1-013 1-019 1:026 1-038 1-051 1-064 1-076 
| 5:398 3-932 3:674 3-491 3-233 3-050 2-908 2-792 | 
| 1-102 127 1-191 1-255 1-382 1-509 1-561 
| 2-609 2-468 Paes 2°033 1-785 1-616 1:561 


Since In U,/q and @—7 are conjugate functions of x and y an le 
network is plotted on joining up the points given in Table 2 to make lines 0 
constant g and constant 6. This facilitates the interpolation which is required 
later when drawing particle trajectories. 
Using equation (5), the distance d of the orifice from the plate was calculated) 
as a function of «, the ratio of the velocities of the deflected and incident streams4, 
This is plotted in Figure 5. i 
Equations (10), (11) and (12) were then used to calculate the velocities round) 
the boundaries of the flow fields taking « as 1-0, 1-6, 2-5, 5-0 and 10-0. These! 
curves are plotted on Figures 6, 7 and 8. Figure 9 shows the free streamlines; 
for the same values of «. t 


02 04 O68 O8 0 TF 
2a/h 
Figure 5. Ratio of velocities of deflected and 
incident streams as a function of the 
distance between orifice and plate. 
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Figure 6. Velocity along wall of jet. Tigure 8. Velocity along the x axis. 
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_ Guided by the values of U and V round the boundaries, the flow fields for « 
to 1-6 and 2-5 were calculated by the method of relaxation. This part 
the work was performed by Scientific Computing Services Ltd. Their 
results are shown in Table 4(a) and (4), U and V, and Table 5, free streamline. 
Corresponding figures for «=2-5 appear in Tables 6 (a), (b) and 7. This system 
of computation gives a rectangular grid of points. 

For «=5 the theoretical formulae were again employed ; the flow field is given 
in Table 8 and the free streamline in Table 9. Orthogonal curves are once more 
obtained for g and 6 and have been drawn in Figure 10. 

This completed the establishment of the flow fields and made it possible to 


proceed to study the trajectories of particles carried along by the fluid leaving 
the jet orifice. 


0 02 04 06 08 10 12 14 
x 


= 
6 04 os iia e ¥ 26 24 28 


Figure 9. Free streamlines, Figure 10. py, curves 
of g and @ for «= 


§4. THE EQUATIONS OF MOTION OF THE PARTICLES 
Let u’, v’ be the resolutes of particle velocity and U’, V’ as hitherto the 
resolutes of stream velocity. Assuming Stokes’ law of resistance, the equations 
of motion of a particle are 


& 


m du’ oes m dv ete 
eae and dee 

here r is the radius of the particle, m its mass and 7 the viscosity of the fluid. 

The equations can be rendered dimensionless by measuring the velocities in 

units of U,, the distances in units of $h, the half width of the jet and the times 

im units of h/2U,. Hence u=u'/U,, x=x/4h, t=t'U,/th, etc. The equations 


transform to eS : oe 


where P=2mU,/6zrnh is the dimensionless particle parameter. It is actually 
equal to the distance (measured in units of $h) that a particle will travel through 
fluid at rest when its initial velocity is U5. 


§5. COMPUTATION OF PARTICLE TRAJECTORIES 
The step-by-step method of computing particle trajectories can be followed 
by referring to Table 10 in which the complete derivation of a trajectory for 
P=0-046 and «=5 is set out. The particle is started at a given value of y with 
_ PROC. PHYS. SOC: LXIV, 1oO—B 3N 
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x—>0o. This is interrupted as the value for « for which the stream velocity is about | 
0-99. The particle of Table 10, was launched at x, = 1-4, y, =0-1; the suffix n | 
refers to the number of the step and is given in column 1 of the Table. Opposite | 
the row n are entered the particle velocities u,, v,, in columns 2 and 3, and | 
coordinates «,, y,, in columns 10 and 11. The length of the step resolute dx is | 
then decided and entered in column 8. Estimates are now made of the velocity — 


changes 5,u, 5,v during the step (columns 4 and 5); the suffix denotes that these | 
are approximations. 
Table 5. Free Streamline for « =1-6 

x y 1,000U =1,000 V x y 1,000 U 1,000 V 
1:2741 1:0 1600 0) 0:6476 2:2 90°8 1597-4 
25 ~ 1:0026 1580 249 0:6426 2s 70:7 = =1598-4 
te? 1:0141 1533 457 00-6387 2:4 55:0 1599-1 
1°15 1-0320 1476 617 0-6356 2°5 42:7 1599-4 
1-:1113 1-05 1424 730 0-6333 2°6 Boi 1599-7 
1:1 1:0559 1407 761 0-6314 2:7 25:9 1599-8 
1:05 1:0865 1325 897 0-6300 2°8 20-1 1599-9 
1:0306 13 | 1289 948 0-6289 2-9 15:6 1599-9 
1:0 1:1240 1227 1026 0-6280 3-0 12-2 1600-0. 
0-9708 1-5 1163 1099 0:6274 3°1 9-5 1600-0 
0-9230 12 1045 1211 0-6268 3-2 7-4 1600-0 
0:9 1-2281 983 1263 0:6264 3°3 5:7 1600-0 
0:8504 1:3 835 1365 0:6261 3-4 4-4 1600-0 
0:8 1:3948 . 669 1453 0:6259 3°5 3:5 1600-0 
0:7976 1-4 661 1457 0:-6257 3-6 2:7. 1600-0 
0:7581 1-5 520 1513 . 0-6255 3°7 2-1 1600-0 
0-7280 1-6 406°8 1547-4 00-6254 3°8 1-6 1600-0 
0-7048 1P37) 317°6 ©1568-2 00-6253 3-9 1-3.) 160020 
0:7 1-7247 299 L572 0-6252 4-0 1:0 1600-0 
0-6870 1:8 247:6 1580-7 0-6252 4-1 0-8 1600-0 
0-6732 1-9 192:9 1588-3 0:-6251 4-2 0-6 1600-0 
0-6624 2:0 150:2 1592-9 0:-6251 4:3 0°5 1600-0 
0-6541 2-1 116°8 1595-7 0°6251 4-4 0-4 1600-0 


Suppose ot is the time taken by the particle to cover the first step. Then 
dt ~dx/(u,, + 35,u) —Sy/(v, +45,v). The value of Sy, corresponding to the 
resolute 5x, is calculated from 


oy = dx(v, + 30,v)/(u,, Ze 30,U), AOS: ore eS .9 (14) 


using columns 6 and 7, and is entered in column 9. With Sx and dy known, the | 
diagram of the flow field can be consulted to interpolate U,, U,,4,, and V,, Vay 
the fluid velocity resolutes at the beginning and end of the step (columns 12 and 
18). Then, from equation (13), we can calculate a more accurate value of Sw as 


follows: 


ot dx (4(U,+U. | 
b= BLU n+ Onin) (ln + Bad} = B {eos ah os: (15) || 


u, +46,u 
This is entered in column 17 while columns 13-16 facilitate the arithmetic. | 


S larly 5) oa ] sae 
1m1 ) oy 1(] 2 J 1) 16 
3) aN nee, eee ( ) 


gol 


Impingement of Particles 


9e01 
FEOL 
8cOl 
6L0L 
6001 


ISOL 
LyOL 
8col 
9COL 
CLOL 
L66 
£86 
1L6 
£96 
686 
LS6 


$1 


o) 


0 


COLT 
£601 
+01 
LyOl 
Stor 
166 
996 
9b6 
(4%0) 
£76 
076 


(An 


(nm annntt mT eantpa)\ TIACQMIYVPvVPIAnr ka NOWUPINAIYA C7 —N WMT HNIAT.T AOTT 


0 
cS 
16 
SIL 
€cl 
oll 
SOL 
+8 
6S 
Of 
0 


EL 


Jo sanjea) uonexejoy Aq poye[noyeo ¢-7=% IOF Play MOLY 


SPIT 
Cell 
coll 
e901 
1ZOt 
£86 
1S6 
976 
806 
L68 
+68 


i+ 


0 
c8 
LE 
p91 
OLT 
6ST 
Sel 
601 
SZ 
8e 
0 


0-1 


60CL 
O6TL 
6ELL 
8LOL 
OCOL 
OL6 
6C6 
668 
LL8 
$98 
098 


O-1 


0 
Sel 
LIZ 
LET 
€fC 
Lie 
SLT 
Sel 
+6 
8b 
0 


6:0 


OcEL 
OLE 
98Tl 
coOl 
cLOL 
8h6 
868 
798 
9€8 
18 
918 


6:0 


0 
8SC 
9EE 
Ive 
ble 
$Le 
ScC 
ILI 
Stl 
8S 
0 


8-0 


6EST 
80+r1 
CETL 
£601 
686 
C16 
$88 
£18 
S8L 
69L 
£9L 


8-0 


L:0 


Lert 
9067 
OST 
Sect 
c90L 
I+6 
9S8 
FOL 
Lee 
€cL 
SOL 
669 


4:0 


SOET 
LC6 
Orl 
€l9 
60S 
9Ih 
8CE 
Eve 
191 
08 
0 


9-0 


FOOL 
O0TI 
8891 
SLL 
cell 
1L6 
6S8 
LLL 
SIZ 
SL9 
9+9 
0e9 
+09 


9-0 


09C7 
9881 
L9+1 
6cTL 
106 
EL 
86S 
C84 
OLE 
LLC 
cst 
16 
0 


S-0 


L98 
OLOL 
elect 
S60T 
956 
Les 
Obl 
929 
£c9 
S8s 
6SS 
brs 
oes 


gO 


0OSZ 
664C 
86¢7 
96¢7 
88h7 
bLyT 
EVV 
18&7¢ 
$9CT 
OL0C 
86L1 
66+1 
LCL 
FOOT 
€¢8 
699 
LES 
Liv 
90€ 
10c 
OOT 
0 


+-0 


86 

crt 
807 
£0€ 
Ofr 
$8S 
Chl 
sss 
6L8 
O€8 
FSL 
8L9 
cl9 
9SS 
+IS 
C8P 
Lor 
Sth 
ber 


+0 


6877 T8r7 
E84TC OLY 
VLYT PST 
O97 CEPT 
SEPT 667 
COT TSET 
SPET £877 
¥9T7C 68172 
LAG 1907 
L961 8681 
IvLt SOLT 
COSt S641 
CLGL 0671 
+901 S601 
£88 L16 
EGE SSL 
6LS 809 
OS+ CLy 
O£E LYE 
9I~ Lt 
LOL CLI 
0 0) 
£:0 c:0 
(9) 9 21921, 
Ss c9 
£Tl 68 
OLT Get 
OFT TLT 
CLE S77 
(Ac +87 
97S SEE 
£6S LLE 
919 £6€e 
86S 98¢ 
LSS 99¢ 
OLS 6£e 
Sor GES 
Str ESC 
FOL 997 
OLE Sz 
HSE OFT 
Ere €£T 
O+rE O£T~ 
£-0 c:0 


efnm\ A ATOAD T 


Ce 

9+ 

$9 

88 

ell 
Ol 
99T 
+81 
lol 
Oot 
c8T 
OLT 
9ST 
trl 
bel 
9CT 
1cl 
Sil 
Lit 


1-0 


oO cocococooc“neocococooeco 


A 


Deo oarceoeqQosse 4S 2aq 


AAMPMODOHDOAAMTNORADOG 


o 


Z 


ool lo lokolo loo ee 


AAMNFNOCRDDHDOGANMTMNOND 


go2 C. N. Davies and Mary Aylward 


is evaluated from equation (14) (column 9) with the aid of columns 19-22. év is | 
entered in column 23. | 

If the values of 5u and Sv in columns 17 and 23 are not in fair agreement | 
with the estimates 5,u, 5,v the step is repeated using intermediate values. Steps10 | 
and 13 in Table 10 had to be repeated for this reason. 

As the particle approaches the plate the size chosen for dx is decreased because | 
the path turns round and tends to run parallel with the y axis. In the example 
illustrated the particle actually struck the plate with a residual velocity com- 
ponent u equal to 0-029. This was for P=0-046; when the process was repeated 


Table 7. Free Streamline. «=2:5 


53 y GS Va x y (Uy V 
0-713 1:0 DS 0) 0:7 1-002 2:432 0-578 
0-614 1:05 1-893 1-633 0:65 1:024 2-137 1-297 
0:567 1:1 1°535 1:973 0-6 1-063 1-792 1-744 
0-533 1:15 1-255 2:162 0:55 1-123 1:397 2-074 
0-508 1:2 1:029 2-279 0:5 1-218 0-958 2-309 
0-488 1:25 0-845 2°353 0:45 1-389 0-488 2-452 
0:471 1:3 0-693 2:402 0-4 ra) 0 2-500 
0-448 1:4 0-468 2-456 
0-432 1:5 0-316 2-480 
0-422 1:6 0-213 2-491 
0-415 1:7 0-144 2-496 
0-410 1:8 
0-407 1:9 
0-405 2:0 


for P=0-025, u,, became equal to zero before the plate was encountered and a 
negative value was obtained on subtracting du from the last value of w,. By 
constructing a graph of P against uw, (positive or negative) at x =0 it was possible 
to interpolate and find accurately the value of P, for which the particle just met 
the plate with zero velocity component uw. This is the critical value P, and the 
corresponding trajectory is the critical trajectory. P, depends on the starting 
point of the particle y,, and a number of trajectories were calculated for different 
values of . 

The interpolation of P, was linear except in one case for which a quadratic 
was required. When P= P, impact with the plate takes place at y= oo. 

If P, was assessed accurately in the first instance, it became difficult to tell 
if the particle would just hit or just miss as it was swept along close to and nearly 
parallel with the plate. The uncertainty was resolved as follows: Suppose at 
the end of the 7th step at (y;, «,) the particle velocity was (u,, v,), whereas the stream 
velocity was practically (0, Vy). ‘The equations of motion of the particle then 


d 
gi oe P du/dt=u or u—Uu,;=(x;,—x)/P 
} Vt Pee ee SY <a coi (La 
and — yet f= Stn (1+ “ i 
Ls es ye aap Pu, 


The first expression gives u>0 when x=0 for a hit and u<0 for a miss, 
remembering that the particle velocity is always negative. The criterion for a 
critical value of P is u=0. The second equation shows that y=co at x=0 
forP=P- 
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§6. DISCUSSION OF PARTICLE TRAJECTORIES 
The critical trajectories of particles starting at a number of different positions 
in the incident air stream (0 <y, <1) are plotted in Figures 11 to 14. All the curves 
meet the plate x=0 at y-+oo. When P>P, a curve lying below the critical path 
is obtained which meets the plate at a finite value of y. For P<P, the particles 


fail to strike the plate and travel parallel to it with the emergent fluid stream 
at a finite value of x. 


0-1 
== a as 10-02 a — _——— 
0 OAT 0-8 As L:2 F620) 2024-228 32° 3:6. 6-4-0 0. 0% 050 O75 100 125 150 75 
7 x 
Figure 11. Critical trajectories. o«=—1. Figure 12. Critical trajectories. 
a=1°6. 
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Figure 13, Critical trajectories. x=2:°5. 
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Table 9. Free Streamline. «=5 


tan@d «x y tanO «x y tan 0 x y 

0-1 0:334 1:0 0-6 0:313 1-009 2:5 0:248 1:097 
0:15 0°333 1:0 0-8 0:303 1-015 3-0 0-240 1-116 
O=2 023320 1:0 0-292 1:024 4-0 0:231 1-149 
0:3 0:328 1-002 1:05 0 -O Sse Osan 5:0: 02225, st i777, 
0-4 0:324 1:004 1-5 0-272 ‘1-050 70> .Oe20Se es deoaGh 
0:5 0-319 1-005 20s 0: 253ee 1075 10:0 0:213 1-267 
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0 0:2 04 06 08 1-0 12 
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Figure 14. Critical trajectories. «=5. 


In Figure 15, for the same set of values of «, P, is shown as a function of y,, | 
the starting coordinate. The smallest value of P, resulting in impingement is © 
achieved in the case of particles launched near the axis. Conditions on the axis 
require further consideration and will be studied later. 
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Figure 15. Variation of the critical particle size parameter Pe with the starting coordinate 4,. 
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For the two greatest values of « there is a value of giving minimum efficiency j 
of impingement since P, rises to maximal values at ), =0-65 (a=2-5) andi 
4, = 0-45 (a =5). 

If particles of the same size and density are uniformly distributed in the = 
approaching air stream the proportion striking the plate depends on the = 
relationship between P, and y,. For «=1-0 and 1-6, where P, increases steadily ; 
with y,, the graph of the efficiency of impingement, or the fraction of the incident | 
cloud of particles which hits the plate, is the same as the (y,, P,) curve; 100% , 
efficiency corresponds to y,;=1, and so on. With «=2-5 or 5, owing to the: 
maximum of P,, particles from a certain region of y, may miss the plate, while : 
those starting at smaller and greater values strike it. The efficiency curve in this s 
event is constructed by making due allowance for this above a certain value of P,. \ 

Graphs of efficiency of impingement as a function of P are plotted in Figure 16) 
for different values of «. These curves can be interpreted as expressing the» 
performance of an ideal dust-sampling instrument on the assumption that all | 
particles striking the plate actually adhere to it, whatever the angle of impact. J 
Collisions between airborne particle have been neglected. 


| 


Efficiency (%) 


0 01 0-2 0-3 0-4 05 0-6 07 0:8 09 10 
P=2mlh/6urnh 


Figure 16. Efficiency of impingement of particles. 


It is clear that the impingement of small particles is facilitated by diminishing d, i 
the distance of the orifice from the plate; this also improves the selectivity of t 
the system since the range of values of P over which impingement changes from 1 
zero to 100% narrows considerably. Figure 17 illustrates the last point. When i 
a is less than 2:5 the range of P spreads out, while a sharp cut-off is obtained | 
for«=3,. ‘There would seem to be no advantage in increasing « above 3 since the: 
accompanying reduction in particle size can also be obtained by increasing the: 
velocity Uy or decreasing the width h of slit. | 


| 


Calculation of P for different particle radii (microns) and different values of { 
the ratio Up/h (sec +) is rapid with the help of Figure 18. . 


| 
i 


§7, PARTICLES MOVING ALONG THE X AXIS 


The velocities of particles moving along the x axis were computed by the 
step method and also by approximating the flow field to a simple form allowing 
integration of the equation of motion of the particle. Inspection of Figure 8 
shows that the stream velocity along the x axis tends linearly to zero with 2.., 
Only some distance away from the plate is the deceleration more gradual; it 


a 
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was thus thought that a reasonable approximation to the motion of the particles 
could be made by assuming U=cx; U=U), x=x, at t=0. The solution of the 
equation of motion then becomes 


N= {cos ptt iF (1-2P,) sinpt} Cer ee (18) 
for P>1/4c, and 


Shi See aby MLE ie —atj2 | ,—tl2P 
w= 3'{(p+4 2c) e! € et OA WS ac ie A Nt las as ere (19) 


for P<1/4c, where p =(4cP—1)"?/2P and gq =(1—4cP)"?/P. 


Efficiency of Impingement (%) 2 


0 0-1 0-2 03 0-4 05 06 07 0:8 TT, 
Iz 


Figure 17. Illustrating how the sharpness of cut-off depends on «. 


Particle Radius 7 (microns) 


P=2mlh/6urqh 


Figure 18. Diagram for calculating P in air at 20°c. and 760mm.Hg. The parallel, diagonal 
lines give the relationship between 7 and P for different values of Uo/h. 


In the first case the particle strikes the plate with a finite velocity, and in the 
second with zero velocity. However, it is not travelling with the stream unless 
Ptendsto0. The critical case is P,=1/4c when x = («y+ Upt) exp (—t/2P,). 

This method of calculation was quite reliable near the plate but made P. rather 
larger than the stepwise method, which was used finally for all values of «. 

In Table 11 values of P, calculated by both methods are given for particles on 
the x axis. 
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Table 11. Critical Values of the Particle Parameter P, on and near the Axis 


a Stepwise method Pe— VAC 
Vi—0 4y,=0-02 
1 0-60 0-345 0-64 | 
1-6 0-31 0-25 0-325 a 
2:5 0-195 0-16 0-225 | 
5-0 0-140 0-025 0-1875 i 
oO a a 0-1625 


Figure 19 illustrates the velocity of impact as a function of P for « =1, 1-6, 2-55 
and 5, while Figure 20, with « =5, shows the velocity as a function of x. 

P,, has not the same significance for particles on the axis as it had for those start- - 
ing with non-zero values of y, since, clearly, all particles on the axis must strike» 
the plate. Examination of Figures 11 to 14 will show that the curve of y, plotted | 


2 | a i ee 
18 Stepwise Calculation 
ah ----- Assuming U= cx 


i es 0-4 0-5 06 
(Velocity of Impact)/% 


Figure 19. Impact velocities of particles on the x axis. 
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Figure 20. Velocities of particles travelling along the x axis for «=5. 


Impingement of Particles QI 


against P, seems to meet the axis y, =0 at a value of P, which is substantially less 
than that found for a particle on the axis. The smallest finite value of y, for which 
trajectories were computed was 0-02, and what happens to the curve between here 
and zero is conjectural. It seems likely that there is a discontinuity, so that, in 
effect, a definite value of P exists below which no impingement takes place. 
Probably the cut-off occurs close to the values of P, for y, =0-02, and these are 
shown in Table 11 for comparison with the figures for the axis. 
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Fluctuations in the Energy-Loss of Fast Electrons 
in a Proportional Counter 
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ABSTRACT. The fluctuations in the energy lost by fast electrons over a fixed small 
fraction of their range have been studied by measuring the pulses the electrons produce in 
passing through a proportional counter. The energy loss distribution calculated theoreti- 
cally by Landau is approximately gaussian with a high-energy tail. The experimental 
distribution is in agreement with the theoretical on the high-energy side of the peak, but 
is wider than Landau predicts on the low-energy side; this discrepancy is greater in krypton 
than in argon. The experimental! values of the most probable energy loss agree with the 
theoretical values. 


$i. INTRODUCTION 

HEN a fast particle passes through a layer of matter, it loses energy 

chiefly by ionizing collisions. "The amount of energy lost by an indivi- 

dual particle over a fixed small fraction of its range is not constant, but 
varies considerably. There are fluctuations in the amount of energy lost at 
each collision, fluctuations in the number of collisions, and a fairly high chance, 
which increases as the number of collisions increases, that a very high energy 
loss, due to a close collision, may occur. 

Landau (1944) has calculated theoretically the probability function W(A, x) 
of a single charged particle such as an electron or meson losing an amount of energy 
Ainasmall distance x. Curve a of Figure 2 shows Landau’s distribution, which 
is approximately gaussian with a high-energy tail. It is given, for convenience, 
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in terms of a parameter (A—A,)/é, and the probability function W(A, x) may be 
obtained from this universal curve at once. Ay is the most probable energy | 
in the particular layer of matter considered, and 


27Ne*p XZ 
f= 8 ot SA’ 
is proportional to the number of electrons in the layer of matter. (Nis Avogadro’ Ss 
number, Z, A, the atomic number and atomic weight of the material, and p its; 
density; m, v are the mass and velocity of the ionizing particle, and e is the? 
electronic charge in E.S.U.) 
[2 
2mv? 


2 - 
Nee (108 +037) where «’= (1 = =) exp (v2/c2). | 


I is an average ionization potential of the atom which Landau takes as JZ, , 


where? =13-5.eve 
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Beam , 7 
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Figure 1. Block diagram of apparatus. Figure 2. Energy-loss distributions for minimum i 
ionizing electrons: curve a, theoretical distri- 
bution according to Landau; curve b, average i 
experimental distribution in argon; curve ¢, , 
average experimental distribution in krypton. , 
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Two assumptions are made in his calculation: (i) the most probable energy 
loss is small compared with the energy of the primary particle; (ii) the observed - 


energy losses are large compared with the average atomic binding energy. 
These large fluctuations are of interest in proportional counter work, since they | 


set a limit to the precision with which the momenta or energy of fast particles may | 


be measured. 
Two interesting points arise from the theoretical calculations when the Landau 
distribution is plotted in terms of energies: for any given material, the relative 


width of the energy loss distribution does not become much narrower as the most | 
probable energy loss is increased, and for a given thickness of the layer of matter, | 
the distribution is narrower in materials of lower Z (even though the most probable | 


energy loss is lower) than in materials of higher ZL. 
These points are illustrated numerically in Table 1, in which are given the | 


theoretical values of the most probable energy loss, and of the widths of the energ@ | | 


distributions for 1-5 Mev. electrons. 
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Table 1 

x (cm. at °% width at 

ene N.T.P.) é (kev.) me beg aaliaah 2 
1-0 0-116 1-14 42-6 
5-0 0-58 6-65 36°6 
Argon 10-0 1-16 14-4 34-6 
100-0 11-6 62:5 29-9 
Hydrogen 10-0 0-143 0-232 25°9 
Fite 10-0 0-143 0-212 28-3 
eon 10-0 0-715 9-4 32-0 
ean Oe 10-0 1-16 14-1 34-6 
Kigpton 10-0 2-32 26-6 36-0 
She 10-0 3-90 44-3 37-4 


§2. EXPERIMENTAL 
The shape of the energy-loss distribution was investigated experimentally by 
firing a narrow beam of mono-energetic electrons from a beta-spectrometer through 
aproportional counter. It is assumed that the ionization produced by the primary 
particle and the secondary electrons in the counter gives a measure of the energy 
‘lost by the primary particle. A block diagram of the apparatus used is given in 
Figure 1. 
: The proportional counter used was a Maze type (Maze 1946) glass counter, 
: 5-5 cm. in diameter, with two thin mica windows, one on either end of a diameter. 
: 


The counter was placed with one window as close as possible to the slit of a beta- 
spectrometer, and an end-window Geiger counter was placed directly in front of 

the other window. The two counters were connected in coincidence, and the 
coincidence pulse opened a gate which allowed the delayed proportional counter 
pulse into the four-channel kicksorter. In this way only those electrons which 
passed right through the counter were recorded. ‘Thus errors due to scattering of 
electrons in the counter were minimized, and background from cosmic rays and 
other sources was practically eliminated. 

The counter was calibrated with K-capture x-rays from radioactive sources 
(Rothwell and West 1950); the calibration lines used were Zn (8-1 kev.), 
Se (10-5kev.) and Pd (20-2kev.). The pulse size distribution from these x-ray 
lines is gaussian in shape. ‘The spread is caused by fluctuations in the number of 
collisions made by the conversion electrons, and fluctuations in the multiplication 
process. This spread becomes very much wider if there are slight impurities in 
the counter; the half-widths at half height of the distributions from the x-ray lines 
were therefore measured at each calibration to test whether the counter had 
deteriorated in any way. 

Two points arise in connection with the mica windows: 

(i) To test the shape of the distribution fairly the energies of the electrons used 
must be greater than 1Mev. This is because the beam of electrons actually 
entering the counter is not quite mono-energetic, owing to fluctuations in the 
amount of energy lost in the mica window (the average energy lost is of the order 
of 100key.). For particles of energies greater than 1 Mev., the specific ionization 
is practically constant, and for fluctuations in energy of the order of those produced 
in the window, all electrons have the same ionizing power. For energies less than 
imev., the specific ionization increases sharply as the energy of the ionizing 
electron decreases, and so the electron beam contains electrons with a considerable 
range in ionizing power. 
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(ii) The additional ionization caused by secondary electrons, formed in the: 
first mica window, entering the counter, is assumed to be equal to the ionization: 
lost in the walls by secondary electrons formed in the counter. S | 


§3. RESULTS 


The results of a series of measurements are recorded in Table 2. ‘The experi-- 
mental values of A, are slightly lower than the theoretical values, except with thet 
1-0 Mev. electrons. In this case the experimental values of Ay are slightly higher! 
than the theoretical. The latter effect is no doubt due to the energy loss of the 
particles in the mica windows, which begins to be important at about 1 Mev.) 
(see column (7) of Table 2). | 

Table 2 | 

(1) (2) (3) (4) (5) (6) (7) 

1 Kr=- CH, 52D 355-2, 2-0 10-2 10:8 

2 ire (Claly 52 oO Soe 2-0 10-4 10:8 

3 Kr+CH, 5D) Se 8 0= 55-2 1°5 10-0 10:6 

4 Kr+CH, 52-2-— 3:-0= 55-2 1-5 9-75 10-6 

5 Kr+CH,  60:0+10-0=70-0 2-0 12°5 13-1 

6 Kr+CH, 60:0+ 10:0=70-0 2:0 12-5, 13-1 

7 Kr+CH, 60:0+10:0=70-0 2:0 12252 13-1 

8 Kr+CH, 60:0+10:0=70-0 1-0 13-4 12:9 13-3 
9 Kr+CH, 60:0+10:0=70-0 pe@ 13-3 12-9 1353 
10 A+CO, 52:0+ 4:8=56°8 1125) 6°3 6:4 
15 A+CO, 56:0+ 6:4=62°4 15 72 dee, 
12 A+CH, 52:0 dQ 5ef5 5:6 5:8 
13 A+CH, 52-0 1) 5°8 5:6 5:8 


(1) run; (2) gas; (3) pressure (cm. Hg); (4) energy of incident electron | 


(Mev.); (5) Ay (experimental); (6) Ay (theoretical); (7) Ay (theoretical, 4 
corrected for mica window). : 
4 
t 


Figures 3 and 4 show the pulse size distributions for the argon and the krypton : 
counters respectively, plotted in terms of the parameter (A—A,)/é. Each set of | 
points is normalized so that the area under the curve is unity, as is that under? 
Landau’s. The ordinate of the resultant curve is then represented by ¢[(A — A,)/€].|| 
The averages of the argon and krypton distributions are drawn in curves 5 and cl 
of Figure 2. / 


0:15 enn a 
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Figure 3. Energy-loss distribution for Figure 4. Energy-loss distribution for minimunf|) 


minimum ionizing electrons in argon. ionizing electrons in krypton. 
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There is a considerable discrepancy between the experimental curves and the 
theoretical curve on the low-energy side, and the distribution in argon is narrower 
than the distribution in krypton. When plotted in terms of the parameter 
(A—A,)/€ the distributions should, in theory, both lie on the same curve. 

Corrections were made for: 


(a) Fluctuations in the multiplication process. A spread of approximately 


_ the width of the x-ray line corresponding to A was assigned to each point on the 
theoretical curve. The shape of the distribution was hardly altered. 


(b) The average lengthening of path due to multiple scattering. This was 
calculated from a formula due to Rossi and Greisen (1941); it is found to be about 
10%, with the krypton fillings and about 5-5 % with the argon fillings. (The value 
is probably somewhat lower in this experiment, since particles scattered through 


angles greater than 8° are not counted.) If several Landau distributions with 


slightly different values of A, are superimposed, the result is a distribution broadened 
on the low-energy side, but a 15°% variation in A, is needed to account for the 
discrepancy between the theoretical and the mean experimental argon curve, and 


an even greater variation (of the order of 25%) for the krypton curve. 


CONCLUSION 


It seems possible that there may be a genuine discrepancy between the 
theoretical and the experimental curves which cannot be accounted for by defects in 
the apparatus alone. In his calculation, Landau assumes an average value for J, 
the ionization potential of an atomic electron, i.e. all the atomic -electrons are 
treated as though they had the same average ionization potential and the same 
chances of collision; this assumption becomes farther from the truth as the atomic 
number increases. Maybe the energy-loss distribution in argon fits Landau’s 
curve slightly better than in krypton because the electrons in the argon counter 
behave more like those in Landau’s theoretical picture than the electrons in the 
krypton counter, as well as because the scattering is less in argon than in krypton. 
It follows that the distribution from a helium counter should fit Landau’s curve 
more closely still. 


ACKNOWLEDGMENTS 
The author would like to thank Dr. B. Pontecorvo for his encouragement and 
many helpful discussions in the earlier stages of the experiment and Mr. D. West 
for his interest and advice. She is grateful to Dr. P. Cavanagh and Mr. W. H. 
Taylor for the use of the beta-spectrometer and associated electronic equipment. 
This work is published by permission of the Director, Atomic Energy Research 
Establishment, Harwell. 


REFERENCES 
LANDAU, L. \1044, J. Phys.; U.S.S.R., 8, 201. 
Maze, R., 1946, 7. Phys. Radium, 6, 164. 
Rossi, B., and GREISEN, K., 1941, Rev. Mod. Phys., 13, 241. 
ROTHWELL, P., and West, D., 1950, Proc. Phys. Soc. A, 63, 541. 


3 0-2 


g16 


The Mechanism of the Low-Frequency Electrodeless 
Discharge in Chlorine and the Influence of Irradiation 


By W. L. HARRIES anp A. von ENGEL 
Clarendon Laboratory, Oxford 


MS. received 2nd February 1951, and in amended form 13th Fune 1951 


ABSTRACT. The current flowing through a 50 c/s. electrodeless discharge in a short 
cylindrical glass vessel with plane electrodes at the ends filled with chlorine at 5-50 mm. Hg 
was investigated with an oscillograph. It consists essentially of between one and fifty or 
more distinct pulses per half cycle, the number increasing with the voltage. At sufficiently - 
large voltage these pulses cover the greater part of a half cycle including the instant of zero | 
voltage. By irradiating the discharge with light from an incandescent lamp the average | 
pulse height decreases, radiation below 4,800 A. being strongly active. By irradiating only | 
the central portion of the vessel the pulse height does not change. Irradiating one wall 
reduces the height of the pulses in every second half cycle, namely those pulses whose | 
electron avalanches start from that wall. | 

Once the discharge starts, electric charges collect on the inner surface of the vessel. Each 
current pulse consists of a series of avalanches which develop between small areas of the : 
inner glass walls. One such area may discharge within a half cycle in several steps. The 
time variation of the applied field, the wall-charge field and the field in the gas are discussed 
and a new mechanism is developed which appears to apply to discharges between solid |) 
dielectrics. 

Atroom temperatures, chlorine molecules form an adsorbed layer, several molecules thick, , 
on the glass walls. Irradiation produces photo-dissociation of Cl,, the atoms being effective » 
in capturing electrons in this layer where slow electrons and many-body collisions are likely. . 
This reduction of the number of secondary electrons reduces the pulse height. Adsorbed 
molecular gas layers also explain observations at higher temperatures. At higher voltages | 
the interval between two successive pulses becomes so small that one pulse reduces the height + 
of the following one. This is probably due to metastable Cl atoms and resonance radiation 1 
diffusing to the wall where they cause dissociation of molecules in the adsorbed layers, the * 
chlorine atoms capturing electrons as before. It also explains why at high voltages the effect 
of irradiation becomes insignificant. 


§1. INTRODUCTION 

‘HENalow-frequency discharge of between 50 and 500c/s. is established in 1 

a glass vessel with external electrodes, filled with a halogen, it has been 1 

found by Joshi and others (Joshi and Deshmukh 1941, Joshi and Deo } 

1943, Deo 1944) that the r.m.s. value of the discharge current decreases when the : 
vessel is irradiated with light from an incandescent lamp. The gas pressures were : 
from 1 to 600 mm. of mercury. 
We have investigated with a double-beam oscillograph the discharge current 
and the applied voltage in chlorine at pressures from 5 to 50 mm. Hg, in cylindrical I 
hard glass vessels of length 1cm. and 3cm. and diameter 6cm., with external | 
electrodes fixed on the ends. The alternating electric fields were up to 6,000 volts : 
per cm. r.m.s. at 50 c/s. parallel to the axis of the vessel. Figure 1 shows a diagram | 
of the circuit. Distortion of the wave form and phase changes were avoided by| 
using suitable carbon resistances. | 
The chlorine was prepared from concentrated hydrochloric acid and potassium|) 
permanganate. After washing with water and drying with calcium chloride, it 
was frozen in a side tube immersed in liquid air. This was then joined to ther 
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apparatus, the whole evacuated and the vessel cleaned by passing a discharge 
through Cl, and re-evacuating, etc. The pressure of the chlorine was controlled 
by altering the temperature of the side tube by cooling it above the liquid air surface. 
The temperature of the tube was measured By a thermocouple and the pressure 
obtained from vapour pressure curves. 
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Figure 1. Circuit for measuring current and applied voltage with a cathode-ray oscillograph. 


§2. GENERAL OBSERVATIONS 

At small voltages a capacitative sine current flowed. When the field in the gas 
exceeded a critical value, a number of pulses or damped oscillations appeared 
every half cycle superimposed on the sine wave. The height of these pulses was 
always large compared with the amplitude of the sine wave. 

As the applied voltage was increased, the number of pulses on each half cycle 
increased, but their average height remained practically unchanged; at the same 
time the pulses started at an earlier phase of the cycle and took up an ever 
increasing part of it. "Though there were continuous fluctuations in individual 
pulse heights, any changes in the average pulse height could easily be observed. 

Figure 2 shows oscillograms of the current for various applied voltages in 
complete darkness and withirradiation. When the vessel was irradiated with light 
from an ordinary 60 watt bulb at about half a metre distance, with the voltage just 
above the starting voltage of the discharge, it was found that the pulses were reduced 
in height (Joshi 1944). At higher voltages the reduction in pulse height became 
smaller and at sufficiently high voltages there was no change in pulse height, 
except for the first pulse in each half cycle. This first pulse was reduced in height 
by irradiation up to the highest voltages applied (about ten times the starting 
voltage). 

We shall first discuss the nature of the pulses and afterwards the effect of light 
on the pulse height. 


§3. NATURE OF THE DISCHARGE PULSES 
(i) Description and Origin of Pulses 


By using a small high-voltage transformer (= 100 watts) critically damped 
pulses were obtained (Figure 3). Their peaks were up to about 200 ja., the time 
from the start to the peak of the pulse was about 2 x 10~* sec. and the total pulse 
time about 2 x 10-*sec. 

In the circuit shown in Figure 1 the glass walls act as condensers which are 
charged by the ionization current flowing through the gas. When the field in the 
gas is sufficiently strong, an initial electron starting at the negative wall will cause 
an electron avalanche in the gas; thus electrons are driven to the positive wall and 
positive ions to the negative one. Positive ions, light and metastable atoms will 
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release further electrons from the negative wall (as in counters with insulated walls) ) 
and a series of avalanches develops. This increase in number of ions and electrons 3 
and the corresponding rise in current would only continue if the field in the gas 3 
remained unchanged. Because unlike charges collect on opposite walls a voltage » 


drop develops across the glass walls. Since the applied potential is sensibly ” 


epics] WITHOUT IRRADIATION WITH IRRADIATION 


(a) 


(b) 


Figure 2. Oscillograms of the discharge current and the applied voltage (r.m.s. values) with 
and without irradiation. : 


~ 2«10amp +}--------- 


¢—-—— 2.10°* see» | 
~ 2x10" sec. t 
Figure 3. ‘The current waveform of a pulse. 


constant during the pulse time, the potential difference across, and the field in, the } 
gas are thus gradually lowered. Consequently the rate of ionization and the }) 
current decrease and approach zero value. (Field distortion through space } 
charges would tend to increase the rate of ionization here.) 
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There is also another cause for the decrease of the field in the gas : when the 
current rises a voltage drop occurs in the unavoidable inductances thereby 
reducing the electrode voltage. This can be seen from the voltage oscillogram 
(Figure 2): each current pulse is accompanied by a kick in the voltage oscillogram. 


(11) Effect of Wall Charges 

Pulses can occur at an instant when the voltage between the external 
electrodes is zero (Figure 2). Since a pulse can only be initiated when the 
field in the gas exceeds a critical value estimated of order 300 v/cm. at p=5 mm. 
Hg, it can then only be due to wall charges produced during the previous half 
cycle if we neglect space-charge fields (cf. point corona in Cl,, Thornton 1939). 

At still higher voltages the current pulses can be seen to occur during almost 
the whole time interval between two successive peaks of the applied voltage 
(B to A of Figure 2(d)). However the polarity of the current pulses remains the 
same although the applied voltage changes sign. This again follows from the 
existence of a wall-charge field; between A and B the field in the gas which is the 
difference between the instantaneous values of the applied and wall-charge fields is 
always in the same direction (see also Figure 4 (5)), and so the current pulses must 
have the same polarity. Since similar results have been found in atmospheric 
air (not given here), we believe that the arguments discussed later in §4 apply 
to other gases and higher pressures than those used above (cf. Gross 1950). 

At high voltage the number of pulses per half cycle is large. Since the average 
pulse height and duration do not change appreciably with voltage, the quantity 
of charge passing through the discharge tube per half cycle is large, and thus 
intense wall-charge fields develop. Their magnitude can be seen from the fact 
that at the instant after the last pulse of a group (B in Figure 2 (d)) the maximum 
applied field is approximately balanced by the wall-charge field. 


§4. POSSIBLE MECHANISMS OF THE DISCHARGE 

(i) The charged inner glass surface of the vessel which lies under the electrodes 
is considered to be an equipotential surface at any instant. Each pulse then 
discharges the whole of these surfaces uniformly but not necessarily completely. 

‘Discharging a surface’ does not necessarily mean that a quantity of the 
adsorbed charges from opposite walls of the vessel are released and travel through 
the gas to the other wall. It is rather the charges produced by electron multipli- 
cation in the gas that move to the walls, first neutralizing the charges already there 
and eventually building up charges of opposite polarity. Hence, during a pulse 
the field in the gas decreases from breakdown to nearly zero value, while the 
current first rises and then ceases. After that the applied field has to rise again 
until the field in the gas reaches breakdown value, when another pulse occurs, and 
so on (Figure 4). 

When the peak of the field in the gas is just above breakdown value (Figure 
4 (a)) only one pulse per half cycle can develop. This is because the field in the gas, 
which is equal to the sum of the external field and the wall charge field, exceeds the 
breakdown value (dotted lines) only once per half cycle. In the steady state the 
wall charge field must have the same magnitude in both half cycles and opposite 
sign. It must be equal to the peak applied field; also Figure 4 shows that 
if it were larger the pulse would occur too soon before the peak, contrary to 
observations. 
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At applied voltages well above breakdown voltage (Figure 4 (4)), because of the + 
large number of pulses, the wall-charge field at the instant before the first pulse of a 
group is large. Since the wall-charge field is then opposite to the applied field, , 
the field in the gas can reach breakdown value while the external field decreases 3 
towards zero. This is why pulses can occur at decreasing parts of the applied 4 
voltage curve. The current pulses do not change direction when the applied | 
voltage changes sign, because the polarity of a current pulse depends on the sign i) 
of the change of charge, which remains the same for each group of pulses. 

When a pulse occurs, the wall-charge field is suddenly reduced by an amount ;) 
which is only a fraction of the initial value (Figure 4(b)). ‘This causes the field 
in the gas to drop to zero at the end of a pulse. Then the applied field, and hence 
the field in the gas, rise again until the latter reaches breakdown value, when | 
another pulse occurs and so on. The wall-charge field thus changes in steps, the ; 
last pulse occurring just before the peak applied voltage. 


ll 
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Field applied 
by electrodes 


SSS Breakdown fields 


—-—- Wall-charge field 
——— Field in the gas 


Figure 4. (a) Electric fields at low applied voltages. This diagram applies either for the whole of © 
the inner glass surface being discharged once per half cycle or for one small area discharging. 
(b) Electric fields at large applied voltages. This diagram applies either for the whole of the »/ 
inner glass surface being discharged several times per half cycle or for one small area discharging |) 
in steps (mechanisms (i) and (iii) respectively). (c) Fields for one elementary area with a pulse “ 
occurring near zero applied field (mechanism (ii)). ! 
The following objections can be raised to this mechanism : 
At low voltage, i.e. when the number of pulses per half cycle is small, the pulses 
occur very close together. From the above mechanism the pulses should be : 
spaced far apart as a large time interval is required near peak voltage for the field | 
in the gas to reach breakdown value again. 4 
For large applied fields (Figure 4(b)) the pulses should occur closer together | 
when the applied voltage passes through zero as the interval between the pulses | 
depends on the rate of change of the applied field. Such a clustering has not been 


observed. 


For the largest applied fields of the order of ten times the starting field, we should | 
expect to have approximately ten pulses per half cycle, whereas three to five times 
this number have always been observed. | 
(ii) A second view is that a pulse discharges only small areas of the charged 
inner glass surface under the electrodes. Thus, contrary to (i), the inner surfaces 
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are here in general non-equipotential. When electron avalanches develop, the 
discharge is confined to a narrow cylinder extending between two opposite areas. 
After a pulse has occurred the field in the gas in front of these areas is reduced, 
preventing further discharge ; for other regions of the vessel the field in the gas is 
effectively unchanged, so other areas can still discharge. As the wall-charge 
density, the nature of the glass surface etc., vary over the inner surfaces, a pulse 
starts from another small area at a different value of the applied field, and so a group 
of pulses develops. Now each pulse originates from a different small area and 
reverses the charge on it from +g to —q or vice versa once in a half cycle. If, 
however, a pulse reduces the field in the gas over adjacent areas, then the following 
_ pulse could only occur some distance away from the first where the field has not 
been reduced. ‘This was observed at low voltages, where the pulses were few; the 
luminous paths of the discharge appeared to originate from small areas of the glass 
and jump about from point to point inside the vessel. 

Figures 4 (a) and (c) show the fields for a single elementary area. If the applied 
field is just above breakdown value (Figure 4(a)) a single pulse will result near its 
peak. ‘The wall-charge fields are assumed constant between pulses, implying 
that recombination in the gas is neglected. 

At large applied fields the pulse is shifted to the leftin the diagram. A pulse at 
an earlier phase is only reconcilable with a larger wall-charge and applied field. 
Figure 4(c) shows this for a pulse occurring near zero applied field. These 
conditions are, however, incompatible; when the wall-charge reverses, the applied 
field available would be much too weak to build up a strong wall charge in the 
reverse direction. So this view cannot explain pulses which occur near zero 
applied field. 

Other objections to this mechanism are that according to Figure 4(a) and (c) 
the pulse carries charges which would increase with the applied field, contrary to 
observation. Also it would be difficult to understand how the two elementary 
areas which provide the first and last pulse of a group (Figure 2(d) A and B) 
should have such vast differences of starting potential. Also large potential 
differences would develop between adjacent areas of opposite polarity and a 
discharge occurs along the glass surface. 

A third mechanism is proposed which is based on both (i) and (ii). 

(iii) Each pulse is associated with a discharge covering a small elementary 
area. In addition each pulse is accompanied by a change of wall charge which may 
be partial only, and need not necessarily result in a change of polarity. 

For low applied fields where there is only one pulse per half cycle, as in Figure 
4(a), this single pulse originates from that elementary area which has the lowest 
starting potential at that time. ‘The next pulse can come from a different area 
because their starting potentials may change (e.g. if adsorbed gas is removed by the 
pulse). 

At large applied fields the pulses would only partially discharge the areas, and 
the wall-charge field would change in steps (Figure 4(5)). This would explain 
why the pulses have the same sign whether they occur before or after zero applied 
field. It also accounts for the constancy of the average height of the pulses. 
Taking the starting potential and the average size of an area as being constant, the 
change of wall charge is proportional to the starting potential. Also for an 
elementary area the number of pulses per half cycle cannot exceed peak applied 
voltage divided by starting voltage. However, the number of pulses observed 
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was very much larger than this, which suggests that a number of elementary areas 
was being discharged. 

Also the wall-charge field for an area will not remain constant between two 
pulses, as shown in Figure 4 (4), because discharges in the neighbourhood as well 


as surface conduction may change its value. Further discharges along the inner — 
surface would occur if the potential difference between two adjacent sites exceeds _ 


the corresponding breakdown value. 


It can be seen that the objections raised to mechanism (i) and (ii) can all be i 


overcome by mechanism (iii). "The same mechanism may be applicable to dis- 
charges in cavities (voids) in insulators (see Austen and Whitehead 1941, Mason 
L957). 
§5. FURTHER EXPERIMENTS AND THEIR DISCUSSION 
(i) Experiments with Split Electrodes 


In order to test whether a discharge takes place between small areas on opposite — 
walls, one of the external electrodes was replaced by two semicircular plates witha | 


narrow earthed electrode between them. The currents to the two electrodes were 
observed with a double-beam cathode-ray oscillograph. At voltages just above 
starting potential, a considerable number of the pulses to the two halves occurred 
at different instances. This means that discharges originate from areas each 
small compared with that of the electrode. 


(ii) Experiments with a Point Electrode 
One of the external electrodes was replaced by a point electrode (wire of 1 mm. 
radius) touching the centre of the glass wall. One pulse per half cycle appeared 
up to voltages corresponding to roughly twice the starting potential, and at four 


times the starting potential the number of pulses per half cycle varied between — 


three and five. Thus the average number of pulses per half cycle seems to be 
proportional to the ratio of the peak applied voltage to the starting voltage. Also 
the time interval between these pulses was much greater than that between pulses 
with plane electrodes. ‘The number and height of the pulses in two consecutive 
half cycles were always the same. 

The following conclusions can be drawn from this experiment: Since the 
field of a point electrode is only strong within a narrow region, discharges can only 


develop from a small area and so a smaller number of pulses per half cycle occurred. | 


One might have expected a larger number of pulses as the field under the point 
electrode is much higher than the uniform field for the same voltage. Thus, 
probably only one pair of small areas was involved in the discharge process; the 
fact that the number of pulses was roughly equal to the number of times the peak 
field is greater than the starting field (see above) supports this view. It also 
explains why the interval between the pulses here is much larger than with a 
uniform field: with only one small area sufficient time had to elapse until the applied 
field changed by an amount which corresponded to the starting potential. 

It should be mentioned that at high voltages a corona discharge developed 
externally in air under the wire electrode producing its characteristic current 
pattern on the screen. ‘The corona current flowing to the plane electrodes was 
negligible at all voltages. 

Although a point—plane gap has a non-symmetrical field, the number and height 
of the pulses in two consecutive half cycles were always about the same. It would 
seem as if the current, the period of oscillation, and hence the charge per pulse, are 
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essentially controlled by the properties of the electric circuit, whereas the discharge 
acts as a switch triggered by the field. Apparently at low pressures and in the 
presence of a glass wall which evened out the field, polarity effects were not 
pronounced or covered by fluctuations. ‘This is in contrast to the breakdown 
with dissimilar (internal) electrodes of molecular gases at atmospheric pressure. 


(ui) The Number of Electron Avalanches in a Pulse 

A pulse appears to consist of a large number of electron avalanches; this can 
be seen by the following arguments. 

The current at the peak of a pulse has been found to be of the order of 10-4 amp. 
Assuming the pulse of Figure 3 (6), its total charge is approximately 10-* coulomb, 
hence the number of electrons is of the order of 1011. Assuming an average 
field in the gas of about 300 v/cm. (Figure 2(a)) and a pressure of 5mm. Hg, the 
value of X/p would be about 60. Taking the corresponding value of « for HCl 
as approximately 3 (unknown for Cl,) (Townsend 1915, von Engel and Steenbeck 
1932) and the length of the vessel as 3cm., the total number of electrons in an 
avalanche produced by one electron would only be about 104. Hence a very large 
number of avalanches seem necessary to transfer the pulse charge assuming a 
small number of initial electrons at the instant of the first avalanche. 

The time from the start of an individual pulse up to its peak has been found to 
be of the order of 10-* sec. (Figure 3), its total] duration of the order of 10~4sec. 
Estimating the mobility of an electron in chlorine at 1 mm. Hg to be of the order of 
10®cm/sec. per volt/cm., then at 5mm. Hg and 300v/cm., the velocity of the 
electrons would be of the order of 108cm/sec. The transit time between the 
glass walls is about 10-*sec., which is small compared with 10~°sec., and thus a 


_ large number of avalanches have occurred by the time the peak current is reached. 


(iv) Approximate Size of the Area discharged by one Pulse 


Each pulse causes a quantity of charge Aq to be carried over to the opposite 
wall. Assume that it forms a uniform layer over an area A. Let the capacity 
between this area and the electrode be C, then the change of field in the gas for a 
vessel of length d will be of the order 
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At the end of the pulse the field in the narrow volume considered would be zero. 
Hence E will be of the order of the breakdown field, i.e. 300 v/cm. at afew mm. Hg. 
From (1) and with the pulse of Figure 3 (¢=10-§ coulomb, d=3cm.), C is about 
10pr. Since the glass wall is about 0-1cm. thick and its dielectric constant 
about 4, A is of order 1cm®. Though this value represents the upper limit, as Ag 
was taken for one of the largest pulses, it is small compared with the area under an 
electrode. 
(v) The Height and Length of the Pulse 

Assume that V,, the starting potential of the discharge, is given and that only 
one discharge exists atatime. Ifthe shape and height of the pulse were controlled 
by the external circuit only, the maximum current flowing through the gas would 
be OTe (TAG trae, NS (2) 
With V,~1,000 volts, L~50 henries (the leakage inductance of the transformer) 
and C~10-"F., we havez,,~5 x 10-4amp. which agrees in order with observations 
(Figure 3(5)); the measured values should be smaller because the voltage drop 
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across the discharge was neglected in (2). When a longer vessel was used, 7, 
was larger because V’, increases with the length of the vessel for p = const. | 

Again the length of the pulse, neglecting damping by the discharge, can be? 
obtained from w~1/(LC)!2=5 x 104 per sec. ; hence the time of the quarter wave : 
is about 3 x 10-° sec. which agrees with the observed time up to the peak of approxi- - 
mately 2x10-5sec. The damping influence of the discharge and the leakage: 
capacity of the transformer would tend to lengthen the period of oscillation. | 


(vi) Fluctuations in Height and Phase of the Current Pulses 


The cathode-ray oscillograph patterns showed that the peak values and the: 
instant of occurrence of the pulses were fluctuating continuously (Joshi 1944). . 
One major cause seems to be the variation of secondary electron emission from the: 
glass wall, produced by positive ions. ‘The emission coefficient y (number of { 
electrons per positive ion) depends, as we know, on the state of the glass surface, , 
on the type of glass used, its preparation, and particularly on the adsorption of { 
gases. y varies from point to point on a wall and, in addition, it varies with time. , 
Also the distribution of charges on the walls will not be uniform, as the pulses: 
occur in narrow channels. ‘These phenomena cause fluctuations in phase as well |, 
as in amplitude. 


(vii) Separation of the Pulses 


Although the pulses fluctuate in height and phase, they have always been 
observed to be separated, i.e. the peaks never merge into one another. The? 
main reason seems to be that the current pulse produces a voltage drop in the leak- - 
age inductance of the transformer, and stray inductance of the circuit (neglecting » 
the resistance), so that the applied field decreases while a pulse develops. This 7 
stops other areas from developing avalanches (see voltage curves, Figure 2). — 
When the pulse ceases, the applied field returns to its original value and a* 
discharge can develop from other areas. 


; 


§6. THE EFFECT OF IRRADIATION ON THE DISCHARGE 
(1) Photo-dissociation of Chlorine 


By irradiating the discharge vessel using incandescent lamps with various: 
filters, we again verified (Joshi and Deshmukh 1941, Joshi and Deo 1943, Deo: 
1944) that light of A<4,800a. caused the reduction in pulse heights. This: 
wavelength corresponds to quanta of energy 2:58ev., which is the energy of | 
dissociation of a chlorine molecule (Bowen 1946). Thus irradiation seems to) 
cause photo-dissociation of the chlorine molecules into two atoms (one of which? 
is in a metastable state, 0-1 ev. above the ground state). 

Chlorine atoms attach electrons readily forming CI when (a) the electron 
velocity is low, (b) the attachment occurs in the presence of a third body (Cl). } 
In this case the probability of attachment is large because the affinity energy; 
(3-5 ev.) can be transformed either into potential energy leading to dissociation 1 
of Cl,, or into kinetic energy of the particles. In the absence of a third body, , 
the affinity energy has to be emitted as a quantum of radiation, which is a rare 
process. | 
The electrons so lost from the avalanches no longer produce electron) 
multiplication and the height of the pulse drops. The negative ions formed | 
have a negligible mobility compared with the electrons, and do not contribute} 
to the current. The probability of electron attachment to the molecule to form 
Cl, is very small indeed (Loeb 1939). | 
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Electrons are also lost by first dissociating a molecule and then becoming 
attached to one of the atoms. However this process goes on whenever a discharge 
occurs, and represents a loss of electrons independent of irradiation. The next 
step is to show numerically the effect of photo-dissociation by an incandescent 
lamp. 

(11) The Rate of Photo-dissoctation 

The tungsten filament of a 40-watt bulb was considered a black body radiator 
at constant temperature 7 =2,500°K. ‘The number of quanta emitted per second 
in an elementary wavelength range is given by Planck’s distribution law. At 
low p only few quanta are absorbed by the gas, the absorption depending on the 
wavelength. Multiplication of these two factors and integration over the 
wavelength range bounded by the band convergence limit of chlorine and by the 
transmission limit of the glass vessel, gives the rate of absorption of quanta, and 
thus the rate of photo-dissociation occurring. 

If E is the energy per unit volume in a wavelength range \ to A+ dX then 
cE, is the radiant flow in erg cm™ sec"! where c is the velocity of light. The 
number of photons 7, emitted per unit area of surface per second in an elementary 


Tange is ny=c,Ejhv=Sac/{e"—1), anes (3) 


mwhere x=hc/ART, h=6:55x10-*" erg sec., R=1-377x10-™ erg deg! and 


2 =7,500°K. 

The number n, of quanta absorbed in passing through dcm. of gas in the 
wavelength range dA, expressed as a fraction of the number n, entering is 

(gn ueade = CSD cd) —~o,Cd,, (6,Cd <1), Ais (4) 

where C is a concentration. «, as a function of A was taken from the 
measurements of von Halban and Siedentopf (1922). By graphic integration 
n,(n,/M),, the number of quanta absorbed per second in a square glass vessel 
of 6cm. side and 3cm. wide (the beam traversing 6cm. of gas) at 5mm. Hg, 
assuming a 40-watt bulb 50cm. away with an emitting surface area of 0-4 cm? 
(emissivity of tungsten =0-43) was found to be 7:1 x 10 quanta per sec. The 
number of chlorine atoms formed per second is twice this in a volume of 
(6 x 6 x 3) x 100 cm’, or the rate of formation is dN/dt =1-4 x 10" atoms/cm? sec. 


(iii) The Equilibrium Concentration of Atoms 
The equilibrium concentration of chlorine atoms can be found from an 
estimated value of the life of a chlorine atom. Bodenstein and Liitkemeyer (1924) 
found that about one collision in one thousand for bromine atoms results in 
recombination, which means as we know now that | collision in every thousand 
is a three-body collision. We assume that this figure applies approximately 
to chlorine. Thus the equilibrium concentration N per cm? of chlorine atoms is 


found from GNidi= Nir AOOe) Eesti Me (5) 


where 7 is the average life of a chlorine atom. At room temperature, since the 
collision interval is //v, but only a fraction N/M collisions occur between Cl and 
Cl, and every thousandth collision leads to recombination, we have 
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where / is the mean free path of chlorine atoms, 4-4 x 10-4cm., v their velocity, 


3:8 x 104cm/sec., and M the number of chlorine molecules present, all at 
5mm. Hg pressure, viz. 2°7 x 10° x 5/760. ‘Thus we have from (6) Nr=2 x 10". 
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From this, the rate of formation given at the end of § 6 (ii), and equation (5) we 
have N~5 x 10"atoms/cm® and r~ 4sec. 
This shows that chlorine atoms do not recombine appreciably in the gas phase. 
Another way of obtaining N and 7 is by using the equation for three-body 


recombination : dNidt=— KMN2, 0 alo oe eee (7) 


where K is the coefficient of recombination, M is the concentration of | 


‘third-bodies’, and N that of Cl atoms. 
Taking the value of K as 0-8 x 10-*? for chlorine (Noyes and Leighton 1947), 


we have dNijdt=—1-4x10-=N2 > 2 ee (8) 


Equating (8) and the above rate of formation, we obtain the equilibrium 
concentration N~ 2 x 1012 atoms/cm’ and 7~ 8 sec. ‘Thus from these two results 
N is of order 10!2 or 101? and r~ 4-8 sec., assuming that the atoms do not strike 
the walls. 

If recombination occurred entirely at the wall (where the rate would be 
exceedingly high, as the wall can act as a third-body), then the atoms could not 
disappear in less time than they would take to diffuse to the wall. This time is 


found from the equation (x)2~(2Dt)¥2 where x2 is the mean square of the 


displacement, D the coefficient of diffusion and ¢ the time during which the atoms 
travel. Taking the coefficient of diffusion of chlorine atoms as 30cm. sec! at 
1mm. Hg, at 5mm.Hg D would be 6. In our case the wall was about 1cm. 
away from the centre of the tube, sot ~0-1 sec. Thus since the time ¢ for diffusion 
to the walls is small compared with the lifetime 7 required for recombination in the 
volume, we deduce that recombination takes place almost entirely on the walls. 


§7. THE EPRFECT OF IRRADIATING DIFFERENT REGIONS OR Vibe 
DISCHARGE 


To investigate this, a narrow parallel beam of light from a 40-watt bulb was 
passed through the discharge (b~ 10mm. Hg) as in Figure 5. No measurable 


Figure 5. Irradiation with narrow beam applied to different parts of the vessel. 


effect on the average pulse height was observed when the central portion B of — 


the vessel was irradiated. When the region A near the wall was irradiated the 
average heights of the pulses in one half cycle only were reduced, and when the 
wall C was illuminated pulses in the other half cycle decreased. Further it was 


found that the illumination of the gas near the wall reduced those pulses which 


occurred in the half cycle in which the external electrode at that wall was negative. 
This means that the height of a pulse was reduced when the region from which 
the electrons started was irradiated. 

As the heights of the pulses in two successive half cycles were unequal, the 
shorter pulses must either have been greater in number or of longer duration 
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because the net total charge transferred across the vessel in two half cycles has to 
be zero. We did not observe an increase in number of the pulses and the 
sensitivity of our apparatus was not great enough to show any increase in their 
duration. 

Light falling on the central portion of the gas dissociates the chlorine 
molecules into atoms. Now the probability of attachment of electrons to atoms 
is high for small and low for high energies. Thus few electrons are lost from 
the avalanche near the centre, because they are too fast. Electrons have small 
velocities only near the wall at which they start, and here the probability of 
attachment is high. ‘The removal of a few electrons here reduces the 
multiplication and hence the peak current considerably. 

Although in the central portion of the tube atoms are formed which diffuse 
to the wall and few recombine in the gas, no reduction in pulse height was 
observed. This suggests that the conditions near the wall are more favourable 
for the attachment of electrons to atoms. This leads to the assumption of a gas 
layer adsorbed on the wall as distinct from the boundary layer consisting of 
ions and excited molecules as suggested before (Joshi 1946). 

The likelihood of electrons becoming attached to excited atoms is small 
(Massey 1950) because of the Auger effect ; the excited ion would be stable if the 
excitation energy were smaller than the affinity, which is rarely the case. 


§8. INFLUENCE OF ADSORBED GAS LAYERS ON THE PULSE 

There is little doubt that gas molecules and atoms can become adsorbed on 
glass walls where they form one or more layers. The atoms will adhere mainly 
by an adsorption process described by Langmuir (Miller 1949) whereby the 
adsorbed atom can be regarded as being held on the glass by a chemical bond. 
The molecules on the other hand are held by van der Waals forces. Thus the 
atoms are much more firmly held, their energy of adsorption being twenty or 
more times larger than the link produced by van der Waals forces. An atom 
coming from the gas which strikes an adsorbed atom will recombine in 
approximately 10-1*sec., i.e. the time of oscillation of a diatomic molecule, and 
the molecule so formed will return to the gas. The likelihood of this process 
taking place is large because the recombining energy balances the energy of 
atomic adsorption, any excess energy being transferred to the wall which acts 
as a third body. : 

Though there are no measurements available about adsorption of chlorine 
on glass at room temperature, one can deduce from the experimental evidence 
with other gases and walls that in equilibrium there will be a sparsely distributed 
population of atoms on the glass with molecules on adjacent sites, and on top of 
the whole several layers of molecules (Miller 1949). In this system of layers 
the density is high; it resembles a phase which is more akin to a liquid than a 
gas. Chlorine atoms from the gas can recombine on the surface of these layers 
without penetrating to the wall surface. 

When the vessel is irradiated, light is strongly absorbed in the adsorbed gas 
layers where the concentration of molecules is high, and hence a large number of 
atoms is formed. An equilibrium process is set up in the adsorbed layers of 
molecules dissociating and atoms recombining, and an equilibrium concentration 
of atoms is established. ‘Thus secondary electrons released from the glass walls 
have a good chance of becoming attached to these atoms, because near the wall 
the electron energy is low and, moreover, the energy of electron affinity can be 
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dissipated easily in such densely populated layers because a ° third body’ is always 
present. 

In addition to external irradiation, light pulses are emitted from the discharge 
producing photo-dissociation. Again electrons of sufficient speed can dissociate 
molecules by collision. ‘The atoms so formed diffuse towards the walls where 
they recombine quickly; the increase in the number of atoms in the adsorbed 
phase is expected to be very small for the reasons given below. Hence there 
will be only an inappreciable change in the number of negative chlorine ions 
and in the average pulse height. 

We have estimated the light produced by the discharge on the basis that 
up to ten times more quanta are emitted per pulse then there are electrons. 
With the largest value of the pulse height and number of pulses per second, we 
found that the number of active quanta produced by the discharge is several 
orders of magnitude smaller than that from external irradiation. 


§9. REDUCTION IN THE PULSE HEIGHT BY THE PRECEDING PULSE 

It has been pointed out earlier that at higher voltages and in the absence of 
irradiation the height of the pulses in a group decreased gradually (Figure 2 (d)). 
It was also found that the first pulse of a group remained large up to the highest 
applied voltages, and was only reduced in height when the vessel was externally 
irradiated. The explanation is that at increasing voltages the interval between 
two successive pulses becomes shorter. It is evident that the basic factor is 
the time interval at which the influence of one pulse on the following one becomes 
apparent. This was estimated (Figure 2) to be of the order of 10-*second. We 
must therefore look for a process which remains active for that time. It could 
not influence the first pulse of a group because the interval between the last 
pulse of one group and the first of the rest was greater than ten times this 
(Figure 2), even at the highest applied voltage. 

The light emitted from the previous pulse cannot be effective, as it would 
reach the wall almost instantaneously and any atoms formed would recombine 
in 10-'* second. Again the effect cannot be due to positive ions: these positive 
ions would be those few from a previous pulse which during 10-4sec. would 
have moved to the region of the next pulse under the influence of lateral electric 
fields which are due to the non-uniform distribution of wall charges; further, 
positive ions might increase the electron emission from the wall, because they 
have sufhcient kinetic energy to penetrate the gas layers: this would make the 
pulse larger and not smaller. 

One mechanism which seems to explain the above facts is associated with 
metastable atoms. Chlorine atoms have not only a metastable level at about 
0-1 ev. above the ground state, but also one at about 8 ev., just below the resonance 
level suggested by H. Kuhn from the term scheme. These metastable atoms are 
produced by electron collision in the gas. ‘They could release further electrons 
from the glass unless a layer of gas was present, adsorbed to the wall. The 
energy of the metastable atoms is used to dissociate molecules in this layer into 
atoms (dissociation energy less than 2:5 ev. which is the dissociation energy of 
the free molecule), to which any secondary electrons released from the wall can 
become attached. As the lifetime of metastable chlorine atoms at these pressures 
is probably of the order of 10-4second or more, they would continue to arrive 
at the wall for this time after any pulse had ceased. Hence dissociation of 
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molecules into atoms could continue for this length of time and thus a reduction 
of any pulse occur sufficiently soon after a previous one. 

A similar argument would also apply to resonance radiation which is emitted 
from excited gas atoms. ‘This radiation diffuses through the gas being emitted 
from one atom, reabsorbed by another and so on, until it reaches the dense 
molecular phase adsorbed on the glass walls in which it dissociates molecules. 

Either of these two processes would satisfy the condition that it increases 
the concentration of atoms on the wall and continues to do so for an interval of 
10sec. or more. 

It can now be seen why at higher voltages external irradiation has practically 
no effect on the pulse height. Since the pulses follow each other very closely, 
a large number of metastable atoms and resonance quanta and a large 
concentration of unexcited atoms is available in the adsorbed layers at the instant 
of the next pulse. As the irradiation reduces the height of the first pulse to about 
that of the next, yet does not further reduce the following ones, one has to conclude 
that with a given light source the concentration of atoms in the adsorbed layer 
cannot exceed a saturation value. Hence at large voltage the heights of the 
pulses remain unchanged when irradiation is applied (apart from the first pulse 
of a group). Other workers have found that the change in r.m.s. current due to 
irradiation became smaller the higher the voltage. 


§10. EFFECT OF TEMPERATURE 
It has been shown (Joshi and Deshmukh 1941, Visvanathan and Kuppuswamy 
1950) that increase of temperature decreases the change of the average current 
due to irradiation. ‘This may be so because an increase of temperature would 
reduce the number of molecules in the adsorbed layer on the wall, on which 
the reduction of pulse height seems to depend. 
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The Saturation Magneto-Resistance of Iron-Aluminium Alloys 


In a recent publication (Parker 1951) the author suggested that the magneto ae 
coefficient at technical saturation magnetization of iron alloys containing up to 9°5 /o silicon ( 
can be considered to depend upon two distinct terms ; these are associated with the | 
temperature-dependent and independent scattering of conduction electrons by the crystal | 


lattice of the alloy. Based upon this hypothesis, the following equation was derived : 
Ap/p=[Aps/pst+ (Apel pr)po(1 + Bt)/psl/[1+po(1+Ft)/ps], «+ + - (1) 


where Ap/p and Ap,/py are the saturation magneto-resistance coefficients of the alloy | 
and of pure iron respectively, Aps/pg is a quantity characteristic of the silicon atoms present 
in the alloy, po(1-+t) the resistivity of iron at ¢°c., and pg the temperature-independent | 
contribution to the resistivity of the silicon solid solution. It was further postulated that | 


Apy/py and Aps/pg were proportional to the square of the spontaneous magnetization 
of the alloy. Values for (Ap/p),, the longitudinal saturation magneto-resistance coefficient 
(magnetization parallel to the direction of current flow) calculated from equation (1), are 


in good agreement with experimental results reported by Shirakawa (1939) if the value | 


of (Apg/pg), at low temperatures is chosen as —1:34 x 10. 


Measurements of resistivity and saturation magneto-resistance have now been made 


on two samples of iron containing 2:68°% and 3-60°% aluminium respectively, in the 


temperature range from —78°c. to +300°c., a range in which a comparison with | 


equation (1) can most readily be made. ‘The chemical analyses of the samples were 
respectively : Al 2:68%, Si0:06%, S 0-:027%, P 0:03%, C 0:07%, Mn 0-08%, Ni 0-085%, 


Cu 0-18%, Fe balance, and Al 360%, Si 0-15%, P 0:028%, C 0-013%, Mn 0-109a5 i 
Ni 0-08%, Cu 0:15%, Fe balance. The measured values of the resistivities of the samples — 


are in good agreement with those already reported for iron—aluminium alloys by Yensen 


and Gatward (1917) and Masumoto (1936). ‘The new values of (Ap/p),, are recorded in i 


Figure 1 together with the curves calculated from equation (1) in the same manner as 


100 0 0 200 ~~ ~—~—«300 2 oo Se 300 
Temperature (°C.) Temperature (°C) . 
Figure 1. (Ap/p)u for 360% Al iron @, Figure 2. (Ap/p), for 360% Aliron @, 
and for two different specimens of and for 2:-68% Aliron x. 


2:68% Al iron x and A; curves, 
calculated from equation (1). 


was previously done for iron-silicon alloys. ‘The value of (Apg/pg),, at low temperatures, 
calculated from (Ap/p), at 10° c., is —0-41x10-*. Agreement between the experimental 


and calculated values of (Ap/p), for both samples is about the same as that reported for ' 
iron-silicon alloys. It is well within the limits of accuracy expected, due to the uncertainty © 
in choosing the correct value of (Ap,/p,)y, the contribution to (Ap/p), of the impurities, 
and the difficulty in correcting the experimental results for the variation of resistivity | 


with spontaneous magnetization. 

The transverse saturation magneto-resistance coefficient, (Ap/p), (magnetization 
perpendicular to the direction of current flow), of the samples was measured between 
15°c. and 300°c., the result being recorded in Figure 2. Unfortunately, direct 
comparison with equation (1) is here not possible owing to the lack of data of (App/pr)y 
for pure iron, thus virtually resulting in two unknown parameters in the equation. It is 
almost certain, however (Webster 1927), that the latter lies between —0-3%10-% and 


—1-6x10~* at room temperature ; assuming any value within this range, a calculation || 


analogous to the one for the longitudinal coefficient predicts a variation of (Ap/p), between 


10°c. and 300°c. of less than 0-3 x 10-8 for both samples, It can be seen from Figure 2 } 


that this is the case, 
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It may be concluded from the above measurements that the longitudinal and transverse 
magneto-resistance coefficients at technical saturation of both samples obey equation (1) 
to within the expected limits of accuracy. 

The author’s thanks are due to Professor L. F. Bates and to Dr. R. Street for their 
interest in this work. 


Department of Physics, R. PARKER. 
The University, Nottingham. 
13th July 1951. 
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The Magnetic Properties of a-Ferric Oxide 


In a comprehensive survey, presented at the Grenoble Colloquy last year, of the 
observed magnetic properties of the various forms, both natural and synthetic, of 
thombohedral ferric oxide a-Fe,O3, Chevallier (1951) justifiably contends that the theory 
due to Néel (1949) to account for the weak ferromagnetism often observed in this substance 
eannot hold when the conditions of chemical formation, or subsequent prolonged 
annealing in oxygen at 1,100° c. or above, preclude the existence of any ferrous ions in 
the structure. 

The essence of Néel’s theory is that the ferromagnetism is due to an ‘impurity’ in the 
form of localized thin layers of what is sensibly magnetite, but with a somewhat deformed 
(cubic to rhombohedral) lattice, the ternary axes of both structures being coincident, as 


_are the common planes of close-packed oxygen atoms normal to this common axis. In the 


discussion at Grenoble, as well as in an earlier published note (1950), Snoek made the 
alternative suggestion that the underlying cause might merely be slight lattice displacements, 
leading to local lack of balance of electron spins; such lattice distortions would be expected 
to be in part relieved by prolonged heating at‘high temperature, leading to the reduction 
in observed magnetic moment which Néel attributed to transformation of magnetite by 
oxidation. 

There would appear to be a fairly obvious way of combining these two basic theories, 
although so far no reference to such a compromise has been seen. Suppose that the local 
dislocations suggested by Snoek were such as to lead to localized thin layers, not of 
quasi-magnetite, but of quasi-cubic ferric oxide y-Fe,O;. The unstrained oxygen lattice 
of the latter is almost identical with that of magnetite, but of slightly smaller dimensions 
(4 ~ 8:3. compared with 8:39a.). Hence the distance between oxygen nearest neighbours 
in the (111) cubic plane for y-ferric oxide is at most 2:94 a., compared with 2-97 a. for 
magnetite. The amount of distortion required to make this plane coincide with the 
(0001) rhombohedral plane of a-ferric oxide, in which the oxygen nearest neighbours are 
2:90a. apart, is thus only one half of that needed to bring about the conjunction with 
magnetite envisaged by Néel. 

In this hypothetical picture there are no unwanted ferrous ions, and we have a basic 
close-packed oxygen lattice, the packing being mainly hexagonal but locally changing 
towards cubic in thin transition zones, a type of varying structure which is already well known 
in pure metallic cobalt after suitable treatment (Lipson and Edwards 1942). The saturation 
magnetization per unit volume of y-ferric oxide being ~0°8 times that of magnetite, the 
conditions for ferromagnetism are much the same as in Néel’s theory. Furthermore, the 
temperature at which y-ferric oxide transforms to the « form is = 650°c., so that it is less 
surprising that the ferromagnetic Curie point found for weakly ferromagnetic w-ferric oxide 
should apparently coincide with the antiferromagnetic transition point 675°c. On heating 
to 1,100°c. or above, the reduction of lattice strains suggested by Snoek would lead to the 
partial destruction of the y-ferric oxide ‘islands’; but it is possible that carrying out this 
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annealing in oxygen impels a tendency to acquire surplus oxygen, with some consequent 4 
diffusion of ferric ions to a distribution more nearly of the y type than of the a, and that b 
this may be the most stable arrangement after cooling, so that the ferromagnetism has not H 
been wholly destroyed, as noted by Chevallier. If this is indeed so, a similar building up of [ 
oxygen on a lattice heavily dislocated by grinding may also be a major cause of the marked q 
increase in magnetic properties on heating in oxygen at 1,100°c., also noted by Chevallier + 
for a powder subjected to such treatment. 


Magnetic Tapes Division, W. P. OSMOND. |, 
Minnesota Mining and Manufacturing Company, Ltd., 
Slough, Bucks. 
26th July 1951. 
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An Introduction to Laboratory Technique, by A. J. ANSLEY. Pp. xv +288. 
2nd Edition. (London: Macmillan, 1950.) 16s. 


This is a second edition of a book first published in 1938, and written for the guidance © 
of laboratory stewards and for those teachers of physics who lack trained laboratory assis- - 
tants. The author, an experienced technical laboratory steward at a university college, 
has added new material to this edition, e.g. the properties of some plastics, the use of the 
epidiascope and film projectors, etc., in an endeavour to bring the volume up to date. The 
chapter headings and much of the contents, however, remain the same as in the first edition, — 
and it is a pity that the opportunity was not taken to omit certain parts of the old text to” 
allow the inclusion of more useful information. For instance, the description of the Kelvin © 
method of measuring the resistance of a galvanometer occupies half a page and yet noi’ 
mention is made of any modern universal measuring instrument. 

In many ways the author does not appear to have kept up with the rapid advance in‘) 
laboratory techniques and apparatus of recent years, e.g. copper-oxide rectifiers are men- | 
tioned with no reference to the selenium variety, and one looks in vain for such important, 
items of modern laboratory apparatus as the cathode-ray oscillograph, the power pack,: 
the photocell, gas discharge lamps, the thermocouple, etc. Again, the references to vacuum! 
techniques are rather scanty, and although no detailed information can be expected in ai) 
book of this size, modern developments in laboratory techniques such as electrolytic polishing 
might at least be recorded, with the sources of further information. 

Several useful tables are included, although the author appears to have missed again a‘ 
point raised with regard to the first edition, viz. the absence of units as in Table V, and nov 
indication is given of the temperature scale in quoting linear coefficient of expansion ini 
Table XX. The book is well produced, but would be the better for the omission of certainr 
illustrations which have little instructive value, e.g. Figure 40. 


R. W. B. 8.3) 


Einfihrung in das Studium der Physik, by W. FINKELNBURG. Pp. 128. (Heidel-| 
berg: Carl Winter, 1950.) DM 4.50; bound in cloth DM 6.30. | 


This interesting booklet is meant to be a guide for prospective German physicists.) 
It contains a little of everything : the aims and problems in physics, the history of physical 
thought, the methods of research, the University syllabus, the main sections in physics and), 
the type of person likely to be an ideal physicist. There is a chapte: on the physicist in! 
various professions which is, alas, the weakest chapter of all. It contains contributions 4 


several sub-authors about the physicist as academic and school teacher, industrial researcher 
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employee in the electrical and chemical industries and about women physicists. A chapter 
is included on physical literature listing textbooks, handbooks and periodicals, and giving an 
introduction into the international decimal classification. Finally there are some notes about 
the German learned societies and professional bodies. 

In Germany physicists are trained at universities and technical academies. During the 
first two years (—4 semesters) 14 to 18 hours per week are spent on lectures on compulsory 
subjects and about three hours on optional subjects apart from the time for practical work 
and seminars. At the end of that period a preliminary diploma examination is taken. During 
the following two years the number of lectures is reduced and the time for practical and 
seminar work increased. At the end of this period the diploma examination (or teaching 
examination) is taken; this requires a miniature thesis but not original research. A further 
two years spent on original research are necessary to complete a dissertation and after an 
examination in the form of our viva, the degree of Dr.Phil. or Dr.Ing. is awarded 

One cannot help feeling that in those chapters in which the author ‘ advises ’ his future 
colleagues, the dogmatic element becomes sometimes so strongly pronounced that it makes 
rather amusing reading; for example he says, ‘“‘ anyone who wants to stick to fixed working 
hours and wants to be left alone in the evenings or spend his spare time on a hobby, should 
abandon the idea of ever becoming a good physicist”. The author’s only support for this 
opinion is that of a well-known mathematician who once said that ‘‘ modern physics is 
actually much too difficult for physicists ”’. 

May I suggest that the author mentions in a future edition the value of personal contact 
by attending meetings at home and abroad, of a knowledge of foreign languages, and last 
but not least, the need of a revival of the travelling scholar. A. VON ENGEL. 


The Industrial Applications of Gas-filled Triodes (Thyratrons), by R. C. WALKER. 
Pp. ix+325. (London: Chapman and Hall, 1950.) 40s. 


This book is an introduction to the many applications of gas triodes and similar valves. 
It commences with a chapter giving a clear and detailed account of the fundamental 
characteristics of the hot-cathode gas-filled valve without entering the realm of fundamental 
theoretical physics. The practical aspects such as cathode temperature, gas pressure, 
ionization and de-ionization times, grid current, etc., are explained so that the reader will 
appreciate the technical problems when they are discussed in the later portions of the book. 
Following this first chapter is a brief chapter on the gas tetrode explaining the improved 
characteristics that have been obtained compared with the triode. Here the author could 
with advantage have given some details of the type 2D21 valve which has become, in recent 
years, a very popular valve for small electronic devices. ‘This valve, although illustrated in 
a manner which does not indicate its small size, is not mentionéd in the text. 

The third chapter deals with the basic principles of control circuits in which A.c., D.c. 
and pulses are applied to the grid combined with a.c. and D.c. on the anode, while the next 
chapter deals with typical examples of simple switching circuits and their applications to a 
wide range of devices such as oscilloscopes, counting circuits, thermostats, printing 
machinery, and control mechanisms employing a photocell to provide the controlling 
potential. 

The fifth chapter deals with heavier current devices such as are used for controlling the 
heating time in electric welding and for the automatic stabilization of the frequency of rotary 
generators, while the following chapter discusses the application of gas valves to the control 
of current and voltage in generators and motors. This is followed by an account of com- 
mutating devices such as current rectifiers and D.C. to A.C, converters. 

The concluding chapter contains details of other types of grid-controlled valves such as 
the ignitron, the cold-cathode triode and the strobotron. Each chapter has its own refer- 
ences and in addition there are comprehensive references to associated subjects not 
mentioned in the text. 

The book is precisely written and at the same time is easy to read, while the diagrams are 
clear and well explained. It should prove to be a valuable addition to the bookshelves of all 
physicists and engineers who are not actively engaged in work on the subject, but who 


encounter various problems of thyratron control in the course of their researches. 
G. E, ASHWELL. 
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Experiments on the Production of Very Low Temperatures by Two-Stage Demag- 
netization, by J. Darsy, J. Harron, B. V. Ro.uin, E. F. W. Seymour and 
H. B. SILsBeEe. 


ABSTRACT. A two-stage demagnetization apparatus employing a superconducting wite 
as the thermal switch between the stages is described. Temperatures in the region of 
0-003° kK. have been obtained using a magnetic field of only 4,200 gauss. The problem of 
thermal insulation was investigated and the residual heat inflow after demagnetization was 
reduced to oneerg/min. The possibility of employing a superconducting switch in cyclical 
magnetic cooling systems and in nuclear demagnetization is discussed. 


A Note on Frohlich’s Theory of Superconductivity, by KuN Huane. 


ABSTRACT. It is shown that the physical mechanisms underlying the two energy terms 
E, and E, in Froéhlich’s theory are very different: EF, is largely the energy change due to a 
modification of the lattice vibration frequencies caused by the adiabatic deformation of the 
electrons; E, is on the other hand a dynamical term, representing the energy change due to 
virtual collisions between lattice oscillations and electrons. A proof is given that certain 
inadmissible consequences of the theory, pointed out recently by Wentzel, derive from the 
energy term £,; the superconductive behaviour discussed by Frohlich arises however 
entirely from the term F,. Although both terms follow from the same perturbation treat- 
ment, the analysis given makes it seem not unlikely that whatever mechanism actually 
inhibits FE, will leave the superconductive behaviour deduced from EF, not substantially 
affected. 


Scintillations from Organic Crystals: Specific Fluorescence and Relative Response 
to Different Radiations, by J. B. Birks. 
ABSTRACT. The scintillation response S of organic crystals depends on the nature and 


energy E of the incident ionizing particle, of residual range r. The specific fluorescence 
dS/dr is not in general proportional to the specific energy loss dE/dr. By considering the 


~ quenching effect of the molecules damaged by the particle on the ‘ excitons ’ produced by it, 


it is shown that dS/dr=(A dE/dr)/(1+-kB dE/dr). Aand kB are constants, which have been 
evaluated for anthracene from observations of S and FE, and the range—-energy data. Curves 
are computed for the relative response |S of anthracene to electrons, protons, deuterons and 
a-particles of E up to 15 Mev., and these are shown to agree closely with the available 
experimental results. The method used for evaluating the relative response is applicable 
to ionizing particles of any nature or energy, and also to the other organic scintillation 


crystals. 


Effect of the Change in Lattice Parameter on the Width of the Forbidden Energy 
Zone according to Kronig and Penney’s One-Dimensional Lattice Model, by 
E. BILuic. 


ABSTRACT. ‘The effect of changing the lattice parameter in Kronig and Penney’s linear 
lattice model is investigated theoretically. Méglich and Rompe have predicted a decrease 
in the width of the forbidden energy zone as the lattice parameter is increased. It is shown 
that this is true only if the width of the potential troughs is increased whilst the thickness of 
the barriers is maintained constant. However, the opposite effect, i.e. an increase in the 
width of the forbidden zone, is expected when the potential troughs are kept constant and 
the width of the barriers is increased instead. In the intermediate case of a uniform expan- 
sion of potential hills and troughs, a small net decrease in the width of the forbidden zone is 


expected, 
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Influence of a Black Sphere and of a Black Cylinder upon the Neutron Density in an 
Infinite Non-Capturing Medium, by B. Davison. $ 


ABSTRACT. For a black body inserted into a homogeneous scattering medium, the 
neutron density at large distances from the body is determined using the perturbation 

method for the case when the scattering in the scattering medium is assumed elastic and 
isotropic capture is absent and the following four types of the black body are consid am 
(i) a sphere of radius small compared with the mean free path in the scattering medium, 
(ii) a sphere of radius large compared with the mean free path; (iii) an infinitely long cylinder 
of radius small compared with the mean free path; (iv) an infinitely long cylinder of radius 
large compared with the mean free path. Combining the results for the above limiting 
cases with the results obtained by other methods the distribution is also estimated for a 
sphere and an infinitely long cylinder of intermediate radii. 5 


The Excitation Function for the Production of ’Be by the Bombardment of °C by 
Protons, by J. M. Dickson and T’. C. RANDLE. 


ABSTRACT. The excitation function for the production of *Be from C has been | 
determined by an activation method, using the internal proton beam of the 110-inch 
Harwell cyclotron. The threshold energy has been found to be 32 Mev. and the cross ; 


section at 156 mev. is 11:0+1:3 x 10-?? cm?. 
Z 


Gamma-Rays from the Deuteron Bombardment of Boron and the Proton Bombard- : 
ment of Aluminium, by J. G. RUTHERGLEN, E. R. Rak and R. D. SMITH. 


ABSTRACT. The gamma-ray spectra from the nuclear reactions B+d, "B+d, and ~ 
27Al11+p have been measured with a pair spectrometer. Lines are found from !°B+d at 
8:88+ 0-06 mev., 6°740:15 Mev., 6-4+0-15 mev., and 4:38+0-05 mev.; from "B+d at 
4-444+0:05 mev.; and from */Al+p at 12:12+0:1 mMev., 10:-46+0-07 Mev., and 
7:62+0:1 Mev. The correlation of these gamma-rays with the energy levels of the nuclei | 
involved is discussed. 


The Density of the Penetrating Particles in Extensive Cosmic-Ray Air-Showers, 
by C. B. A. McCusker and D. D. MiLtar. 


ABSTRACT. The ratio of penetrating particles to electrons in extensive cosmic-ray 
showers has been investigated experimentally at sea level. For the great majority of showers | 
this ratio, in agreement with previous results, is found to be 25%. A small number of * 
showers was detected in which the proportion of penetrating particles was much greater. 


The Diamagnetic Anisotropy of Large Aromatic Systems : II1I—Structures with \ 
Hexagonal Symmetry, by R. McWerny. | 


ABSTRACT. The molecular orbital method is employed in calculations of the - 
diamagnetic anisotropies of a series of highly condensed aromatic systems having hexagonal | 
symmetry and containing up to fifty-four carbon atoms. A new method of calculation | 
is introduced in which the z-electron susceptibility of the whole molecule is determined | 
by an integration, no knowledge of the unperturbed energy levels being necessary : the } 
inclusion of overlap between adjacent atomic orbitals is shown generally to have little } 
effect on the calculations. 

The results are compared with those calculated using the method of Pauling ; this } 
semi-classical method is found to be quite misleading when applied to large molecules }) 
or to molecules possessing stray C-C bonds. Qualitative valence bond considerations are / 


advanced in support of these conclusions, and certain general results are tentatively |) 
suggested, 
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